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SECTION  I 


INTRODUCTION  AND  BACKGROUND 

A.  FLUID  PROPERTIES  AND  THE  EHD  LUBRICATION  PROBLEM 

In  the  elastohydrodynamic  (EHD)  lubrication  problem  the  fluid  undergoes 
a complicated  stress/deformation  history:  the  stresses  are  transient,  and 
have  both  isotropic  (pressure)  and  shear  components.  (See,  for  example, 
reviews  by  Dowson  and  Higginson  [1]  and  Tallian  [2] . ) The  resultant  compres- 
sional  and  shearing  deformations  must  depend  on  the  viscoelastic  properties 
of  the  fluid,  but  in  a rather  complicated,  and  unknown  way.  One  coule  pose 
the  question:  Given  a particular  stress  history,  what  are  the  physical  prop- 
erties of  the  fluid  which  must  be  known  in  order  to  determine  the  resulting 
deformation  history  (or  vice  versa)? 

Figure  1 illustrates  the  approach  taken  by  Trachman  [3] . The  pressure 
distribution  (peak  pressure  pQ)  and  geometry  in  the  Hertzian  contact  zone  (of 
width  b)  is  indicated  in  this  figure;  the  fluid  undergoes  a time  varying  pres- 
sure history  as  it  flows  through  this  contact  region.  By  assuming  a simple 
viscoelastic  model  for  the  liquid,  and  a specific  pressure  dependent  viscosity 
function,  the  "effective"  or  "time  dependent"  viscosity  for  the  liquid  in  the 
contact  zone  is  calculated  and  shown  in  Figure  lb.  By  comparison,  one  sees 
that  the  "effective  viscosity",  hence  shear  traction,  is  much  less  than  by 
assuming  an  equilibrium  fluid  with  no  relaxational  part  to  the  viscosity. 

Thus  while  one  expects  to  enhance  the  viscosity,  the  transient  contribution  is 
less  than  the  steady-state  value,  therefore  partially  compensates.  For  such 
calculations  however,  a more  complete  knowledge  of  the  shear  and  volume  relaxa- 
tional behavior  of  the  liquids  as  well  as  appropriate  time-pressure  equivalence 
principles  are  required.  One  would  expect  that  two  liquids  of  different  chem- 
ical character,  although  of  the  same  atmospheric  pressure  steady-state  viscosity, 
would  exhibit  a dramatically  different  response  in  the  contact  zone.  Thus,  for 
practical  as  well  as  theoretical  importance,  a fuller  understanding  of  the 
viscoelastic  response  is  required. 


B.  VISCOELASTIC  PROPERTIES  IN  THERMODYNAMIC  EQUILIBRIUM- PHENOMENOLOGICAL 
AND  MOLECULAR  DESCRIPTION 

1.  Behavior  of  the  Linear  Viscoelastic  Parameters 
T](u>),  G ' (id) , k'((d),  k"  (cu),  m'(ou),  m"  (cu),  rg, 

In  order  to  phenomenologically  fully  describe  the  mechanical  response  of 
a liquid  in  thermodynamic  equilibrium  at  a pressure  and  temperature  P,  T we  need 
the  shear  viscosity  function  T]s(uu),  the  shear  modulus  function  G^uu),  the  volume 
viscosity  function  (u>)  [or  equivalent  K*  (ou)  since  ouT^  = K."  (u))  ] , the  bulk 
modulus  function  k'(uu).  It  is  sometimes  convenient  to  approximate  these  func- 
tions by  their  limiting  values,  and  refer  to  a shear  relaxation  time  Tg  and 
volume  relaxation  time  Tv(t  = 1/(d)  at  which  a transition  occurs  between  limiting 
values. 

To  develop  a physical  feeling  for  the  useful  parameters  in  describing  ma- 
terial properties,  it  is  convenient  to  collectively  examine  Figures  3 through  6. 
In  Figure  3 (see  Plazek  [4] ) the  specific  volume  (reciprocal  density)  of  the  ma- 
terial is  shown  to  vary  with  temperature;  we  shall  allude  to  analogous  variation 
with  pressure.  We  ignore  the  possibility  of  crystallization.  If  we  are  in  the 
" liquid"  regime  of  this  curve,  and  decrease  the  temperature  slightly,  then  one 
expects  the  volume  to  decrease  to  its  new  equilibrium  value  " instantaneously" ; 
certainly  faster  than  our  11  time  scale  of  measurement”  T,  for  which  we  will  assume 
t = 1 second.  Thus  one  could  say  that  the  structural  relaxation  time  t , is 
« t.  The  same  would  be  true  for  an  increment  in  pressure.  By  contrast,  in 
the  "glassy"  regime  of  this  curve,  in  response  to  a decrement  in  temperature 
(or  increase  in  pressure)  the  volume  would  change  very  slowly.  Thus  the  struc- 
tural relaxation  time  t , the  time  for  the  sample  to  reach  a new  equilibrium 
volume  is  » T;  indeed  could  be  years.  Hence  notice  that  the  glass  tempera- 
ture is  not  unique  but  depends  on  whether  » t or  « T,  namely  depends  on 
our  choice  for  the  time  scale  of  the  measurement. 

Consider  the  example  of  pressure  dependence  of  the  bulk  dynamic  compliance, 
B (B  = 1/K  , where  K is  the  bulk  modulus),  the  curve  of  Marvin  [5],  Figure  5b, 
in  which  the  measurement  time  scale  is  taken  as  T ~ 1/1000  sec:  at  P=-  2000  bars 
(or  at  high  temperatures),  then  « T and  we  obtain  the  liquid-like  compli- 
ance B ; whereas  at  ±2000  bars,  t » t and  then  we  obtain  the  glass-like  com- 
pliance  BOT.  In  these  examples,  the  T of  the  measurement  is  kept  fixed  and  Tv 
is  changed  by  pressure  or  temperature!  Alternatively,  if  the  P and  T are  held 


2 
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constant  as  in  Figure  5a  and  only  T is  changed,  then  one  measures  a liquid-like 
compliance  Bc  at  T » tv,  and  glass-like  B for  T « t . 7his  figure  illus- 
trates that  Tv  can  be  determined  at  each  fixed  T and  P if  t is  variable  over 
a wide  time  (frequency)  range.  We  now  conclude  that  even  the  " liquid- like" 
regime  of  Figure  3 can  behave  "glass-like"  if  t is  short  enough. 

Hence  the  viscoelastic  properties  of  the  lubricant  Be(P),  (P)  and  Tv(P) 
must  be  determined  as  a function  of  pressure.  And  more'. 

Turning  next  to  the  viscoelastic  response  of  a lubricant  to  (dynamic)  shear 

deformation,  Figures  4a  and  b indicate  additional  parameters  required  to  describe 

the  fluid.  There  are  two  characteristic  relaxation  times*  t and  T if  the 

c g 

lubricant  contains  polymers,  through  molecular  models  these  are  associated  with 
conformational  change  of  the  polymers  T , and  orientational  changes  of  the  sol- 
vent molecules  t . Thus,  for  t » T > T (or  u)  « cu  < cu  ) , the  lubricant  be- 
s 7 eg  eg* 

haves  liquid-like  with  a shear  viscosity  and  "zero'  shear  modulus;  for 
Tc  > T > t , the  viscosity  drops  to  a new  level  of  and  the  shear  modulus  g' 
increases  to  rubber-like  values;  while  for  t < the  viscosity  drops  further 
and  G/  reaches  a glass-like  value. 

Hence  the  shear  viscoelastic  properties  of  the  lubricant  7]^(P),  TU(P), 

Ga(P),  Tc  and  must  also  be  determined  as  a function  of  pressure. 

As  an  aside,  it  should  be  mentioned  that  a polymeric  fluid  being  sheared 
will  also  generate  normal  stresses  which  tend  to  push  the  contact  surfaces  apart. 
The  magnitude  of  these  normal  stresses  are  related  to  T|o  and  T . 

2.  The  Free-Volume  Picture:  Concepts  of  Time-Temperature  Superposition 

One  unifying  model  of  liquids  which  has  had  a large  degree  of  success  in 
tying  together  tempersture  and  pressure-dependent  properties  has  been  the  free 
volume  model  of  Doolittle  [6] . The  idea  is  that  the  actual  volume  occupied  by 
the  molecules  is  less  than  the  actual  sample  volume  v,  so  that  the  difference 
is  the  free  volume  v^.  The  larger  the  fractional  free  volume  f = vf/v,  the 
easier  it  is  for  molecules  to  slide  past  each  other  or  to  rearrange  conformations. 


In  every  case,  by  "characteristic  relaxation  time"  we  mean  the  mean  value  of 
the  relaxation  time  distribution  function. 
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Thus  Doolittle  postulated  that  the  temperature  and  pressure  dependence  of  f 
would  be  linear: 


f = f + ot  (T  - T ) + P (P  - P ) 
o f o f o 


(1) 


and,  further:  that  at  two  different  temperatures  T and  T , each  characteristic 


relaxation  time  would  change  in  the  same  ratio  as  the  zero  shear  viscosity, 
then  a temperature  shift  factor  is  defined: 


s Xc/T  _ ^o/T 

aT  T Tl 

c/T0  o/Tq 


(2) 


from  which  he  obtains 


aT  = B(?  - f) 


(3) 


This  equation  was  extended  by  Ferry  [7]  to  the  WLF  equation 


-C.  (T  - T ) 
1 o 


aT  C„  + T - T 


(4) 


which,  in  combination  with  Eq. (2)  has  formed  a time- temperature  superposition 
principle  that  has  been  successful  when  applied  to  conformational  relaxation 
times.  The  inference  of  such  an  equation  is  that  the  character  of  the  relaxa- 
tion mechanism  does  not  change  with  temperature  (pressure),  but  only  the  mag- 
nitude and  characteristic  time  scale  will  change.  This  is  not  true  for  all 
relaxation  mechanisms.  The  utility  of  a time- temperature  superposition  principle 
is  that  predictions  can  be  extrapolated  to  time  (or  temperature)  ranges  other 
than  those  measured. 


3.  Hie  Free-Volume  Picture:  Concepts  of  Time-Pressure  Superposition 


Time-pressure  superposition  is  far  less  developed  for  large  pressure  ranges. 
Ferry  [7]  has  extended  the  concepts  above  to  obtain: 


(B/2.303  f ) (P-  P ) 
i o o 

109  ap  * — t 

o f 


(P  - p ) 
o 


(5) 


The  free-volume  model  would  say  that  only  the  final  volume  is  important  regard- 
less of  the  temperature  or  pressure  history  route.  This  idea  needs  further 
testing.  Other  approaches  have  been  developed  by  Cohen  [8]  and  by  Gibbs  [9] . 
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C.  VISCOELASTIC  RESPONSE  TO  THE  TRANSIENT  (NONEQUILIBRIUM)  EHD  SITUATION 

This  subject  of  the  relation  of  equilibrium  properties  to  transient  response 
has  not  thus  far  been  studied  extensively,  and  therefore  is  part  of  our  experi- 
mental goals.  One  useful  view  taken  by  Harrison  [10]  and  by  Trachman  [3]  uti- 
lizes free  volume  concepts  discussed  above.  The  central  argument  can  be  viewed 
as  if  volume  is  the  unifying  concept,  then  pressure  (and  temperature)  history 
may  be  put  aside,  and  the  response  may  be  determined  by  integration  over  volume 
history,  with  the  assumption  of  equilibrium  properties  at  each  segment  of  the 
history. 

D.  EXPERIMENTAL  TECHNIQUES 

1.  Mechanical  Spectroscopy:  Steady  State,  Low  Frequency,  Ultrasonic 

The  mechanical  relaxation  spectra  of  interest  for  lubricating  liquids  can 
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reach  from  steady  state  (D.C.)  to  10  Hz.  The  part  of  the  spectrum  of  interest 

will  of  course  change  for  each  pressure  (and  temperature).  As  indicated  earlier 

the  functions  which  must  be  known  are  the  shear  viscosity  and  modulus,  and  the 

bulk  viscosity  and  modulus.  While  a variety  of  mechanical  techniques  extending 
0 

from  1 Hz  to  10  Hz  have  been  used  in  this  laboratory  (Figure  2)  for  liquids  at 
atmospheric  pressure,  not  all  of  them  are  directly  applicable  to  high  pressure 
techniques.  Some  exceptions  are  applications  of  similar  methods  to  high  pressure 
problems  by  Barlow  [11],  Marvin  [5]  and  by  Lamb  [12]. 

The  primary  approach  taken  in  our  studies,  is  to  construct  a practical  high 
pressure  cell  optimized  so  that  the  widest  frequency  range  of  measurements  can 
be  made  in  the  same  cell  - both  under  equilibrium  and  transient  pressures.  Be- 
cause of  the  limit  of  small  sample  volume,  we  emphasize  the  techniques  which 
utilize  light  scattering,  and  those  which  couple  mechanical  viscoelastic  tech- 
niques with  light  scattering  - rheo-optic  techniques. 

2.  Light  Scattering  - Twelve  Years  of  Progress 

With  the  advent  of  lasers,  fast  photomultipliers,  highest  quality  Fabry- 
Perot  interferometers,  high  speed  (computerized)  autocorrelators  - material 
property  information  - extracted  by  light  scattering  has  developed  very  far  in- 
deed these  past  twelve  years.  Here  we  will  only  list  a few  techniques  relevant 
to  lubrication,  a few  developed  in  out  laboratory  for  the  past  eight  years,  many 
developed  elsewhere  (see  Fleury  f*13] , Cummins  [14]). 
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The  general  features  of  the  light  scattering  experiment  is  found  in 
Figure  8.  Light  from  a laser  is  focussed  into  a narrow  beam  (circa  0.1  mm) 
at  the  liquid  sample,  and  the  light  scattered  at  any  set  scattering  angle  0 
is  collected  by  a photomultiplier.  The  "liquid"  scattering  cell  may  be  one 
of  the  new  high  pressure  cells  of  Figures  9 or  11,  or  our  older  high  pressure 
cell  [15]  of  Appendix  A.  The  scanning  photometer  [16]  paper  provides  a more 
detailed  description.  Optical  detection  can  be  direct  as  shown  in  Figure  8 
for  autocorrelation  and  cross-correlation  function  analysis  of  scattered  in- 
tensity; or  a Fabry-Perot  Interferometer  (FP)  may  be  inserted  between  liquid 
and  photomultiplier  (PM)  for  Brillouin  peak  and  Rayleigh  Wing  analysis;  filters 
for  fluorescence  studies  and  analyzers  for  depolarized  scattering  are  also 
introduced  between  liquid  and  PM. 

(i)  Rayleigh  line.  Figure  7 shows  a light  scattering  spectrum  from 
Benzene  at  0 = 90°  taken  by  Cummins  [17]  using  a FP.  The  central  line  is 
called  the  Rayleigh  line  because  the  center  frequency  is  the  same  as  that  of 
the  laser.  However,  whereas  the  laser  line  width  is  a few  Hz,  this  Rayleigh 
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line  is  ~ 10  Hz  and  is  due  to  nonpropagating  entropy  fluctuations,  i.e.,  changes 
in  local  ordering  of  the  molecules  or  rotational  ordering  of  the  molecules 
(when  measured  in  HV) . Notice  W refers  to  the  polarized  component  of  the  light, 
HV  to  the  depolarized  component.  Speed  of  molecular  motion  is  related  to  line 
width. 

(ii)  Brillouin  lines  W polarization.  The  lines  on  either  side  of  the 
central  line  in  Figure  7 are  due  to  pressure  fluctuations  inherent  in  all  ma- 
terials [18,19];  this  becomes  equivalent  to  an  ultrasonically  generated  pressure 
disturbance  which  of  course  propagates  as  a wave.  Such  waves  are  spontoneously 
present  at  all  frequencies  (and  wavelengths),  so  that  one  selects  the  acoustic 
wavelength  by  setting  0,  the  scattering  angle,  and  measures  the  frequency  shift 
with  respect  to  the  central  line  to  obtain  the  velocity.  Together  with  the 
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density  of  the  fluid,  this  yields  the  longitudinal  modulus  M,  where  M=  K + — G, 
see  for  example  Dill  [20]  or  McSkimin  [21] . A more  detailed  discussion  of  the 
nuances  in  extracting  K and  G,  as  required  for  lubricant  properties  and  of 
structural  relaxation  [22]  is  omitted  here.  The  width  of  these  Brillouin  peaks 
yields  the  wave  attenuation,  or  dynamic  viscosity  T1  . 
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(iii)  Brillouin  lines  HV  polarization.  If  the  frequency  and  pressure  are 
such  that  the  material  behaves  glass-like,  and  can  sustain  a shear  stress,  then 
a transverse  velocity  and  attenuation  corresponding  to  and  Tjo,  can  also  be 
determined. 

(iv)  Quasi-elastic  scattering.  When  polystyrene  spheres  of  size  < X 
are  added  in  minute  concentration  to  the  liquid,  they  become  the  primary  light 
scattering  source,  and  their  relative  motions  give  rise  to  time  varying  inter- 
ference of  the  light  reaching  the  PM  (no  PP) . Hie  resulting  photocount  fluc- 
tuations are  analyzed  by  the  computer/autocorrelator  as  seen  in  Ref. 16  and 
Appendix  A.  Hie  autocorrelation  function  yields  diffusion  coefficient  of  the 
spheres,  hence  one  obtains  the  viscosity  (low  frequency  limit)  of  the  liquid. 

(v)  Quasi-elastic  scattering.  A new  technique  under  development  is  to 
replace  the  polystyrene  spheres  used  above  with  colloidal  magnets  and  force 
mechanical  oscillation  at  medium  frequencies  (1  Hz  to  10  KHz)  by  an  external 
magnetic  field.  The  motion  is  analyzed  by  a cross-correlation  technique  re- 
lated to  the  one  used  above. 

E.  OBJECTIVES  OF  THIS  RESEARCH  PROGRAM 

In  elastohydrodynamic  (EHD)  lubrication  the  liquid  lubricant  is  subjected 
to  a complex  transient  stress  history  in  which  the  isotropic  pressuie  components 
may  reach  beyond  200,000  psi  in  durations  of  the  order  of  1 millisecond.  One 
may  ask:  How  do  the  viscoelastic  characteristics  of  a fluid  influence  the  stress 
built  up  in  this  lubrication  situation:  Therefore  our  lonq  term  objectives  are: 

(i)  To  measure  the  rheological  properties  of  a variety  of  liquid  lubri- 
cants over  the  full  frequency  ranqe  of  viscoelai.tic  respvinse  from 
1 Hz  to  10^  Hz;  over  a pressure  range  to  300, TOO  psi;  and  over  a 
temperature  range  from  -100°C  to  +300° C.  Know  edge  of  the  visco- 
elastic response  is  important  both  under  equilibrium  and  transient 
pressure  and  temperature  histories,  as  well  as  in  uniform  and  non- 
uniform  pressure/temperature  fields. 

(ii)  Phenomenological  theories  are  to  be  explored  which  can  relate  the 

above  equilibrium  state  viscoelastic  measurements  to  the  stress ^strain 
response  of  a liquid  under  transient  elastohydrodynamic  lubrication 
conditions.  In  other  words,  the  studies  should  lead  to  time- temperature 
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and  time-pressure  superposition  principles  applicable  under 
transient  and  nonuniform  conditions. 


(iii)  Molecular  theories  are  to  be  examined  which  relate  molecular  struc- 
ture to  the  measured  macroscopic  response,  thereby  establishing 
possible  relations  between  lubricant  properties  and  molecular 
structure. 

(iv)  Experimental  procedures  in  this  laboratory  are  based  primarily  on 
rheo-optic  techniques,  that  is,  application  of  stresses  by  external 
mechanical  or  by  spontaneous  means,  and  measurement  of  response  by 
optical  (laser  light  scattering)  means. 

(v)  Selection  of  materials:  After  consultation  with  Mr.  F.C.  Brooks 
and  Dr.  J.F.  Dill  it  was  established  that  the  fluid  PR-143  AC  would 
be  most  suitable  for  initial  studies,  and  Emery  2958  and  Stauffer  704 
would  be  appropriate  follow-on  candidate  fluids.  In  this  discussion 
it  was  established  that  the  most  meaningful  progress  could  be  made 
if  coordination  between  laboratories  was  established  on  the  basis 
that  each  measured  some  identical  fluids . Of  course,  an  under- 
standing of  the  relationship  between  chemical  structure  and  me- 
chanical properties  of  lubricants  remains  a continuing  objective. 
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SECTION  II 


EXPERIMENTAL  PROGRESS  AND  RESULTS 

A.  TOE  LIGHT  SCATTERING  SPECTROMETER 

Our  light  scattering  spectrometer  is  in  the  process  of  being  upgraded  to 
add  the  capability  of : 

1)  Brillouin  line  width  measurements  for  polarized  and  depolarized 
components. 

2)  Scattering  from  magnetic  spheres  in  linear  motion  in  the  liquid. 

3)  Scattering  from  magnetic  spheres  in  oscillatory  motion  in  the 
liquid. 

These  are  additional  to  the  capability  hitherto  available  of : 

1)  Rayleigh  (quasi-elastic)  scattering. 

2)  Scattering  from  suspended  polystyrene  spheres  in  Brownian  motion. 

In  addition,  for  measurement  of  pressure,  we  must  add  the  capability  to 
measure : 

1)  Fluorescence  spectral  shift  from  ruby  crystals. 

1.  The  Light  Scattering  Photometer 

The  restrictions  imposed  by  a small  volume  high  pressure  cell  (see  dis- 
cussion of  cell),  together  with  the  requirement  for  several  experiments  have 
led  us  to  the  design  cf  a multiport  microscope  light  scattering  photometer  to 
replace  our  current  apparatus.  The  microscope  will  have  three  ports  as  seen 
in  Figure  12:  one  for  direct  visual  or  video  scan  observation  of  liquid-to- 
glass/crystalline  transformations  (and  kinetics  thereof)  under  high  pressure; 
one  port  for  quasi-elastic  light  scattering  (Rayleigh  line)  data  analysis  in- 
cluding oscillating  magnetic  spheres  and  for  Brillouin  and  Rayleigh  wing  spec- 
tral analysis;  and  one  for  determining  local  pressure  by  the  fluorescent  spec- 
tral shift  of  ruby  as  measured  with  an  optical  scan  spectrometer. 

For  the  diamond  anvil  high  pressure  cell,  surface  scattering  can  be  a 
serious  problem.  As  will  be  discussed  in  the  next  section,  ON  THE  HIGH  PRESSURE 
DIAMOND  ANVIL  CELL,  we  have  been  able  to  adjust  geometry  to  reduce  the  surface 
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scattering  contribution  by  several  orders  of  magnitude.  In  addition,  the  triple 
pass  feature  of  our  Fabry-Perot  reduces  overlap  between  Rayleigh  and  Brillouin 
spectra. 

2.  The  Fabry-Perot  Interferometer 

A major  system  component  is  the  Fabry-Perot  interferometer  for  Brillouin 
shift  and  line  width  measurements.  After  comparison  of  specifications  of 
Tropel  and  Burleigh  units,  the  Burleigh  Fabry-Perot  Interferometer  Model  RC10 
with  multipass  option  was  selected.  Aside  from  stability,  we  concluded  that 
the  multipass  feature  of  this  unit  will  be  very  helpful  in  the  reduction  of  the 
contribution  of  scattering  from  the  surfaces  of  the  diamond  cell. 

The  unit  has  been  set  up  and  interfaced  to  the  computer.  The  electronics 
for  individual  adjustment  of  the  three  piezoelectric  drivers  and  the  high 
voltage  sweep  circuit  have  been  constructed  and  interfaced  to  the  digital-to- 
analog  drives  from  the  computer  (Figure  10) . The  present  computerized  data 
acquisition  and  analysis  (autocorrelation)  system  as  seen  in  Figure  12,  con- 
tinues to  work  well.  A computer  algorithm  for  tuning  the  interferometer  and 
for  automatic  frequency  drift  correlation  has  been  written,  but  has  yet  to  be 
implemented.  Next  we  discuss  details  of  the  theory  of  the  Fabry-Perot  inter- 
ferometer and  some  preliminary  results  obtained  with  it. 

(a)  Theory  of  Fabry-Perot 

A discussion  of  the  theory  of  the  Fabry-Perot  interferometer  has  been 
included  as  Appendix  C. 

(b)  Experimental  Results  by  Photographic  Method 

Our  first  Brillouin  scattering  results  have  been  obtained  for  several 
liquids  of  interest  at  atmospheric  pressure  and  with  the  photographic  mode  of 
operation  of  the  Fabry-Perot  interferometer.  It  was  deemed  useful  to  survey 
several  liquids  with  this  less  sensitive  and  less  precise  photographic  mode  of 
operation  prior  to  use  of  the  computerized  photoelectric  detection  mode  in 
order  to  obtain  an  overview  of  the  comparative  liquid  characteristics. 


The  following  liquids  were  examined: 

1)  Benzene 

2 ) Ethanol 

3)  Methanol 

4 ) Cyclohexane 

5)  Acetone 

6 ) Toluene 

7)  Lubricant  Oil  ML071-6 

8)  Lubricant  ATL5066 

9 ) Lubricant  ATL506  7 . 

The  results  for  the  first  seven  fluids  are  illustrated  in  Figures  1 through  7. 
Detailed  analyses  on  some  of  these  results  have  been  performed  and  summarized 
in  Table  1.  The  following  cursory  observations  can  be  made:  (1)  Lubricants 
ATL5066  and  5067  are  highly  fluorescent  (yellow)  for  incident  wavelengths 
of  4880A  and  5145A,  so  that  no  Brillouin  spectrum  was  obtained.  The  experiment 
will  be  repeated  after  we  obtain  an  appropriate  narrow  bandpass  filter. 

(2)  Lubricant  MLO-71-6  exhibits  moderately  weak  Brillouin  lines  and  may  give  us 
some  problems  for  high  pressure  measurements.  (3)  In  each  photograph  several 
orders  are  seen.  The  Rayleigh  line  (more  intense)  for  each  order  has  a Stokes 
and  anti-Stokes  line  associated  with  it,  apparent  at  somewhat  smaller  and 
somewhat  larger  radii.  The  triplet  is  repeated  for  several  orders.  (4)  The 
overall  apparatus  stability  is  now  shown  to  be  very  good,  as  evidenced  by  the 
fact  that  some  of  the  exposures  extended  for  one  hour. 

B.  HIGH  PRESSURE  DIAMOND  ANVIL  CELL 
1.  Description 

The  basic  design  [23]  of  the  rheo-optic  light  scattering  high  pressure 
cell  instrument  is  centered  on  microscopic  observation  of  the  fluid  under  high 
pressure.  The  use  of  a tiny  volume  as  in  the  diamond  cell  (see  Figure  11) 
seems  essential  if  one  is  to  achieve  pressures  in  the  300,000  psi  range,  since 
the  transition  from  liquid-like  to  glass-like  behavior  of  the  sample  would 
resist  the  pressure  build-up.  Furthermore,  the  cell  volume  becomes  orders  of 
magnitude  closer  to  the  volume  in  EHD  contacts,  so  that  one  brings  the  labora- 
tory experiment  much  closer  to  the  real  contact  situation.  Hie  process  of 
generating  high  pressures  (and  impulses)  becomes  simplified,  a single  band 
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TABLE  1 

BRILLOUIN  FREQUENCY  SHIFTS  FOR  SEVERAL  LIQUIDS 
AT  22° C USING  PHOTOGRAPHIC  METHOD 


Sample 

Observed  Brillouin  Shift 
at  90°  in  GHz 

Published  Values 

Reference 

Toluene 

6.04  (4880  A) 

6.08  (4880  A) 

Kato  et  al. 

Ethanol 

4.90  (4880  A) 

5.10  (4416  A) 

Gustafsson  et  al. 

Benzene 

6.37  (4880  A) 

6.16  (5145  A) 

6.59  (4880  A) 

Kato  et  al. 

Cyclohexane 

5.60  (4880  A) 

r ^ ™ 

screw  is  sufficient  to  generate  the  desired  pressures.  The  additional  complex- 
ities that  arise  from  the  microscopic  volume  are  connected  with  sample  loading, 
measurement  of  pressure,  and  surface  scattering.  Furthermore,  microscope 
optics  are  required  for  the  small  sample  volume,  and  were  described  above. 

A diamond  anvil  cell  designed  by  Dr.  A.  Van  Valkenberg,  while  at  the 

National  Bureau  of  Standards,  was  ordered  and  delivered  from  High  Pressure 

Diamond  Optics,  incorporated.  Both  optical  and  sample  handling  procedures  are 
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currently  being  developed  for  these  samples  of  volume  circa  10  cc. 

Pressure  and  pressure  distribution  measurement  can  be  determined  by  the 
shift  in  the  fluorescence  spectra  of  ruby  crystals  suspended  in  the  fluid  [24,25] 
and  locally  monitored  by  the  microscope  optics.  Design  is  underway  for  the  fast 
pressure  pulse  triggering  device. 

Fabrication  of  Washers 

As  can  be  seen  in  Figure  1,  the  liquid  sample  is  contained  in  a small 
hole  (shaded  area)  within  the  washer  (or  spacer)  B.  The  hole  diameter  must  be 

sufficiently  smaller  than  the  hexagonal  surfaces  of  the  diamond  anvils  in  order 

to  provide  adequate  interface  contact.  This  interface  overlap  region  provides 
the  normal  stress  o , to  force  the  washer  material  into  a plastic  domain,  and 
provide  the  shearing  traction  a required  for  a liquid  seal.  Thus,  typically 
the  diamond  surfaces  are  circa  1 mm  diameter,  while  the  sample  hole  is  about 
0.5  mm  diameter.  The  initial  thickness  of  the  spacer  depends  on  choice  of  ma- 
terial but  is  approximately  0. 5 mm. 

Hie  procedure  for  fabricating  the  spacers  is  as  follows:  discs  of 

diameter  from  1/4"  to  1"  were  punched  out  of  the  material,  then  placed  into  the 

cell  between  the  anvils.  Sufficient  pressure  was  applied  to  leave  an  impression 
of  the  diamond  surfaces.  With  the  aid  of  a binocular  microscope,  holes  of  1/2  mm 
diameter  could  be  punched  by  hand  at  the  center  of  the  diamond  impression  using 
a dental  tool.  Alternatively,  the  spacers  were  cemented  to  a block  and  a high 
speed  steel  or  carbide  circuit  board  drill  were  used  to  drill  the  1/2  mm  hole. 

Several  materials  were  tested  to  see  if  they  could  withstand  the  high 
pressures.  Liquids  used  for  this  test  were  usually  water,  benzene  or  cyclo- 
hexane. These  were  chosen  because  high  pressure  crystallization  could  readily 
be  observed.  Thus  far  the  most  satisfactory  materials  for  the  spacers  were 


13 


Inconel  and  stainless  steel.  These  could  be  used  repeatedly  with  little  de- 
gradation. Berylium  copper  was  also  found  to  withstand  the  high  pressures, 
but  could  not  withstand  repeated  cycling  without  bursting.  An  aluminum  alloy 
was  tested,  but  could  not  withstand  the  high  pressure. 

Some  nonmetals  were  also  tested.  Plexiglass,  polymethylmethacrylate, 
were  able  to  withstand  pressures  sufficient  to  freeze  cyclohexane  even  though 
crazing  and  cracking  developed  external  to  the  "sandwiched"  region.  Epoxy 
(Eccobond)  spacers  shattered,  and  a more  flexible  Epoxy  was  unable  to  contain 
the  high  pressures.  Mica  spacers  have  not  yet  been  examined.  It  would  be 
useful  to  have  spacers  of  transparent  material  so  that  the  scattered  light 
could  also  be  observed  at  90°. 

2.  Microscopic  Measurements 
Experimental  Results 

Several  fluids  were  examined  in  the  high  pressure  diamond  cell  under 
the  microscope  in  order  to  study  the  morphological  properties  and  to  gain 
visual  experience  in  the  use  of  the  cell.  The  results  for  crystal  transitions 
in  water,  benzene  and  cyclohexane  are  shown  in  the  photographs  of  Figures  2, 

3 and  4.  The  phase  diagram  of  Kennedy  and  Holser  [26],  Figure  5,  for  water 
provided  identification  of  the  pressures  at  which  the  phase  transitions  occur. 
The  lubricant,  a perf luoroalkylpolyether,  Krytox  143AC  (MLO-71-6)  was  examined 
in  this  way;  no  crystallization  was  observed,  rather  the  sample  underwent  a 
transition  to  the  glassy  state. 

3.  Resolution  of  the  Surface  Scattering  Problem 

The  primary  technical  limitation  of  the  diamond  anvil  cell  for  light 
scattering  measurements  has  been  the  severe  surface  scattering  at  the  diamond 
interfaces.  We  have  successfully  resolved  this  problem  by  bringing  the  inci- 
dent light  into  one  of  the  side  facets  of  the  diamond.  The  refracted  laser 
beam  will  then  enter  and  exit  the  sample  area  (the  hole  in  the  washer)  at  points 
on  the  diamonds  that  can  readily  be  spatially  resolved;  hence  these  image  points 
are  readily  blocked  off  by  an  adjustable  aperture.  Full  discussion  of  the  cell 
geometry  will  follow.  Fluorescent  sodium  fluorescein  solutions  have  been 
utilized  to  demonstrate  this  effect  visually  - the  point  surface  scatters  are 
seen  spatially  separated  from  the  sample  fluoroescence. 
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Description  of  Cell  Geometry 


The  geometry  of  the  diamond  can  be  seen  in  Figures  1 and  2,  which  show 
the  cross-sectional  and  top  and  bottom  views,  respectively.  Notice  that  the 
area  of  the  anvil  end  (gasket  end)  of  the  diamonds  are  smaller  than  the  support 
end  so  that  the  pressure  developed  in  the  support  will  not  force  plastic  de- 
formation of  the  stainless  steel  holder. 

Limiting  examples  of  refraction  at  the  diamond  surfaces  are  shown  in 
Figure  1.  The  refractive  index  of  the  lubricant  is  taken  to  be  1.4,  and  the 
indices  of  diamond  and  air  are  2.42  and  1.0,  respectively.  Rays  pass  through 
the  center  of  the  liquid  sample  (center  of  the  gasket).  Thus,  in  principle, 
rays  incident  at  70°  with  respect  to  the  normal  of  the  major  surface,  enter 
the  liquid  at  42° . By  symmetry,  the  scattered  light  observed  through  the 
second  diamond,  would  appear  externally  at  140°  with  respect  to  the  incident 
beam,  and  the  scattering  angle  in  the  liquid  would  be  84° . For  a ray  striking 
a 30°  facet  (see  also  facet  A of  Figure  2),  that  is  at  normal  incidence  to  the 
facet  but  30°  with  respect  to.  The  window  axis,  the  scattered  beam  could  be 
observed  at  60°  with  respect  to  incident,  and  the  scattering  angle  could  then 
be  twice  59°  or  108°  in  the  liquid.  Thus  in  principle,  a range  of  scattering 
angles  is  possible. 

In  addition  to  refraction,  the  diffraction  limitation  of  beam  size  must 
also  be  considered.  The  laser  beam  diameter  is  1.5  mm,  and  assuming  that  we 
use  a 15  cm  focal  length  lens  to  focus  the  beam,  then  the  diameter  of  the  beam 
at  the  geometrical  focal  point  is  given  by  the  Airy  diffraction  equation 

d = 2.44  (6) 

where  d is  beam  diameter  at  the  geometric  focus  - called  the  "waist"  diameter, 

D is  the  laser  beam  diameter,  f is  the  focal  length  and  X is  the  wavelength 
of  light.  Thus  the  beam  diameter  would  be  100  microns  (about  0. 005" ) . This 
is  a small  restriction  on  entrance  (and  exit)  position.  The  length  of  the 
"waist"  of  the  beam  is  given  by 

f2 

i = 4.88  -r  X (7) 

D 


Thus  the  beam  diameter  would  be  2.4  cm,  sufficiently  small  >ud  long  to  make 
it  through  the  diamonds. 

Notice  that  with  this  choice  of  a 15  cm  lens,  the  incident  beam  is  no 
longer  parallel  but  converges  in  a core  of  angle  0.01  radians  (0.57°).  This 
becomes  the  "uncertainty"  in  the  geometrically  defined  scattering  angle. 

In  practice,  the  choice  of  incident  beam  angle  is  limited  by  mechanical 
constraints.  As  seen  in  Figure  3,  the  diamond  must  be  mounted  in  a support 
plate  A,  which  in  turn  is  held  in  alignment  chuck  B.  The  diamond  is  cemented 
into  a small  recess  (not  shown)  in  the  support  plate.  The  diamond  support 
area,  and  the  thickness  of  the  backing  of  the  support  plate  must  be  sufficient 
to  accommodate  a few  thousand  pounds  of  force  without  excessive  deformation  or 
plastic  flow.  Hitherto  we  have  used  a 1 mm  hole  through  this  plate  to  pass 
the  light  beam.  This  has  restricted  the  angle  of  incidence  to  about  10°  with 
respect  to  the  normal. 

In  the  diamond-anvil  cell  modification  that  is  currently  being  machined 
we  will  slit  the  support  plate,  so  that  a rectangular  aperture  the  full  width 
of  the  diamond  surface  will  be  exposed.  We  believe  that  the  area  of  support 
will  remain  adequate.  Thereby  the  range  of  entrance  angle  for  the  incident 

beam  will  be  extended  to  about  40° . The  mechanical  obstacle  then  becomes,  the 

configuration  of  the  alignment  chuck.  This  configuration  can  be  changed,  but 
this  will  be  done  at  a later  time. 

The  alignment  chuck  B,  is  seen  in  cross-section  in  Figure  3,  and  in 

top  view  in  Figure  4.  As  seen  in  these  figures,  there  are  two  screws  which 

protrude  from  the  chuck  B to  act  as  pivots  for  the  support  plate  A.  Two  addi- 
tional screws  adjust  the  angle  of  tilt  about  the  pivot  axis  and  lock  the  plate 
to  the  chuck.  The  tilt  axes  of  the  opposing  anvil  assemblies  are  perpendicular. 
The  two  anvil  assemblies  fit  into  the  holder  part  E of  Figure  5.  Not  shown  in 
this  figure  is  the  pressure  plate,  and  the  spring  and  lever  assembly  by  which 
the  pressure  is  actually  applied. 

The  purpose  of  these  changes  is  to  extend  the  range  of  scattering  angle, 
and  thereby  extend  the  effective  frequency  range  of  the  Brillouin  light  scat- 
tering measurements.  In  addition,  these  changes  will  provide  greater  flexi- 
bility in  separation  of  transmitted  and  scattered  beams. 
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C.  MAGNETO-MECHANICAL  VISCOMETRY 


In  order  to  cover  the  extensive  frequency  range  of  interest  in  viscoelastic 
response,  techniques  must  also  be  developed  for  measurement  of  steady  state  and 
low  frequency  dynamic  shear  viscosity  measurements  of  lubricants  at  high  pres- 
sures. Inasmuch  as  a diamond  anvil  cell  is  required  to  generate  pressures  in 

excess  of  200,000  psi;  the  working  volume  is  restricted  to  a liquid  cylinder 

-4 

of  0.5  mm  diameter  by  0.5  mm  long  - a volume  of  10  cc.  The  following  tech- 
nique is  being  developed  here  to  accommodate  viscosity  measurement  in  this  small 
volume:  a tiny  magnetic  sphere,  diameter  0.1  mm,  is  placed  into  the  sample  fluid. 
An  external  magnetic  field  is  then  applied  to  force  steady  state  or  oscillatory 
motion  of  the  sphere.  The  motion  is  then  measured  by  optical  means  through  the 
diamond  windows.  This  method  is  analogous  to  our  previous  one  in  which  nonmag- 
netic spheres  were  subjected  to  Brownian  forces  and  observed  optically.  The 
primary  difference  is  that  we  will  be  able  to  supply  the  large  forces  required 
as  the  viscosity  becomes  very  high  at  high  pressures.  The  viscosity  TJ,  of  the 
liquid  is  determined  from  the  Stokes  relation 

F = 6n11vr  (8) 

where  the  sphere  is  known,  and  the  force  and  velocity  must  be  measured. 

In  order  to  quantitatively  determine  the  force  F,  that  will  produce  a 
translation,  requires  an  understanding  of  the  interaction  between  the  magnetic 
sphere  and  the  external  field. 

— ♦ 

The  effect  of  an  external  magnetic  field  B on  an  iron  sphere  is  two-fold: 

(1)  The  field  induces  a magnetic  polarization  M into  the  sphere.  The  magnitude 

—4  —4 

of  the  polarization  (or  dipole  moment)  M = M(B)  is  a function  of  B and  increases 

—4 

with  B until  saturation.  (2)  The  force  acting  on  a magnetic  dipole  is  propor- 
tional to  the  spatial  gradient  of  magnetic  field.  Thus  the  force  can  be  written 
as 

F =>  aM(B)  VB  (9a) 

or  considering  only  the  x-direction 

SB  (x) 

Fx  = aM(B)  — (9b) 

This  nonlinear  relationship  is  only  meaningful  provided  the  field  generated  by 
the  sphere  can  be  neglected.  For  example,  the  sphere  should  not  come  close  to 


another  magnetic  material  such  as  the  pole  of  the  main  magnet.  Equation  (2.1) 

—4  — • 

is  not  only  nonlinear  but  M (B)  will  exhibit  magnetic  hysteresis.  To  resolve 

these  considerations  and  still  obtain  the  maximum  force  F , one  should  make  B 

x’ 

large  enough  so  that  M approaches  its  saturation  (maximum)  value,  and  where  it 
is  least  sensitive  to  changes  in  B.  Thus  in  a region  where  the  magnetic  field 
is  described  by 


we  then  approximate 


and  Eq. (9b)  becomes 


dB(x) 

B (x)  = B + — Ax 

o dx 


M(B)  = m(b  +|^  Ax')  M(B  ) 
' o ax  / o 


F (t)  = aM(B  ) 
x o 


dB  (x,  t) 


Under  these  conditions  the  force  would  become  a linear  function  of  the  gradient 


dB  (x,  t) 


; this  gradient  could  then  be  made  a time  varying  function. 


In  order  to  test  the  range  of  validity  of  Eq. (12)  and  obtain  an  estimate 
of  the  forces  that  can  be  developed,  the  experimental  arrangement  shown  in 
Figure  1 was  employed.  A small  iron  sphere  of  diameter  0.1  mm  was  used  as  a 
test  probe  in  a plastic  cell  filled  with  oil.  Motion  was  observed  with  a tele- 
microscope. The  magnetic  field  Bq,  between  the  poles  of  a large  electromagnet, 
could  be  set  to  various  strengths  by  adjusting  the  current  iQ.  The  gradient 

^ could  be  adjusted  by  the  current  i,  through  a small  coil  fixed  near  one  pole 

piece.  The  gradient  between  pole  pieces  was  not  uniform,  and  the  first  objective 
was  to  determine  the  dependence  of  the  gradient  as  a function  of  position.  This 
is  accomplished  by  finding  the  current  through  the  gradient  coil  required  to 
"just  lift"  the  sphere  at  each  position.  Lift-off  corresponds  to  a force  = mg, 
m is  the  mass  of  the  sphere. 

The  test  procedure  worked  well  and  the  experimental  results  are  tabulated 
in  Table  1 and  presented  graphically  in  Figures  2 and  3.  The  procedure  was  to 
position  the  cell  with  the  sphere  a measured  distance  above  the  lower  pole  piece, 
to  set  the  current  through  the  electromagnet  to  the  desired  fixed  values  of 


18 


magnetic  field  Bq.  The  field  strength  is  approximately  500  gauss  per  ampere 
of  current  through  the  electromagnetic.  The  pole  pieces  were  3”  diameter  and 
1"  gap.  Next,  the  current  was  adjusted  until  the  iron  sphere  lifted  off  the 
surface  of  the  cell  as  observed  with  the  telemicroscope.  Scatter  in  the  lift- 
off current  is  primarily  due  to  irregularities  of  the  cell-sphere  interface 
associated  with  this  crude  test  cell. 

One  restriction  imposed  by  our  test  apparatus  was  that  we  were  unable  to 
traverse  the  full  magnet  gap.  Because  of  the  physical  restrictions  introduced 
by  the  cell  and  the  gradient  coil  instead  of  the  full  range  of  0 to  1",  we  were 
only  able  to  traverse  the  range  0.125"  to  0.5";  hence  our  conclusions  apply  only 
to  this  range. 

The  first  conclusion  can  be  drawn  from  Figure  32.  That  is,  the  static 
field  B , should  be  greater  than  14  amperes,  that  is,  7000  gauss,  to  bring  the 
iron  spheres  tested  here  into  a saturated  condition  where  the  linearized  Eq. (12) 
becomes  applicable. 

The  second  conclusion  can  be  drawn  from  Figure  33.  Over  the  range  of 
position  0.12"  to  0.5",  the  gradient  is  not  independent  of  x,  but  has  the 
smallest  variation  at  the  position  0. 5" . 

The  third  conclusion  from  Figure  33  is  that  at  the  same  position  of  minimum 
slope  in  the  gradient,  the  magnitude  of  the  gradient  is  greatest,  that  is,  the 
least  current  is  required  for  lift-off.  Inasmuch  as  about  20  ma  through  the 
gradient  coil  is  required  for  lift-off,  as  seen  in  Figure  33,  and  we  are  able 
to  pass  more  than  300  ma  without  heating,  the  force  that  we  can  generate  is  about 
15  times  gravity. 

We  have  also  explored  a variety  of  colloidal  magnetic  particles  and  tiny 
iron  spheres  for  use  in  the  high  pressure  cell  dynamic  viscosity  technique.  The 
colloidal  suspensions  compatible  with  the  lubricants  to  be  studied  have  shown  a 
tendency  to  aggregate.  The  tiny  iron  spheres  seem  to  work  well  in  the  lubricants 
These  were  available  from  Whittaker  Corporation  in  250  it  and  50  H diameters, 
suitable  for  our  experiments.  They  can  be  individually  handled.  We  are  testing 
a variety  of  magnetic  field  drive  geometries. 
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TABLE  2 


LIFT-OFF  CURRENT  THROUGH  GRADIENT  COIL 


Position 

Inches 

2.5  A 

5.0 

Current  Setting  of  Electromagnet 

A 7.5  A 10  A 12.5  A 15 

17.5 

20 

.125 

- 

- 

- 

- 

83 

69 

72 

61 

.150 

- 

- 

- 

- 

61 

60 

50 

48 

.175 

- 

- 

99 

74 

53 

51 

46 

39 

.200 

- 

- 

79 

72 

48 

43 

38 

32 

.225 

- 

- 

77 

58 

43 

43 

34 

31 

.250 

- 

90 

64 

57 

38 

32 

34 

28 

.275 

- 

83 

60 

44 

35 

31 

24 

25 

.300 

- 

74 

51 

42 

31 

28 

29 

24 

.325 

- 

69 

51 

39 

28 

26 

21 

21 

.350 

- 

68 

43 

34 

27 

25 

19 

19 

.375 

- 

57 

42 

32 

23 

21 

19 

21 

.400 

100 

54 

40 

29 

23 

19 

17 

17 

.425 

99 

51 

39 

25 

21 

19 

16 

17 

.450 

90 

54 

37 

29 

20 

17 

16 

16 

.475 

87 

49 

33 

25 

19 

17 

18 

18 

.500 

84 

47 

33 

25 

18 

19 

20 

18 
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D.  ULTRASONIC  SHEAR  WAVE  APPARATUS 

1.  Description  of  Apparatus  and  Results 

There  are  several  atmospheric  pressure  reference  state  measurements  that 
are  required  for  proper  interpretation  of  higher  pressure  data.  One  of  the  more 
important  is  the  high  frequency  limiting  shear  viscosity  and  modulus  (transverse 
wave)  and  relaxation  time(s).  Steady-state  shear  viscosity  (low  frequency  limit) 
data  is  available  for  the  test  fluids,  however,  the  directly  measured  high  fre- 
quency limiting  shear  modulus  appears  to  be  lacking.  Since  Brillouin  scattering 
from  shear  waves  is  difficult  to  observe,  the  high  frequency  limit  and  T]^ 
must  be  calculated  from  the  longitudinal  modulus  with  some  assumptions.  There- 
fore, we  assembled  auxiliary  ultrasonic  delay  line  apparatus  [27,28]  for  direct 
measurement  of  G/  and  T]*  even  though  it  is  only  useful  at  1 atmospheric  pressure. 
This  apparatus  has  been  calibrated  in  the  frequency  range  2 to  50  MHz  and  over 

O O 

a temperature  range  from  -10  C to  +70  C.  Thus  this  apparatus  serves  two  pur- 
poses: (i)  to  provide  the  high  frequency  shear  moduli  and  relaxation  times 

mentioned  above  for  bulk  lubricants  and  (ii)  to  provide  the  limiting  shear  re- 
laxational  behavior  for  lubricants  with  polymer  additives  (i.e.,  a measure  of 
the  configurational  relaxation  contribution  of  the  polymer  additive). 

This  apparatus  has  been  constructed,  tested,  calibrated,  and  we  are  now 
in  the  process  of  measuring  Krytox  143. 

Calibration  curves  have  been  obtained  for  several  fluids  which  we  believe 
to  be  Newtonian  up  through  the  frequencies  of  interest. 

Preliminary  results  have  also  been  obtained  on  solutions  of  polystyrene  in 
a di-2-ethyl-hexyl  phthlate  fluid.  These  results  are  relevant  to  the  role  of 
polymer  additives,  and  indicate  that  at  high  frequencies  the  viscosity  of  the 
solution  does  not  become  that  of  the  solvent.  Namely,  theory  has  assumed  that 
when  conformational  relaxation  of  the  polymer  has  been  "frozen  out",  then  only 
the  viscosity  of  the  solvent  would  remain.  This  picture  is  incomplete.  Exten- 
sive data  and  analysis  of  this  phenomenon,  called  the  problem,  together  with 
apparatus  description  and  calibration  are  detailed  in  Appendix  B. 
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2.  Relevance  to  Lubricants 


The  aforementioned  results  led  us  to  a very  important  conclusion  as  to 
the  roles  of  the  polymer  additive  in  a lubricant.  Addition  of  the  polymer, 
naturally,  enhances  the  viscosity,  but  the  mechanism  of  enhancement  may  be 
through  the  traditionally  accepted  "backbone  conformation  relaxation",  or 
alternately  may  be  through  interaction  with  the  lubricant  (solvent)  such  as 
to  increase  the  solvent  orientation  relaxation  time,  and  thereby  the  viscosity. 
While  these  two  mechanisms  are  indistinguishable  in  steady  state  viscosity 
measurements,  the  effectiveness  of  the  additive  can  be  dramatically  different 
in  a high-speed  lubrication  situation.  The  apparent  viscosity  enhancement 
due  to  the  polymer  additive  operating  by  the  first  mechanism  will  disappear  at 
high  speeds  (frequency),  whereas  the  contribution  of  the  additive  operative  by 
the  second  mechanism  will  be  maintained  at  high  speeds. 

Clearly  this  represents  a new  criterion  for  characterization  and  evalu- 
ation of  polymer  additives. 

E.  CAPILLARY  TUBE  VISCOMETRY 

We  have  also  measured  the  temperature  dependence  of  viscosity  and  density 

of  the  test  sample  MLO-71-6  at  atmospheric  pressure,  in  contrast  to  the  pressure 

dependent  measurements  of  Cameron  [29]  and  Winer  [30,31].  The  viscosity  data 

has  been  plotted  as  log  T)  vs.  —— — , where  the  value  Tq  was  determined  so  as 

o 

to  linearize  the  plot.  The  temperature  Tq  was  found  to  be  -120°C.  On  the  basis 
of  free  volume  concepts,  this  is  the  temperature  at  which  free  volume  disappears, 
and  viscosity  becomes  infinite.  Extrapolation  of  the  curve  to  a "glassy" 
viscosity  of  10  poise,  would  place  the  glass  transition  temperature  at  -80° C. 
These  data  only  extend  down  to  -33°,  the  lower  limit  of  the  bath;  our  revised 
dewar  set-up  is  now  capable  of  temperatures  down  to  -70° C. 
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TABLE  3 


FLUID:  MLO-71-6 

COMPOSITION:  PERFLUOROALKYL  POLYETHER 

PRESSURE : 1 ATMOSPHERE 


Temperature  °C 

Viscosity  (Poises) 

Density  g/cc 

25° 

11.96 

1.89 

20° 

16.60 

1.90 

15° 

23.53 

1.90 

10° 

34.47 

1.91 

* 25.0 

11.43 

1.90 

* 37.8 

5.33 

1.88 

* 98.9 

.488 

1.78 

* 148.9 

.160 

1.69 

This  data  has  been  supplied  by  Mr.  F.C.  Brooks  of  the 
Air  Force  Materials  Laboratory,  but  not  on  the  identical 
sample . 


TABLE  4 

FLUID:  ATL-5066 

(Emery  2958) 

COMPOSITION:  DI-2-ETHYLHEXYL  AZELATE 

PRESSURE : 1 ATMOSPHERE 


Temperature  C 


Viscosity  (Poises) 


Density  g/cc 


0.908 

0.911 

0.915 

0.918 


SECTION  III 


CONCLUSIONS 

1.  Separately,  the  technical  problems  concerning  the  diamond  anvil  cell  and 
the  Fabry- Perot  optics  have  been  solved  and  tested;  we  are  now  prepared 
to  couple  the  two.  But  we  still  feel  that  the  development  of  the  Diamond 
Light  Scattering  Apparatus  will  provide  an  optimum  and  unique  facility 
whereby  extensive  data  about  the  phenomenological  and  molecular  properties 
of  the  lubricant  can  be  obtained  under  laboratory  conditions  that  come 
closest  to  stimulation  of  actual  EHD  bearing  lubrication. 

2.  The  utilization  of  colloidal  spheres  for  measurement  of  viscosity  at  high 
pressures  by  technique  of  quasi-elastic  light  scattering  has  been  demon- 
strated as  in  the  revised  manuscript  of  Appendix  A. 

3.  An  initial  study  for  utilization  of  single  magnetic  spheres  for  steady 
state  and  dynamic  measurements  has  been  completed. 

4.  A high  frequency  ultrasonic  shear  apparatus  has  been  completed,  tested  and 
utilized  to  study  the  role  of  additives. 

5.  A new  criterion  for  characterizing  and  understanding  the  role  of  additives 
to  high  speed  lubricants  has  been  proposed  based  on  the  above  results. 

6.  Steady-state  reference  data  on  the  test  fluids  are  reported. 


FLUID  BEHAVIOR  IN 
ELASTOHYDRODYNAMIC  FLOWS 


Figure  1(a)  Stress/Deformatlon  History  In  the  EHD  Problem 
(Contact  Zone) 
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Figure  1(b)  Stress/Deformation  History  In  the  EHD  Problem 
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Figure  4(a)  The  Dynamic  Shear  Modulus  at  Thermodynamic  Equilibrium 
and  Atmospheric  Pressure 
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Figure  4(b)  The  Dynamic  Shear  Viscosity  at  Thermodynamic 
Equilibrium  and  Atmospheric  Pressure 
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Figure  5(a)  The  Dynamic  Compressibility:  Viscosity 
and  Modulus 


32 


0 


OOO 


2000 


2000  000 

Reduced  pressureyP(bors) 


Figure  5(b)  The  Dynamic  Compressibility:  Viscosity 
and  Modulus 
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Figure  7 Light  Scattering  Spectrum  with  Rayleigh  (Central  and 
Brillouin  Peaks  for  Benzene 


Homodyne  Liaht  Scatterin 


Figure  8 Light  Scattering  Apparatus 
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Figure  9 The  New  High  Preeeure  Cell 
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VISCOELASTICITY  LABORATORY 
DATA  ACQUISITION  SYSTEM 


EXPERIMENTS  DATA  PROCESSOR 
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Figure  10  Viscoelasticity/Light  Scattering  Data  Acquisition  System 
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Figure  12  - Captions 

A - Argon  ion  laser,  B - Single  frequency  etalon,  C - Polarization  rotator, 

D - Spatial  filter,  E - Partial  (5%)  reflector  - for  beam  intensity  monitor 
or  alternately  entrance  point  for  HeNe  (red)  laser  (second  source), 

F - Lens  to  focus  beam  into  "waist"  diameter  of  50  microns,  G - Prism.  Rotate 
180°  to  displace  beam  from  "forward  to  backward  scattering"  illumination, 

H - Diamond  anvil  cell,  I - Cell  clamp  with  cell  orientation  adjustment, 

J - Pressure  adjustment  screw  (only  partly  shown),  K - Pneumatic  impact  relay, 

(not  shown),  L - Primary  goniometer  to  adjust  scattering  angle,  - Mirror 
for  beam  steering  in  forward  scattering  regime,  Mj  - Partially  reflecting  mirror 
for  beam  steering  in  back  scattering  regime,  N - Adjustable  slits  or  apertures, 

0 - Microscope  objective,  P - Three  port  microscope  body  attached  to  Table  L, 

M3  - Rotatable  and  removable  mirror  for  beam  steering  to  R,  S;  R - Eyepiece 
or  Vidicon  port,  S - Lens,  T - Piezoelectrically  driven  Fabry  Perot  interfero- 
meter  of  adjustable  free  spectral  range,  U - Multipass  corner  cube  reflector 
and  an  Iodine  filter  (not  indicated),  V - Photomultiplier,  W - Lens,  X - Double 
grating  spectrometer,  Y - Filter  for  fluorescence. 


Figure  14  Brillouin  Scattering  from  Benzene  at  90°. 
d = 1.545  cms 
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Figure  15  Brillouin  Scattering  by  Photographic  Method  from  Methanol 
at  Scattering  Angle  0 = 90°,  Temperature  22°C,  Laser  Wave- 
length 5145  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 
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Figure  16  Brillouin  Scattering  by  Photographic  Method  from  Ethanol 

at  Scattering  Angle  9 = 90° , Temperature  22°C,  Laser  Wave- 
length 4880  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 
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Pigure  17  Brillouin  Scattering  by  Photographic  Method  from  Cyclohexane 
at  Scattering  Angle  0 = 90° , Temperature  22°C,  Laser  Wave- 
length 4880  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 
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Figure  18  Brillouin  Scattering  by  Photographic  Method  from  Acetone 

at  Scattering  Angle  9 « 90° , Temperature  22°C,  Laser  Wave- 
length 5145  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 


Figure  17  Brillouin  Scattering  by  Photographic  Method  from  Cyclohexane 
at  Scattering  Angle  0 - 90° , Temperature  22°C,  Laser  Wave- 
length 4880  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 
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Figure  18  Brillouin  Scattering  by  Photographic  Method  from  Acetone 

at  Scattering  Angle  0 * 90°,  Temperature  22°C,  Laser  Wave- 
length 5145  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 


Figure  19  Brillouin  Scattering  by  Photographic  Method  from  Toluene 

at  Scattering  Angle  0 ■ 90° , Temperature  22° C,  Laser  Wave- 
length 4880  A.  Free  spectral  range  9.68  GHz.  The  spectral 
order  is  identified  by  symbol  M. 
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Figure  20  Brillouin  Scattering  by  Photographic  Method  from  MLO  71-6  at 
Scattering  Angle  0 =»  90° , Temperature  22°C,  Laser  Wavelength 
5145  K.  Free  spectral  range  9.68  GHz.  The  spectral  order 
is  identified  by  symbol  M. 


Figure  21  Brillouin  Scattering  by  Photographic  Method  from  MLO  71-6  at 
Scattering  Angle  6-90°,  Temperature  22°C,  Laser  Wavelength 
4880  A.  Free  spectral  range  9.68  GHz.  The  spectral  order 
is  identified  by  symbol  M. 
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Figure  22  Microphotographs  of  Water  in  the  Diamond  Anvil  Cell 
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Figure  23  Microphotographs  of  Benzene  in  the  Diamond  Anvil  Cell 
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Figure  24  Microphotographs  of  Cyclohexane  in  the  Diamond  Anvil  Cell 
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Figure  25  Phase  Diagram  of  Water  and  the  Melting  Curve  of  Carbon  Dioxide 
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Figure  29  The  Top  View  (Support  End)  and  Bottom  View  (Anvil  End)  of  the 
Ground  Diamonds.  Symbols  S and  L refer  to  the  smaller  and 
larger  diamonds. 
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ELECTROMAGNET 
POLE  PIECES 


Figure  31  Apparatus  to  Evaluate  Magnetic  Lifting  Force  for  High 
Pressure  Viscoroetry  Measurements 
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Figure  33  Curves  of  Constant  Magnetic  Lifting  Force  F **  mg.  (Parameter 
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A NOVEL  METHOD  FOR  VISCOSITY  MEASUREMENTS  AT  HIGH  PRESSURES 


The  technique  of  quasi-elastic  light  scattering  is  utilized  as  a method 
to  determine  the  viscosity  of  fluids  subjected  to  large  hydrostatic  pressures. 
In  this  method,  a small  concentration  of  insoluble  polystyrene  spheres  of  known 
diameter,  0. 109  microns,  were  added  to  the  fluid.  The  high  pressure  cell  in 
which  the  fluid  is  contained  has  three  optical  windows.  The  laser  light  scat- 
tered from  the  suspended  spheres  is  analyzed  to  determine  the  diffusion  coeffi- 
cient of  the  spheres.  Since,  by  Stokes  law, 


kT 

8rrrD 


(1) 


the  viscosity  of  the  fluid  is  determined  from  the  measured  diffusion  coefficient 
and  known  particle  radius  r.  A small  correction  is  required  for  the  contraction 
of  the  sphere  and  for  refractive  index  change  of  the  water  under  pressure. 
Results  on  the  viscosity  of  water  at  25°c  up  to  pressures  of  2.5  kilobars  is 
reported  and  compared  to  published  results. 

The  high  pressure  system  and  scattering  cell  (1)  is  illustrated  in 

•ft- 

Figure  1.  The  system  is  assembled  from  standard  high  pressure  components. 

The  scattering  cell  consists  of  a high  pressure  cross;  three  blank  conical 
plugs  are  drilled  out  and  fitted  with  sapphire  cylinder  windows  which  are 
epoxied  in  place.  Windows  from  1/8"  to  1/4"  diameter  by  3/4"  long  have  been 
used  successfully  to  pressures  of  35,000  psi.  The  fluid  is  filled  through  the 
reservoir  to  overflow  at  the  window  to  ensure  that  no  air  is  trapped,  then  the 
window  plug  is  tightened.  Pressure  is  generated  by  hand,  by  compression 


American  Instrument  Company,  Silver  Springs,  Maryland. 
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of  fluid  in  the  valves.  The  sequence  is:  close  all  valves,  open  valves  A, 

B.  close  A,  open  C,  close  B,  C,  then  repeat  sequence  from  open  valves  A,  B. 
Three  repetitions  of  the  sequence  will  bring  the  pressure  to  35, 000  psi  for 
the  case  of  a water  sample. 

The  light  scattered  by  the  polystyrene  spheres  is  observed  at  a scattering 
angle  of  90°  to  the  main  beam.  The  intensity  of  scattered  light  is  monitored 
by  a photomultiplier  which  is  operated  in  a photon  counting  mode,  thus  the 
number  of  photoelectric  events  (photocounts)  per  second  n (t)  is  proportional 
to  the  intensity  of  light.  The  intensity  is  not  constant  on  short  time  scales 
but  exhibits  fluctuation  associated  with  the  Brownian  motion  of  the  scattering 
particles.  Hie  sequence  of  "instantaneous"  intensities  are  stored  in  real 
time  in  the  memory  of  a computer,  then  the  autocorrelation  function  is 
calculated. 

The  photocount  rate  autocorrelation  function  c(t)  is  defined  in  terms 
of  the  sequence  n(t);  and  in  turn  c(t)  is  related  (2)  to  the  translational 
diffusion  coefficient  D: 

2 — 2DK^t 

c (t ) = <n(t)n(t  + T))t  = <n(t)>  6(t)  + <n  (t)  + Ae  . (2) 

The  brackets  ( ) denote  the  average  overall  starting  times  t;  6 (t ) =1  for 
T = 0 and  6 (r ) =0  otherwise. 

For  our  computer  method,  the  real  time  t and  delay  time  t are  always  a 
multiple  of  T,  that  is  t - jT  with  j = 0, 1, 2 . . . . ; A is  a coefficient  related 
to  the  intensity  of  light  scattered  from  the  macromolecules;  and  the  magnitude 
of  the  scattering  vector  K is  obtained  from 


L 


where  nQ  and  X.^  are  the  refractive  index  of  the  solution  and  wavelength  of  the 

light  in  vacuum.  In  this  experiment,  the  scattering  angle  0 was  kept  at  90° . 

The  correlation  time  — , and  therefore  D,  is  calculated  from  the  slope  of 

2DK 

2 

jtn[c(T)  - (n(t))  ] vs.  t,  where  the  shot  noise,  the  first  term  of  Eq.  (2)  may 
be  neglected  since  it  is  zero  except  at  T = 0,  and  (n(t))t  is  the  average 
count  rate. 

Hater  was  taken  as  the  first  sample  to  test  the  technique.  An  example 
of  the  autocorrelation  function  is  seen  in  Figure  2,  taken  at  a pressure  of 
15,000  psi.  The  diffusion  coefficient  determined  from  the  slope  of  this  curve 
in  conjunction  with  Eq. (1)  yields  the  viscosity  at  15,000  psi.  Our  results 
for  the  limited  pressure  range  of  0 to  35,000  psi  (2.5  kilobars)  taken  at  a 
temperature  of  23°C  are  shown  by  the  points  in  Fig. 3.  These  points  are 
superposed  on  the  extensive  data  of  Butt  and  Capi  (3)  as  represented  by  the 
solid  lines.  The  results  are  only  intended  as  an  illustrative  example  of 
the  method. 

This  method  is  applicable  to  lubricants  where  the  viscosity  is  many  orders 
of  magnitude  greater  than  water.  The  correlation  time  for  0.1  micron  spheres 
in  water  is  circa  37  X 10  ^ seconds.  The  apparatus  is  capable  of  determining 
correlation  times  some  six  or  seven  orders  of  magnitude  longer,  hence  for 
correspondingly  higher  viscosities.  It  may  also  be  possible  to  use  spheres  of 
about  1 order  of  magnitude  smaller  to  extend  the  method  to  even  higher 
viscosities. 
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Figure  1 Block  Diagram  of  the  Homodyne  Spectrometer. 
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ABSTRACT 


It  has  been  shown  experimentally  by  various  investigations  that 
the  dynamic  viscosity  of  a dilute  polymer  solution  in  the  limit  of  high 
frequency  does  not  become  the  viscosity  of  the  solvent.  To  study  this 
limiting  behavior,  a high  frequency  surface  delay  line  device  has  been 
developed.  The  pulsed  ultrasonic  shear  wave  propagated  into  the  liquid 
over  the  frequency  range  from  3 MHz  to  30  mz  and  has  been  used  to  deter- 
mine the  viscoelastic  properties  of  polymer  solutions.  A computerized  data 
acquisition  system  has  been  interfaced  with  this  device  to  determine  wave 
attenuation  and  phase  shift  even  for  pulsed  radio  frequency  signals.  The 
method  of  data  acquisition  is  based  on  adaptation  of  a conventional  high- 
speed sample  and  hold  circuit,  the  generation  of  precise  synthesizer 
locked  trigger  pulse  train,  and  the  utilization  of  "walking  frequencies"1’23. 
Data  analysis  is  based  on  cross-correlation  among  echoes  to  yield  phase 
and  amplitude.  Differences  in  phase  and  amplitude  with  and  without  solu- 
tions are  converted  to  the  viscoelastic  properties. 

The  results  of  the  high  frequency  limiting  behavior  of  several 
polystyrene  solutions  (5%,  10%,  20%,  40%)  together  with  steady-state 
viscosity  results  are  reported,  with  these  results  the  possible  explana- 
tion of  the  contribution  of  polymer  molecules  to  the  total  steady-state 

viscosity  is  presented  and  compared  to  the  theories  of  Rouse,4’25,  zimm5 

20 

and  others 
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PART  I 


INTRODUCTION 

The  viscoelastic  response  of  dilute  polymer  solution  is  usually 
separated  into  two  parts,  the  polymer  and  solvent  contributions.  The 
theories  of  Rouse,  Zimm  and  others  attribute  the  polymer  contribution  to 
interned,  relaxational  modes  of  the  macromolecules. 

At  low  frequencies,  the  long-range  correlated  segmental  motions 
dominate  the  viscoelastic  response.  The  low-frequency  dynamic  response 
has  been  investigated  by  a number  of  people  using  the  Bimboim  apparatus. 
Thereby,  it  has  shown  that  the  behavior  can  be  described  rather  well  by 

the  theories  of  Rouse,  Zimm  and  others  as  well  as  experimentally  by 

, 6 
Bimboim  and  Ferry  . 

The  high  frequency  response  in  dilute  polymer  solution  has  been 
explored  by  several  investigators  recently.  On  the  basis  of  the  theories 
of  Rouse  and  Zimm  one  might  expect  the  polymer  contribution,  evident  at 
lower  frequencies,  to  drop  to  zero  in  the  limit  of  high  frequencies,  and 
thus  the  viscosity  of  the  solution  to  approach  the  viscosity  of  the 
solvent  in  the  limit  of  high  frequency.  However,  it  has  been  found  ex- 
perimentally that  this  assumption  is  not  borne  out  and  the  difference 
between  the  limiting  solution  viscosity  and  that  of  the  solvent  has  been 
attributed  to  internal  friction  within  the  polymer  molecule.  Since  the 
higher  frequency  modes  of  motion  are  short  range  in  nature,  they  could 

be  greatly  influenced  by  the  internal  viscosity  of  the  submolecules. 

3 4 5 27 

This  interpretation  has  been  brought  into  question  by  Schrag  * ’ ’ and 
the  problem  remains  open.  This  discrepancy  between  theory  and  experiment 
in  the  high  frequency  limit  is  called  the  problem. 

81 


% 


The  experimental  techniques  that  have  been  used  to  investigate 
this  phenomena  have  mainly  been  in  the  frequency  range  from  0.01  Hz  to 
10  KHz.  This  means  that  it  has  been  necessary  to  use  highly  viscous 
solvents  in  order  to  shift  the  T^-phenomenon  into  experimentally  accessible 
regimes.  This  shifting  normally  requires  an  assumption  of  applicability 
of  the  principle  of  time- temperature  superposition. 

We  now  present  a brief  discussion  on  the  historical  development 

of  shear  property  measurements  in  liquids.  A great  deal  of  attention  has 

been  given  to  lower  frequency  experiments  by  using  a quartz  crystal  vi- 
24 

brating  in  torsion  . This  was  used  to  measure  the  complex  shear  vis- 
cosity of  the  liquid  surrounding  the  crystal.  The  change  in  resonant 
frequency  and  half  width  due  to  adding  the  liquid  are  measured  in  terms 
of  impedance  change  by  a bridge.  When  the  shear  impedance  of  the  liquid 
is  quite  high,  a different  pulse  method  employing  torsional  waves  travel- 
ing along  a delay  rod  can  be  used.  Hie  phase  and  amplitude  of  waves 
reflected  back  from  the  bottom  end  or  the  rod  are  then  measured.  A 
liquid  sample  is  then  placed  around  a rod  and  then  the  change  of  the 
phase  and  amplitude  is  determined.  In  another  technique,  the  first 
thickness  shear  modes  of  an  infinite  plate  were  used  in  a nondispersive 

, 9 

delay  line  by  Meitzler  . This  mode  was  propagated  in  a thin  rectangular 
strip  delay  line.  The  edges  were  taped  with  absorbing  tape  to  damp  edge 
generated  modes,  and  the  SSO  mode  was  excited  piezoelectrically.  Also, 
Meitzler  demonstrated  that  the  shear  mode  could  be  propagated  in  a finite 
plate  10  to  20  wavelengths  in  width  if  absorbing  material  is  placed  on 
a minor  surface.  These  lines  are  constructed  of  a metal  strip  with 
piezoelectric  transducers  bonded  to  the  ends.  Hie  crystals  can  excite 
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only  the  thickness-shear  mode,  which  is  zeroth  mode,  and  thus  a delay 

line  free  of  unwanted  inodes  can  be  achieved.  In  designing  the  zeroth 

mode  delay  line,  the  frequency-thickness  product  is  kept  below  that  for 

cut-off  of  the  first  dispersive  mode.  This  cut-off  occurs  when  the 

thickness  equals  one  wavelength  or  V /2f  and  as  the  resonant  frequency 

s 

of  these  transducers  has  no  width  dependence,  the  transducer  width  cam 
be  equal  to  the  strip  thickness. 

Hie  polygonal  delay  material  was  tried  for  high  frequency 
measurement  here  in  the  laboratory,  but  it  has  been  found  experimentally 
that  the  assumption  of  only  shear  wave  propagation  was  not  correct. 

There  existed  longitudinal  wave  contributions  which  were  generated  at  the 
reflecting  boundaries,  since  the  cut  angle  of  the  delay  material  at  the 
reflecting  boundary  did  not  satisfy  the  critical  angle  condition.  Wie 
critical  angle  is  the  angle  at  which  only  the  shear  wave  can  continuously 
propagate.  Experimentally,  it  is  difficult  to  fabricate  this  critical 
angle  in  delay  substrate  with  perfection.  Therefore  this  has  been 
abandoned. 

Hence,  a major  objective  in  this  research  has  been  to  develop 
a device  and  method  for  measuring  the  viscoelastic  response  of  low 
viscosity  liquids  in  the  frequency  range  2MHz  to  30  MHz.  The  device 
consists  of  a glass  plate  in  which  an  RF  ultrasonic  transverse  wave  is 
propagated  down  the  plate,  while  the  liquid  is  in  contact  with  the  major 
surfaces.  The  transverse  mode  of  vibration  which  is  polarized  in  the 
same  direction  as  major  surfaces  is  called  the  Meitzler  mode  as  mentioned 
above.  The  same  mode  was  used  by  Meeker  and  Meitzler  in  an  aluminum 
strip  delay  line  that  was  subsequently  applied  by  Hnuston11  and  Meyers 
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to  polymer  solutions.  Hie  technique  described  here  differs  also  from  the 
high  frequency  technique  of  McSkimin10.  The  wave  propagates  in  the  plane 
of  the  solid/liquid  interface  instead  of  at  a normal  or  oblique  incidence 
to  the  interface  as  in  the  two  methods  of  McSkimin. 

It  should  be  understood  that  all  these  measurements  deal  with 
the  mechanical  response  of  the  sample  under  very  small  deformations 
where  the  phenomenological  theory  of  linear  viscoelasticity  is  valid. 

In  the  apparatus  we  shall  be  measuring  the  difference  in  attenuation  (AA) 
and  difference  in  phase  (Ab)  in  the  absence  of  liquid  and  in  the  presence 
of  liquid.  Hien  the  mechanical  impedance  Z*  of  the  liquid  is  related  to 
the  above  AA  and  AB  by  the  following  formula:  Z*  = 2bpv  (AA  + iAB). 

Also,  the  quantity  G*,  called  the  complex  shear  modulus,  is  related  to  Z* 
by  equation  Z*  = ,/pG*.  Hie  quantities  Z*  and  G*  are  functions  of  frequency. 

For  the  data  acquisition,  a cross-correlation  and  walking  fre- 
quency method  has  been  utilized  to  obtain  the  phase/amplitude  measure- 
ments on  pulsed  RF  echoes  to  yield  AA  and  AB. 
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PART  II 
APPARATUS 

A.  Theory 

In  this  experiment  the  main  interest  is  measuring  the  dynamic 
shear  viscosity  T|#  of  the  liquid  sample  which  is  in  contact  with  major 
surfaces  of  the  delay  line.  The  equations  for  propagation  of  ultrasonic 
waves  in  a plate  are  developed  in  this  section  and  then  restricted  to  the 
case  of  interest,  namely,  the  lowest  order  transverse  mode.  This  mode  is 
often  called  the  Meitzler  mode^ . We  obtain  solutions  for  the  velocity 
and  attenuation  of  this  wave  in  the  presence  and  absence  of  liquid  at  the 
interface.  Equations  for  the  relationship  between  velocity  and  attenua- 
tion change  and  the  mechanical  impedance  and  the  shear  modulus  of  the 
liquid  are  also  obtained. 

Consider  an  infinite  plate  (Fig.l)  of  thickness  2b  in  a coordi- 
nate system,  where  it  is  bounded  by  the  major  planes  at  X = b and  X = - b 
with  a liquid  sample.  When  an  ultrasonic  transverse  wave  propagates  down 
the  plate  in  the  Z-direction,  a wave  is  also  generated  in  the  liquid  be- 
cause of  the  no-slip  condition  between  solid  plate  and  liquid  at  the 
interface. 

Oie  steady-state  sinusoidal  excitation  with  the  appropriate 
boundary  conditions  results  in  a set  of  homogeneous  equations  whose  solu- 
tions have  discrete  modes  of  wave  propagation. 

For  the  delay  line  used  here,  the  dimensions  and  polarization 
of  the  piezoelectric  (P.Z.T.  ceramic  transducer  material)  are  chosen  so 
that  only  the  principal  shear  mode  (Meitzler  mode)  of  propagation  is  ex- 
cited while  higher  shear  modes  are  suppressed  (when  the  plate  is  in 
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vacuum).  It  is  very  useful  to  have  only  the  nondispersive  Meitzler  mode 
SS  (O)  because  its  velocity  is  independent  of  frequency  and  the  complica- 
tion of  having  several  modes  of  different  velocities  can  be  avoided. 


Even  in  the  presence  of  liquid  on  the  major  surfaces  it  will  be 
assumed  that  only  the  same  SS  (O)  mode  is  present  since  the  mechanical 
impedance  of  the  liquid  is  much  less  than  the  solid. 

1 . Equations  for  Shear  Waves 

The  equation  of  wave  propagation  for  an  isotropic  homogeneous 
linear  viscoelastic  material  is  generally  written  as  follows: 

A2u  + ( \*  + p*)A(A  • u)  = -pcA  (1) 

where 

(j.*  - complex  shear  modulus 
X*  - other  Lame  coefficient 

u - total  particle  displacement  U,  = u elwt 

H (x,y,  z,t)  (x,y,  z)C 

u - particle  displacement  vector 

p - density. 

If  one  uses  Hamilton's  decomposition  method  by  defining 

u = A<£  + (A  x W)  (2) 

in  terms  of  scalar  and  vector  potentials  <p  and  W,  then  the  following  wave 
equations  in  <f>  and  W are  obtained  from  (1): 


..2  2 , 2, 

V^A  0 = - (o  0 

(3) 

VgA?'W  = - a>  W 

(4) 

with 


A*  W 


where 


2U*  + >* 


in  which  Vg  and  VD  are  the  transverse  and  dilatational  wave  velocities. 

In  our  delay  line  device,  the  transducers  are  polarized  to  provide 
displacements  with  a component  in  the  y- direction  only,  hence  we  assume 
Ux  = Uz  = °*  Thus»  ^ wave  has  displacement  Uy  = Uy(x,  0,  z,  t)  in  the  y- 
direction  only;  the  displacement  is  assumed  to  be  the  same  everywhere 
in  y.  Thus,  Eqs.  (3)  and  (4)  reduce  to: 

*A2  2 

^*AUy  = - pa,  uy  (6) 


^ 2 

\ dx 


2 

- P<i)  u 


Equation  (7)  is  a shear  wave  equation  with  complex  shear  modulus 
p*  (one  of  Lame's  constants).  By  the  application  of  the  method  of  separa- 
tion of  variables  to  Eq.  (7),  namely  Uy(x,0,z,t)  = U (x.O.zJe1  t,  one 
obtains  for  the  total  solution  of  Eq.  (7) 

Uy.  = A(CemZ  + De'  mZ)(EenX  + Fe"  °X)  (8) 

where  cj  is  the  angular  frequency,  A,  C,  D,  E,  F are  constants  to  be  de- 
termined by  the  boundary  conditions  and  m and  n are  the  complex  propa- 
gation constants  in  the  z-  and  x-directions,  respectively. 
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Displacement  Boundary  Conditions 


* 


This  same  equation  can  now  be  applied  to  each  of  the  three  regie 
of  interest:  the  solid  delay  line  and  the  upper  and  lower  liquid  layers; 
then  the  boundary  conditions  are  matched.  In  steady  state  the  frequencies 
are  the  same  in  all  regions. 

Consider  the  transverse  wave  in  the  viscoelastic  medium  above  and 
below  the  plates.  This  wave  is  traveling  unbounded  in  the  positive  z-  and 
x- directions  and  since  it  is  traveling  in  an  attenuating  medium,  the  dis- 
placement must  approach  zero  as  z and  x approach  infinity.  Then 
C = E = 0 in  the  liquid  above  the  plate  and  similarly  C = F = 0 in  the  liquid 
below  the  plate. 


In  the  plate  itself  the  waves  are  bounded  in  the  x- direction  and  thus 
there  is  only  one  boundary  condition  that  is  C = 0.  The  solutions  can  be 
obtained  with  symmetric  boundary  condition  in  x-direction  as  follows; 


u 


V 

V 


m_z  - n„x 


„ - m,z  . n,x  - n,x4 

He  1 (e  1 + e 1 ) 


H e-  m°Z  en°x 


(9a) 

(9b) 

(9c) 


where  the  subscripts  0,  1, 2 designate  the  region  above  the  plate,  the 
region  within  the  plate  and  the  region  below  the  plate,  respectively.  The 
symbols  Hq  and  are  appropriate  combinations  of  A,  C,  D,  E,  F. 


3.  Stress  Boundary  Conditions 

In  the  presence  of  liquid  the  stress  continuity  across  the  plane 
boundary,  x = ±b,  requires  that 

a =0  (10a) 

xx 
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dU 


xyU,  y»  t)  = p* 


a = 0. 

xz 


a ~ iwt 

Xy(s,y,  z.t)  CT\y(x,y,  z)e  (10b) 


(10c) 


If  the  plate  is  in  vacuum,  all  component  stresses  across  the  plate 

duv 


surface  must  be  zero,  hence  from  a = p*  = 0 at  x = ±b  one  obtains 


xy 


Thus  real  [n^]  = 0 and 


IMfn^ 


iijb 


IT  I 


= e 


-nlb 


, I = 0,1,2, 


(11) 


(12) 


If  we  now  restrict  our  consideration  to  the  solution  with  1 = 0 (the 
Meitzler  mode),  Eq.  (10)  for  the  wave  in  the  plate  becomes 


U = 2H  e -miVZ 
yl  1 


(13) 


where  and  m^  will  refer  to  the  propagation  vector  in  the  solid,  in 

the  z-direction,  in  the  absence  and  presence  of  liquid. 

If  the  plate  is  in  a viscoelastic  liquid,  the  complex  shear  modulus 
of  the  liquid  p*  is  obtained  from  the  stress  continuity  equation  (10b)  com- 
bined with  Eq.  (9b)  to  yield 


nj(  1 “ e 


2njb) 


1 (i  i 2n,b 
n ( 1 + e 1 ) 
o 


(14) 


since  ( a ) = (a  ),,  at  x = b. 

xy  o xy'  1 

Next,  the  following  wave  vector  for  the  liquid  is  obtained  from 
Eqs.  (9z)  and  (7) 

2 2 

m + n + — — = 0 
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(15) 


From  displacement  continuity  across  the  plate,  u = u at  x = +b 

Vo  Vl 

for  all  values  of  z.  Thus,  by  equating  Eqs.  (9a)  and  (9b),  one  conclues 


that 


m = m,  _ 
o 1L 


(16) 


Hence  Eq.  (15)  can  be  combined  with  (16)  to  obtain  an  equation  that  couples 
wave  vectors  in  the  liquid  and  in  the  solid. 


2 2 

m,  + n + 
1L.  o 


at  p 


= 0 


(17) 


Next,  the  following  wave  vector  relation  for  the  plate  (in  liquid  con- 
tact) is  obtained  from  Eqs.  (9b)  and  (7) 


at  p 


m + n + 
1L  1 


i 


= o 


(18) 


While  for  the  plate  m vacuum,  we  have  shown  by  Eq.  (12)  that  n^  = 0 for 
the  Meitzler  mode,  thus,  the  wave  vector  relation  (18)  becomes 


2 

co  p. 


miv  + 


= 0 


(19) 


In  our  experiment  is  measured  in  air  (i.e.,  vacuum)  and  again 


in  liquid.  Thus,  by  combining  Eqs.  (18)  and  (19)  we  obtain 

2 


2 2 
miL-  miV 


= - n 


1 


(20) 


The  complex  modulus  of  the  viscoelastic  liquid  |a*  can  be  obtained 
by  using  the  experimental  values  of  Eq.  (20)  together  with  Eqs.  (14)  and 
(17)  to  give  a quadratic  equation  for  |j,*.  Alternatively,  an  ezplicit  rela- 
tion to  the  mechanical  impedance  of  the  liquid  z*,  can  be  obtained  as  seen 
below. 


Notice  some  assumptions  that  were  made  in  the  preceding  discussion: 
(1)  Motions  are  steady- state  sinusoidal  of  small  amplitude  so  that 
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no  nonlinear  response  of  the  propagating  media  can  be  neglected;  (2)  the 
plate  is  assumed  to  be  isotropic  and  homogeneous  and  (3)  the  transducer 
will  only  excite  the  nondispersive  SS(0)  mode. 


4.  Relation  of  Mechanical  Impedance  of  Liquid  to  Propagation  Vector 
in  Solid 


The  equations  derived  in  the  preceding  section  can  also  be  related 

26 

to  the  mechanical  impedance  Z*  of  the  liquid 

The  mechanical  impedance  of  a material  is  defined  as  the  ratio  of 
the  stress  to  the  particle  velocity  produced  by  that  stress  at  any  surface 
through  the  material.  For  the  shear  stress  components  of  interest  here 
in  the  liquid  we  can  write 


or 


Z* 

o 


gxyo(x,  z) 
io»U  , . 

yo(x,  z) 


a - icoZ*  U 
xyo  o yo 


(21a) 


(21b) 


It  can  be  shown  that  the  characteristic  (shear)  mechanical  impedance 
for  an  unbounded  liquid  is  related  to  the  complex  shear  modulus  and 
density  by 


Z* 

o 


'V?1'2 


(22) 


Across  the  liquid-plate  interface,  x = +b,  we  have  seen  that  dis- 
placement and  stress  must  be  continuous,  hence  velocity  and  stress/ 
velocity  must  also  be  continuous. 


At  x = b 


U 

yo 


(2  3a) 


U 

yo 


(23b) 
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or 


xyl 


gxy  1 

Uyi 


o = Z*  U 
xyo  o yo 


7=2  = i(uZ* 

U o 


(23c) 


(23d) 


yo 


Since  we  had  for  displacement  and  stress  in  the  plate 

Jyi 


, - m,  Z icjt 

U . = 2H^cosh  n^x  e 1 e 


(24) 


and 


, - m,Z  io>t 

a , = - 2H  lj.  n sinh  n x e 1 e 

xyl  11  1 1 


(25) 


the  values  of  these  functions  at  x = b are  substituted  into  Eq.  (23d)  to 
obtain 


Z*  = i tanhnb 

o^l 


(26) 


If  we  expand  tanh  n^b,  and  take  | n^bj  « 1, 


(n1b) 

tanh  n^b  = n^b  + + 


then 


Z*  = 
o 


^"ib 


Ci) 


We  had  Eq.  (20),  that  is, 


2 2 2 
-ni  = miL'mlV 


(27) 


(28) 


(20) 


It  is  convenient  to  write  the  complex  propagation  constants  in  their 
component  form: 

(29a) 


and 


m,  __  = A + iB 

IV  o o 


m = A + iB 
1 L 


(29b) 


where  A , B and  A,  B are  the  attenuation  and  propagation  constants  in 
o o 

the  plate  in  the  absence  and  presence  of  liquid,  respectively.  Hence, 
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from  Eq.  (20) 


= (A  + iB)2  - (A  + iB  )2 
•i  o o 


For  the  usual  case*  where 


Eq.  (30)  becomes 


(B  + B ) » (A  + A ) 
o o 


= 2Bq(-  AB  + iAA) 


where 


A A = A - A and  AB  = B - B. 
o 


and  from  Eqs.  (19)  and  (29a) 


2 1/2 


Now  Eq.  (28)  is  rewritten  as 

Z*  = 2 V b(AA  + iA  B)  (35) 

o s l ' 

This  is  the  basic  experimental  equation  that  relates  the  mechanical 

impedance  of  the  liquid  to  the  measured  differences  in  attenuation  in  the 

plate  delay  line. 


Assume  B^  = B,  and  from  Eqs.  (31), 

(A  + iB)2  - (A  + iB  )2  = A2-  A2  - B2  + B2  + 2i( AB  - A B ) 
o o o o o o 


= (A  - A )(A  + A ) - (B  - B )(B  + B ) + 2i( AB  - A B ) 
o o o o o o 


=(  B + Bo) 


A + A AA  - B 

----  ° - AB  + 2i  — 


(*•  T ■*» 

^bTI2- 

n 


- B \ 

nr 


=*2Bq  (-  AB  + iAA) 
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5.  Relationships  of  Viscoelastic  Constants z 

It  is  convenient  to  summarize  the  interrelation  between  the 
mechanical  impedance  and  the  complex  shear  modulus  parameters  for  descrip- 
tion of  the  liquid.  These  relationships  are  based  on  the  linear  theory 
of  viscoelasticity.  Also  note  that  in  viscoelastic  convention,  the 
symbol  G*  = |i*  for  the  shear  modulus: 


2*  5 ft  + l * (36) 

<**£/<*' = <k'  + LS"  (37) 

l«rl  s c 1381 

z*  = (/*Cr*)X  , (39’ 

R = f.  & ( 1*1  -C^)/x  ,40> 

X = ?, Yx  ( \gl  <411 

G,'  = C R *-X*)/£  <43> 

q"  =.  zrx  /ft  (43) 

Gr*  =T  (44) 

y,*  « , „ (45) 

where  R and  X are  mechanical  resistance,  pQ  the  liquid  density,  G7  and  G* 
are  the  shear  storage  and  loss  moduli,  T]/  and  1\"  are  the  real  and  imag- 
inary components  of  dynamic  viscosity. 


B.  Description  of  Delay  Line 


This  instrument  is  intended  for  the  use  of  low  viscosity  samples, 
particularly  dilute  polymer  solutions,  in  the  frequency  region  from 
about  3 mz  to  30  MHz. 


1.  Transducer.  Delay  Material  and  Equations 

The  piezoelectric  transducer  will  vibrate  when  subjected  to  an 
alternating  electrical  field  whose  frequency  is  within  the  bandwidth  of 
the  transducer.  The  center  frequency  of  the  thickness  shear  transducer 
depends  on  its  thickness  and  choice  of  piezoelectric  material.  The  P. z. T. 
ceramic  transducers,  whose  center  frequency  is  20  MHz,  are  rectangular 
bars  with  the  width  2b  (b  = 0.030"),  the  same  as  the  thickness  of  the 
delay  line,  of  length  one  third  the  width  of  delay  line,  and  of  thickness 
0.007"  as  shown  in  Fig.l.  This  transducer  is  designed  to  drive  only  the 
nondispersive  SS (O)  mode  of  propagation  in  the  delay  line.  The  polari- 
zation of  the  transducer  motion  is  in  the  y-direction,  i.e.,  parallel  to 
the  major  surface  of  the  delay  line.  Therefore,  the  displacement  exists 
only  in  the  y-direction.  The  propagation  of  the  wave  is  in  the  z- 
direction. 


The  delay  medium  of  low  loss  fused  silica  (quartz)  is  three 

inches  long,  two  inches  wide  and  .060"  thick.  The  smoothness  of  the  major 

surface  of  delay  medium  is  extremely  important.  For  example,  consider 

-2 

a Newtonian  liquid  of  viscosity  10  poise.  The  damping  length  in  centi- 

a 


meters  for  a viscoelastic  liquid  is  given  by  (|G*|  -G)J  ; for 


2 1 G 

Newtonian  liquid  the  damping  length  becomes  (U)^p/2G* ) . 

above  liquid  at  30  MHz  the  damping  length  is  1 X 10~5  cm.  At  the  same 


Thus  for  the 
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frequency  the  wavelength  in  the  solid  is  1 X 10  cm.  The  surface  irregu- 


larities should  be  kept  much  smaller  than  10  cm.  Surface  smoothness  of 


the  delay  medium  reported  herein  have  not  been  established  in  spite  of  the 


claims  made  by  manufacturers  that  they  have  done  high  quality  mechanical 


polishing  on  the  major  surfaces  of  the  delay  medium.  The  small  streaks  of 


the  random  lines  are  observed  under  20  X magnification  microscope.  Hie 


streak  of  lines  on  the  major  surface  might  have  been  caused  by  a careless 


polishing  process.  The  nominal  delay  time  of  the  delay  line  is  20  p,  sec. 


The  transducers  are  mounted  at  both  ends  of  the  delay  material. 


A transducer  at  one  end  is  the  driving  transducer  and  a transducer  at  the 


other  end  is  the  receiving  transducer.  Hie  mounting  of  the  transducer  to 


the  delay  material  is  critically  important  since  the  transducer  thickness 


is  only  0.007"  thick.  The  thickness  of  the  bond  between  transducer  and 


cross-section  of  the  delay  material  has  to  be  thin  compared  to  a wave- 


length (0.014")  so  that  wave  reflection  will  not  occur  between  the  inter- 


faces; the  adhesive  has  to  be  spread  uniformly  over  the  whole  area  in 


order  to  produce  uniform  phase  shift  of  the  mechanical  signal.  A thin 


lead  backing  is  also  important  for  reflection  consideration  and  is  neces- 


sary in  order  to  make  the  electrical  connection  to  the  backside  of  the 


transducer.  Hie  connection  to  the  transducer  is  by  a very  thin  flexible 


wire  so  that  it  does  not  exert  stress  on  the  transducer  and  also  minimize 


the  acoustic  wave  reflection.  Hie  details  are  shown  in  Fig. 2. 


The  coupled  wave  solutions  for  displacement  in  the  y-direction  u, 


in  the  solid  and  in  the  liquid  are  respectively  given  by 


usc*,r,*)=  H* exp  [-«£-<■  ■+<<•*]  c.s  (£ 


X(46) 


I 


■>  - ^ 


Ul  (*,*,*)=  HLexp[-i(^'->*'^0^  +c«t] «*p[-  X J 

Through  Eq.  (46),  the  difference  in  attenuation  AA  and  in  phase  AB  of  the 
wave  measured  at  Z **  l in  the  absence  and  in  the  presence  of  liquid  can 
be  related  to  the  mechanical  z*  of  the  liquid  by 


(47) 


AA  + = 


l 


2bf>V 


(48) 


where  p and  p are  the  density  of  the  solid  and  the  liquid  and  V is  the 

Xj 

velocity  in  the  solid  with  liquid.  Hie  mechanical  impedance  Z*=R  + iX 
is  related  to  the  dynamic  shear  moduli  G*  = G*  + xg"  = iouT]7  by 

c/ s i?-*1,  ; Gf"  - z R x. 

For  a Newtonian  liquid  of  viscosity  7),  G*  =*  0,  thus 

r = x = ^ y* 


(49) 


(50) 


2.  Delay  Line  Holder 

The  delay  line  is  mounted  in  a polycarbonate  cell  which  can  be 
easily  inserted  in  a constant  temperature  bath  as  shown  in  Fig. 3.  Hie 
electrical  connections  are  made  of  coaxial  cables  which  have  waterproof 
teflon  insulation,  and  the  cables  are  terminated  in  o-ring  clamps  (Fig. 3) 
to  seal  them  into  the  polycarbonate.  Two  one-inch  thick  polycarbonate 
plates  are  used  to  make  the  delay  line  cell.  There  are  two  recessed  o- 
ring  grooves  in  one  plate,  and  one  groove  in  the  other.  Hie  large  outer 
o-ring  is  for  isolation  of  the  electrical  part  of  the  delay  line  from  the 
bath  water.  The  delay  line  is  suspended  between  the  two  inner  o-rings. 
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These  o-rings  also  serve  to  isolate  the  samples  to  be  tested.  That  is, 
the  delay  line  major  surfaces  and  holder  inside  the  wall  form  a compartment 
for  the  sample  to  be  tested.  To  form  a liquid  tight  seal  with  the  brittle 
fused  silica  delay  line  at  the  lowest  possible  stresses  (o-ring  compres- 
sion), a rather  soft  silicone  o-ring  is  used.  The  thickness  of  the  o-ring 
is  selected  carefully  so  that  it  is  large  enough  for  wave  penetration  and 
does  not  have  any  wave  reflection  in  the  liquid  sample.  Silicone  o-rings 
have  better  durometer  (softness)  than  those  of  other  materials.  Chemical 
silicone  was  satisfactory  for  water,  alcohols,  DOP  and  polys tyrene-DOP 
solution  but  swells  in  the  presence  of  toluene,  cyclohexane  and  benzene. 

Hie  o-rings  were  bought  from  the  Irving  B.  Moore  Company,  Cambridge, 
Massachusetts.  The  sizes  and  durometer  were  as  follows: 


o-ring 

1 

I.D.  = 2«j, 

w = 1/16 

#10 

Durometer  =60 

o-ring 

2 

rf 

II 

Q 

H 

W = 3/16 

# 3 

Durometer  = 70. 

There  are  two  holes  in  each  side  of  the  cell  plate  for  filling 
and  emptying  the  sample  liquid.  Syringes  are  connected  to  the  end  of 
teflon  tubings  which  extend  from  the  sample  compartment  inside  the  poly- 
carbonate cell  to  outside  of  the  temperature  bath.  Therefore,  the  volume 
of  the  sample  injected  into  both  major  surfaces  of  the  delay  line  can  be 
read  off  from  the  syringe  scale.  The  thickness  of  the  sample  cell  is 
0.015"  and  the  diameter  is  2.5".  Hie  volume  of  the  sample  liquid  in  each 
cell  becomes  0.74  m l.  However,  the  reading  of  the  syringe  scale  was 
approximately  2 ml.  Hiis  was  because  the  sample  liquid  was  also  filled  in 
the  length  of  the  teflon  tubing.  The  holes  for  sample  liquid  exchange  have 
to  be  lined  up  vertically  so  that  emptying  and  filling  of  the  sample  can 
be  carried  out  with  the  help  of  gravity.  Gravity  helps  avoid  the  air 
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bubble  forming  and  leaving  sample  liquid  trapped.  Also,  a small  weight 
of  1 lb  was  used  at  the  end  of  the  syringe  to  create  the  gentle  steady 
gradual  motion  for  emptying  and  filling  of  the  cell  with  liquid  sample. 

3.  Temperature  Control 

It  is  necessary  to  control  the  temperature  of  the  delay  line 
not  only  because  the  s simple  (liquid)  properties  are  temperature  dependent 
but  also  because  the  delay  medium  (fused  silica)  is  strongly  temperature 
sensitive. 

A Neslab  RTE-4  unit  was  used  to  control  the  temperature  of  the 
bath.  The  Z-pump  of  the  control  unit  enables  water  to  circulate  between 
the  bath  and  the  control  unit.  It  takes  approximately  two  hours  to  equi- 
librate the  sample  cell  in  the  bath  water  when  temperature  setting  of  the 
control  unit  is  lowered  or  raised  by  5°C.  The  total  circulating  distance 
between  the  bath  and  the  control  unit  is  about  6 ft.  The  volume  of  the 
bath  is  approximately  16,000  cc.  The  bath  is  constructed  of  thermopane 
double  glass  walls  mounted  in  an  aluminum  frame.  Several  coats  of  RTV 
adhesive  cure  used  to  seal  the  glass  sections  into  the  frame.  Hie  double 
glass  walls  provide  an  excellent  thermal  insulation  from  the  fluctuation 
of  ambient  temperature.  The  water  was  used  for  the  bath  liquid  and  the 
range  of  the  temperature  in  the  experiment  was  between  10° c and  25° C.  Hie 
temperature  fluctuation  of  the  bath  was  about  0.005°C  in  an  hour  and 
0.05°C  in  a day. 

4.  Thermometry 

A quartz  crystal  thermometer  was  used  to  measure  the  temperature 
of  the  bath  water.  A crystal  quartz  had  been  cut  in  such  a way  that  the 
resonant  frequency  changes  with  temperature.  When  the  quartz  crystal  is 
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adjoined  with  H.P.  2830A  oscillator  sensor,  the  G.R.  1192B  counter  is 
provided  to  read  the  output  of  this  thermometer.  The  quartz  crystal  has 
a rod-like  shape  which  is  immersed  in  the  bath. 


The  calibration  constant  for  this  thermometer  was  found  to  be 
0.87  KHz/C  when  it  was  calibrated  against  the  mercury  thermometer  as 
shown  in  Fig. 8.  This  value  is  not  in  agreement  with  the  coefficient  of 
0.9861  KHz/C  which  is  the  value  that  the  manufacturers  supplied.  No 
explanation  can  be  given  for  this  at  the  present  time. 


5.  Manufacturer  of  the  Delay  Line 

The  delay  line  was  fabricated  by  Andersen  Laboratories,  Inc., 
Bloomfield,  Connecticut.  The  transducer  and  the  delay  medium  are  P.Z.T. 
and  a fused  quartz,  respectively.  The  method  of  fabrication  was  the 
multilayer  transducer  mounting  with  an  optimum  energy  transduction  tech- 
nique as  shown  in  Fig. 2.  The  tinfoil  was  used  between  conductive  epoxy 
in  order  to  achieve  the  maximum  energy  couple  between  the  transducer  and 
the  delay  medium.  The  following  data  were  supplied  by  Andersen: 

Delay  Medium  - Fused  Quartz,  polished 

Particle  Motion  - Shear,  Parallel  to  Major  Surfaces 

Delay  - 20  p,s  nom. 

Loss,  Terminal-to-Terminal:  1 - 2 db 

Loss,  in  50  Q comparison  system:  - 14  db  (approx. ) 

Reflections:  -10  db  (approx.)  referred  to  previous  signals 

Transducers : 20  f«z  ceramic 

Cp  ~ 1500  - 2000  pfd 
Rp  (dynamic)  ~ 10  Q 

Density  of  Delay  Medium  - 2.2  g/cc 

Acoustic  Shear  Velocity  of  Delay  Medium  - 3.8x10^  cm/sec. 
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. Drive  and  Response  Electronics 


1 


The  method  used  here  is  a pulse-echo  technique,  This  method 
uses  a series  of  R. F.  drive  pulses  generated  at  a relatively  slow  repeti- 
tion rate  so  that  when  a drive  pulse  is  generated,  all  of  the  vibration 
resulting  from  the  previous  drive  pulse  has  been  completely  attenuated. 
Each  of  the  pulses  is  composed  of  a train  of  sinusoidal  vibration.  Its 
nominal  width  is  20  p,sec  long  and  the  repetition  interval  is  1280  p,sec. 

The  pulse  method  has  a very  important  advantage  of  a low  duty 
cycle;  since  the  pulse  rate  is  relatively  slow,  only  a small  amount  of 
energy  is  added  to  the  system.  This  phenomena  is  important  in  order  to 
keep  the  temperature  from  fluctuating  in  the  device;  the  velocity  of  the 
wave  is  strongly  dependent  on  the  temperature. 

These  input  pulses  to  delay  line  are  generated  electrically  by 
an  R. F.  pulser  which  was  designed  to  give  the  proper  repetition  rate  and 
pulse  width.  The  C.W.  input  to  this  pulser  originates  from  a Hewlett 
Packard  5110A/5100A  frequency  synthesizer.  A precise  40  M-sec  clocking 
pulse  for  the  pulser  is  generated  by  the  programmable  clock;  this  interval 
is  based  on  the  same  stable  1 MHz  oscilliscope  which  drives  the  frequency 
synthesizer. 

The  R. F.  pulser  is  designed  in  such  a way  that  the  pulse  width, 
repetition  rate  and  pulse  delay  time  cam  be  adjusted  as  follows : the 
pulse  width  is  set  by  a monostable,  adjusting  time  constant  with  an 
external  potentiometer  R^.  The  repetition  rate  is  set  by  digital  multi- 
plier of  the  40  p,sec  pulse  interval.  The  pulse  delay  is  set  precisely 
to  am  integral  number  of  1 p,sec  pulse.  The  detailed  circuit  diagram  and 
pulse  sequence  are  shown  in  Figs. 5 and  6. 
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These  R.F.  drive  pulses  are  transformed  by  the  input  piezoelectric 
transducer  into  acoustic  waves  which  are  echoed  in  the  delay  line  until 
they  have  been  completely  attenuated.  These  echoes  in  the  delay  line  will 
exhibit  an  exponential  decay  as  shown  in  Fig. 2.  When  driven  at  the  center 
frequency  of  circa  20  MHz,  more  than  ten  echoes  have  been  observed.  The 
sample  and  hold  circuit  is  a Tektronic  586  sampling  unit,  where  the  sampling 
part  of  the  circuit  and  the  external  trigger  mode  were  utilized  in  order 
to  sample  the  incoming  pulsed  radio  frequency  signal  from  the  device.  This 
sample- and-hold  circuit  is  triggered  every  40  n,sec  by  a fast  pulse  which 
is  produced  by  the  combination  of  the  mark  pulse  from  the  mono  (^)  of  the 
R.F.  pulse  and  the  same  1 MHz  reference  oscillator  that  drives  the  fre- 
quency synthesizer.  These  sampled  analog  signals  are  then  fed  into  an 
analog-to-digital  converter  and  stored  in  the  memory  of  GE/PAC-4020  data 
processing  computer. 

In  order  to  ensure  the  precise  time  synchronization,  all  the 
electronics  involved  in  data  acquisition  process  have  to  be  referenced  to 
a single  clock,  i. e. , the  1 MHz  reference  oscillator  in  the  frequency 
synthesizer.  The  phase  change  produced  by  adding  the  sample  is  so  very 
small  that  jitter  in  the  sampling  circuit  trigger  can  lead  to  erroneous 
results.  Therefore  this  experiment  involves  the  use  of  a very  stable 
external  electronic  source.  In  order  to  sample  a high  frequency  train  of 
sine  waves  and  then  reconstruct  the  sampled  data  into  a single  sine  wave, 
the  sampling  circuit  must  satisfy  two  basic  conditions.  First,  the  aperture 
time  window  of  the  sampling  circuit  must  be  very  short  compared  to  the 
period  of  the  high  frequency  signal  in  order  to  measure  an  "instantaneous 
voltage".  Secondly,  the  time  interval  between  sampling  pulses  must  be 
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precisely  known,  and  must  have  little  jitter  in  order  to  maintain  phase 
coherence  between  the  train  of  sampling  pulses  and  the  high  frequency 
sine  wave  of  interest.  In  order  to  achieve  these  requirements,  a snap- 
diode  is  used  to  convert  the  standard  1 MHz  sine  wave  to  a 1 MHz  train  of 
narrow  pulses  and  a gate  selects  each  40th  pulse  as  the  sampling  pulse. 
These  pulses  are  a few  nanoseconds  wide  fast  pulses. 

An  HP  46 1A  broad-band  amplifier  is  used  at  the  output  of  the 
delay  line,  and  either  20  db  or  40  db  of  gain  were  used  according  to  level 
of  the  output  signal  of  the  delay  line.  Hie  output  of  the  amplifier  is 
then  connected  to  B input  of  a model  4S1  plug-in  unit  of  sampling  scope 
with  the  vertical  gain  set  at  10  mv/cm.  The  A input  of  the  plug-in  unit 
is  derived  from  the  output  drive  frequency  of  the  synthesizer  with  vertical 
gain  set  at  100  mv/cm.  This  A input  signal  is  C.W.,  and  is  tapped  at  the 
voltage  divider  just  before  the  mixer. 

Hie  two  "hold"  outputs  A and  B of  the  sampling  scope  are  then 
fed  into  the  analog- to-digital  converter.  The  differential  amplifier  input 
terminals  of  the  analog-to-digital  converter  are  connected  to  the  single 
ended  output  of  the  sampling  scope.  Hie  outer  shield  is  grounded  at  the 
scope  only,  and  a load  resistor  is  connected  between  the  negative  side 
and  ground  of  the  input  terminals  of  the  ADC.  The  detailed  electronics 
sure  shown  in  Fig.  4. 

1.  R.F,  Pulser 

The  program  clock  is  set  to  produce  a train  of  1 nsec  wide  pulses 
every  40  nsec,  upon  computer  program  initialization  signal  to  the  pulser 
16  clock  pulses  cure  counted  in  the  pulser  to  produce  one  gate  pulse  after 
640  nsec,  and  subsequently  every  32  clock  pulses  are  counted  to  produce 
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gate  pulses  at  1280  p.sec  intervals  thereafter.  Monostable  #1  is  adjusted 
to  set  the  gate  pulse  width  to  nominal  20  p,sec. 


The  gate  pulse  is  mixed  with  the  R. F.  signal  to  produce  the 
R. F.  pulse  for  input  to  the  delay  line.  The  same  output  pulse  from  the 
counter  is  used  to  trigger  the  adjustable  monostable  #2  in  order  to  pro- 
duce a variable  delay  time.  In  turn  this  output  triggers  monostable  #3 
to  produce  a 0.9  p,sec  wide  mask  pulse  (P3).  The  width  of  this  mask  pulse 
will  permit  only  one  of  the  1 MHz  train  of  snap  diode  pulses  to  pass  and 
to  trigger  the  sample  and  hold  circuit  every  1280  p,sec.  The  pulse  width 
of  monostable  #3  is  set  by  the  R.C.  combination  of  1.271  kQ  and  10  3 n farad. 

The  purpose  of  the  (variable)  delayed  mask  pulse,  (adjustment  of 
the  50  kQ  potentiometer  of  the  R.C.  combination  in  the  external  timing 
circuit  of  monostable  #2)  is  to  enable  us  to  select  only  one  delay  line 
echo  at  a time;  thus,  the  delay  time  of  the  monostable  #2  is  set  accurately 
so  that  the  mask  pulse  occurs  at  the  center  of  the  selected  output  echo  of 
the  delay  line. 

The  reset  pulse  produced  by  program  initialization  is  used  to 
reset  the  counter.  The  details  are  shown  in  Figs. 5 and  6. 

2.  Detection 

2.1  Amplifier.  A broad-band  amplifier,  H.P.  461  20  dB  - 40  dB  was  used 
to  amplify  the  low  level  response  signal  at  the  output  of  the  delay  line. 

Hie  gain  control  was  set  to  20  dB  and  the  output  signal  of  the  amplifier 
was  greater  than  40  mv.  We  assumed  that  this  amplifier  is  linear  and  that 
the  phase  shift  introduced  by  the  amplifier  is  the  same  for  all  the  pulses 
regardless  of  their  voltage  level  within  our  experimental  range. 
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2.2  Sampling  Scope.  As  indicated  in  the  previous  section,  a sampling 

scope  (Tektronix  model  661  with  model  461  vertical  amplifier  and  model  5T3 

timing  unit)  was  used  to  sample  the  signals,  A channel  and  B channel  of 

the  vertical  amplifiers  were  set  as  follows: 

Gain  control  A;  100  mv/cm 

Gain  control  B:  10  mv/cm 

5t3  timing  unit  was  set  as  follows: 

trigger  mode:  external,  DC,  50  Q 
slope:  negative 
sample  points/cm:  100. 

Under  these  conditions,  the  staircase  generator  of  the  timing 
unit  is  inactive,  and  the  input  pulse  directly  triggers  the  sampling  bridge 
circuit.  Hie  vertical  gain  controls  are  used  to  attenuate  the  driving 
and  response  signals  at  the  input  A and  the  input  B,  respectively.  The 
smoothness  control  is  adjusted  to  minimize  the  sampling  hysteresis.  In 
practice,  it  has  been  found  necessary  to  select  one  echo  at  a time  for 
sampling.  When  the  sampling  circuit  was  constantly  switched  between  one 
echo  and  another  excessive  hysteresis  noise  was  introduced.  This  phenomena 
is  probably  due  to  the  long  discharging  time  of  capacitor  at  the  sampling 
circuit. 

D.  Data  Acquisition  (Hardware  and  Software) 

1.  ADC  (Analog  to  Digital  Converter)  to  CPU 

Hie  output  voltages  from  the  sample  and  hold  circuit  of  the 
sampling  unit  cure  available  for  ADC  conversion.  Hie  ADC  digitizes  the 
analog  signal  and  the  digitally  coded  information  will  be  transferred  into 
the  memory  of  GEPAC-4020  computer.  Hie  same  pulse  recurring  every  40  M-sec, 
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that  is  used  to  trigger  the  sampling  unit,  is  also  used  to  initiate  ADC 
conversion,  to  execute  a data  transfer  to  memory,  and  causes  a computer 
interrupt.  Each  interrupt  initiates  a program  subroutine,  "READ  DATA" 
whose  function  it  is  to  organize  the  incoming  data  into  memory  as  described 
in  the  next  section. 


2.  Selection  and  Development  of  Walking  Frequencies  from 
Synchronous  Frequency 

The  method  of  walking  frequencies  as  shown  in  Fig. 7 has  been 
developed  for  the  purpose  of  transforming  a high  frequency  periodic  signal 
into  a low  frequency  signal  that  can  be  digitized  and  read  into  the  computer 
at  a rate  of  one  point  every  40  |isec. 

In  order  for  a complete  and  separate  sinusoidal  wave  to  be  re- 
constructed in  the  memory  from  each  echo,  the  synthesizer  must  always  be 
set  to  a walking  frequency  f given  by 


N 

where  f = — corresponds  to  any  synchronous  frequency  (e.g.,  f = N X 
NT  N 

X 25  KHz  since  the  sampling  interval  T = 40  p,sec  and  N ■=  1,  2,  3,  ....). 

Hie  integer  32  corresponds  to  the  number  of  samples  taken  following  each 
drive  pulse.  Thus  the  cyclic  repetition  interval  is  1280  p,sec  (32  X 40nsec) . 
Note  that  the  length  of  the  delay  line  was  chosen  so  that  successive  echoes 
are  about  40  p.sec  apart  in  order  that  the  sampling  pulse  overlap  the  echo 
duration.  We  then  choose  n memory  location  to  correspond  to  2n  radians. 

Since  n = Af/32T  = constant,  then  choose  n to  provide  sufficient  phase 
resolution  for  sine  wave  storage  consistent  with  memory  size.  At  the  same 
time,  since  the  frequency  of  the  synthesizer  earn  only  be  set  to  a multiple 
of  0.01  Hz,  then  A f must  be  a multiple  of  0.01  Hz.  The  possible  combinations 
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of  n and  Af  are  either  n =>  125  and  Af  = 6.25,  or  n = 625  and  Af  = 1.25. 

Since  n = 625  gives  the  better  phase  resolution,  this  latter  combination 
has  been  selected.  The  values  of  n larger  than  625  have  been  neglected 
even  though  there  are  solutions  for  n larger  than  625,  because  of  limited 
memory  size.  As  it  is,  we  require  625  locations  for  each  sinusoid,  namely 
each  echo,  to  be  stored. 

3.  Data  Acquisition  Software  outline 

The  memory  region  associated  with  data  is  initially  set  to  zero 
(memory  locations  from  21,000  to  37,777).  in  the  initialization  procedure, 
the  programmable  reset  pulse  must  be  sent  to  reset  the  R.P.  pulser  to  its 
initial  state.  This  was  accomplished  by  an  output  to  the  external  control 
register  "XA"  bit  (i.e.,  B = 2504  1004  with  A = 0000  3000  ).  Next,  the 
programmable  clock  "XC"  in  the  interface  module  is  reset  by  B = 2504  1004 
with  A = 0000  3000,  and  then  set  to  produce  a pulse  every  40  p,sec  by 
B = 2504  1004  with  A = 0014  2400.  The  purpose  of  this  pulse  is:  1)  to 
activate  the  repetition  rate  counter  in  the  R.F.  pulser,  2)  to  enable  a 
trigger  to  the  sampling  scope  every  40  p,sec,  3)  to  initiate  the  ADC  every 
40  p.sec  and  4)  to  send  a program  interrupt  to  computer  every  40  p-sec. 

The  same  instruction  that  sets  the  programmable  clock  also  sets  the  mode 
of  data  transfer  to  the  computer  memory  through  the  DMA,  direct  memory 
access  channel.  By  setting  bit  10=1,  we  have  selected  the  single  address 
mode,  so  that  new  output  from  ADC  is  transferred  directly  to  memory  loca- 
tion 20,000  every  clock  interval,  40  p,sec. 

This  means  that  every  40  p,sec  fresh  data  enter  memory  location 
20,000  and  the  instruction  at  interrupt  location  265  is  executed.  A program 
instruction  of  "Branch  to  read  data"  is  stored  in  the  interrupt  location 
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in  order  to  carry  out  the  sequence  of  reading  data  (i.e.,  1400/3027)  into 
the  sequential  locations  as  each  interrupt  comes  along.  The  first  inter- 
rupt is  in  synchronization  with  the  start  of  the  R.  F.  pulser,  as  is  each 
subsequent  32nd  interrupt. 

By  using  the  initial  interrupt  and  every  32nd  interrupt  or  after, 

as  a synchronization  method,  sequence  of  the  data  that  appears  at  memory 

location  20,000  is  numbered  b , b,  ....  b, ,,  b, b,,  ....  b_,  and  repeated. 

’ o’  1 . 14’  15’  16  31 

When  ADC  channel  1 is  selected  (i.e.,  B = 25041004  with  A = 00142400), 
then  the  data  at  memory  location  20,000  are  the  b (response)  signal;  where- 
as when  channel  2 of  the  ADC  is  selected  (i.e.,  B = 25041004  with  A = 

= 00142400),  then  the  drive  signal  A is  selected.  The  "Read  Data"  program 
is  designed  so  that  the  sequences  at  location  20, 000  correspond  to  bQ,  b^, 

b2  **•  bl4>  ai5>  bl6  b31  311(1  tb-*-s  repeated.  The  entire  "READ  Data" 

routine  is  recycled  456  times  in  order  to  preaverage  the  data  points  to 
reduce  the  noise  in  the  sine  wave  reconstruction  process. 

When  the  "Read  Data"  program  is  finished,  the  digital-to-analog 
converter  (DAC)  is  utilized  to  display  the  stored  data  sequentially  on  the 
CRT.  The  detailed  program  is  listed  in  Appendix  A. 

E.  Data  Analysis 

1.  Cross-Correlation  and  Averaging 

The  determination  of  the  relative  amplitude  and  phase  difference 
of  the  drive  and  response  signal  is  the  most  important  part  of  the  dynamic 
analysis  experiment.  As  indicated  earlier,  the  sinusoidal  drive  signal 
A = a^5  sin  out  is  generated  by  the  frequency  synthesizer  through  the  R. F. 
pulser.  The  response  signals  B^  = b^  sin  (out  + cp^)  + noise  from  the  out- 
put transducer  of  the  delay  line  are  sinusoidal  echoes  with  the  same 
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frequency  but  have  noise  components  in  addition.  The  subscript  i desig- 
nates the  echo  number;  the  first  response  signal  has  traveled  the  length 
of  the  delay  line  without  reflection,  but  it  will  be  still  called  the 
zeroth  echo  for  the  sake  of  convenience. 

The  method  of  walking  frequencies  enables  sampling  of  the  amp- 
litude of  each  analog  signal  in  real  time  at  intervals  of  40  (xsec.  The 
digitized  amplitude  of  the  sampled  signal  is  then  stored  into  the  data 
region  of  the  memory  (i.e,  locations  22200-  36600)  as  mentioned  in  "Read 
Data"  section.  For  each  individual  echo  and  the  drive  signal,  the  sampled 
and  digitized  values  are  stored  sequentially  according  to  their  phase  of 
2TT/n,  where  there  cure  n «*  625  locations  for  each  period  for  each  echo  and 
the  drive  signal.  Thus,  a single  sinusoid  is  reconstructed  in  memory  for 
each  echo;  the  data  cure  accumulated  in  625  X 32  X 40  p.sec. 

The  data  obtained  and  stored  in  this  manner  lends  itself  ideally 
to  the  determination  of  the  phase  difference  6 and  relative  amplitude  b^/a15 
through  the  mathematical  technique  of  cross-correlation. 

Since  the  signals  are  of  known  frequency  uu,  the  following  useful 
integrals  can  be  formed 


.VC 

) Cl^Sin  + C sih  u)t  cM  = 7i  C CoS  <j>^ 

'o 

( Sin  (uJt-t  <£)  c CosujtM  - 7lC  <\>h 

Jt> 

sin  Cwt  * ) c SI  n wt  dw t = 7C  c be 


(51) 


(52) 


(53) 


( be  Sin  ( 60 1 + C c«s  iot  du)t=7L  cbt’smfo  (54) 

•'ft 
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in  which  co  and  cp . are  each  phase  angles  of  a and  b . relative  to  an  arbi- 
a i i 

trarily  chosen  time  t = 0 in  Eqs.  (51)  - (54).  An  ideal  software  constructed 
is  stored  in  memory  and  is  used  as  a reference  signal.  This  helps  to 
eliminate  errors  due  to  nonzero  base  line  in  the  above  integrations. 
However,  the  phase  angles  cp  of  bi  relative  to  a is  just  cp.^  - cpQ,  thus 
independent  of  choice  for  t = 0. 

The  desired  experimental  quantities  b Va  and  Acp  can  be  obtained 
through  the  use  of  Eqs.  (51)  - (54).  This  technique  of  integration  is 
called  the  method  of  cross-correlation.  The  phase  bhift  of  rr/2  required 
for  Eqs.  (52)  and  (54)  is  achieved  by  a circular  rotation  of  the  b^  vector 
regions  through  N/4  locations,  where  N = 625  here.  Therefore,  n/4  becomes 
the  noninteger  number  of  156.25;  but  the  memory  locations  must  be  integer 


numbers. 


In  order  to  resolve  this  difficulty,  the  following  approxima- 


tions for  Eqs. (51)  and  (52)  were  adopted: 


j = 


a,ysih  Si  hurt  dwt  = 4-  2-  CLiCj 

(9  ° 


= 4 7L  aiSC  C0S<J>A 


r in 

4-)  ft, 5 Sin  (wt  + <f>A)  C CoS  wt  dwt 


2.7L 

~ 3 * aJ  Cj  + 4 


0=  624 

&jCj„S7 

6 ^ J-o 


~ 4 7L  a,jr  C.  sin 


Notice  that  the  data  region  is  considered  to  have  circular  boundary  con- 
ditions, so  that  j ■ 625,  626,  ....  is  equivalent  to  j ■ 0,  1, 

This  approximation  is  exact  to  the  extent  that  linear  inter- 
P°^ati°n  is  valid  for  the  value  of  the  cosine  function  within  an  interval 
of  2tt/625  radians.  The  same  procedure  can  be  applied  to  Eqs.  (53)  and  (54). 
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PART  III 


CALIBRATION 


file  calibration  of  the  device  is  based  on  fluids,  OOP  and  ETOH 
(ethanol),  which  we  believe  to  be  Newtonian  up  through  the  frequency  range 
of  interest,  i.e.,  beyond  30  MHz. 

For  a Newtonian  fluid  the  mechanical  resistance  and  reactance 
are  equal  and  given  by 


R=  X = 


r\ 

a 


(59) 


where  p and  7)  are  the  density  and  the  steady-state  viscosity,  respectively. 
The  relation  between  measured  attenuation  per  unit  length  AA  and  phase 
shift  per  unit  length,  Ab,  and  R and  X were  given  by  equation  and  shown 
here  A A - R 


A/\  = 
*0  = 


X 


(60) 


2 Vs  fib 

where  V , and  2b  are  the  wave  velocity  in  the  delay  medium,  the  density 
of  delay  material  and  the  thickness  of  delay  line,  respectively. 

The  raw  experimental  attenuation  and  phase  shift  data  for  ETOH 
and  DOP  is  seen  in  Figs. 16  and  17.  This  data  is  measured  at  the  various 
indicated  temperature  from  25° c to  10° C and  over  the  frequency  range  from 
3.4  MHz  to  27.1  NHz.  These  data  are  then  replotted  as  a function  of  the 
square  root  of  frequency  in  Figs. 18  through  20  in  order  to  test  the 
linearity  predicted  by  Eqs.  (59)  and  (60)  above. 

Notice  the  data  in  Fig. 18  in  which  the  ratio  of  output  to  input 
of  the  entire  delay  line  system  including  amplifiers  is  plotted  as  a 
function  of  JT . The  major  resonance  of  the  transducers  is  evident  as  is 
the  fact  that  the  total  attenuation  for  the  more  viscous  DOP  is  greater 
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than  for  EPOH.  An  important  secondary  effect  can  also  be  observed  from 
this  curve.  That  is,  in  addition  to  the  higher  attenuation  of  the  DOP 
response  curve,  it  can  be  seen  (by  vertical  and  horizontal  shifting)  that 
the  DOP  curve  is  shifted  slightly  towards  lower  frequency. 

Die  primary  effect  of  liquid  loading  on  the  delay  line  is  to 
cause  an  »/f  bias  on  the  magnitude  of  attenuation,  hence  some  apparent 
distortion  in  the  resonance  curve.  Beyond  this,  the  shift  in  the  DOP 
response  curve  is  attributed  to  the  coupling  of  transducer  response  to 
mechanical  impedance  increment  of  the  delay  line.  Diis  coupling  is 
neglected  in  the  analysis. 

Die  data  of  Figs. 16  and  17  are  then  replotted  as  a function  of 
at  fixed  temperatures  as  presented  in  Figs. 19  and  20.  Next,  we  calcu- 
late the  theoretically  predicted  slope  from  delay  line  dimension,  and 
these  are  drawn  onto  the  same  figures. 

Die  theoretical  calculations  sure  as  follows.  From  Eqs. (59)  and 
(60),  we  obtained  the  attenuation  equation: 

(*/«*  = jK  ( &•))  JZ—r  i «n> 

it  Vs  fs 

where  V is  the  acoustic  velocity  in  a fused  quartz  (3.8  X 10^  cm/sec), 

b 

2b  is  the  thickness  of  the  delay  line  (0.154  cm),  p is  the  density  of  a 

s 

fused  quartz  (2.2  g/cc),  and  l is  the  length  of  the  liquid  - delay  line 
interface  region  (6.7  cm). 


In  practice,  the  cell  with  ETOH  is  taken  as  a reference  state 
thus  Eq.  (61)  must  be  modified  when  applied  to  DOP,  and  becomes 


Similarly,  the  phase  difference  of  DOP  with  respect  to  E70H  is 


CAB)i-  3 I 


(63) 


The  coefficient  ^T/2Vgbps  for  a fused  quartz  delay  material  is  1.39  X 10  5. 

The  slope  of  the  calibration  curves  can  be  obtained  from  Eqs. (62) 
and  (63)  together  with  steady-state  viscosity  and  density  data  for  the 
fluids.  These  are  summarized  in  Table  1 and  the  theoretically  predicted 
lines  are  indicated  in  Figs. 19  and  20. 


TABLE  1 

THEORETICAL  CALIBRATION  COEFFICIENTS 


r 


Conclusions  about  calibration.  As  seen  from  Figs. 18  to  20, 
the  experimental  calibration  data  are  consistent  with  the  straight  line 
predicted  theoretically  from  first  principles.  Some  sources  for  the  re- 
maining discrepancy  that  have  not  been  incorporated  into  the  theoretical 
analysis  are  as  follow: 

1)  The  coupling  between  the  transducer  and  the  liquid  via  the 
delay  line  has  not  been  included  in  the  analysis.  The  liquid 
changes  the  mechanical  impedance  of  the  delay  line  and  in 
turn  affects  the  load  seen  by  the  transducer,  and  hence 

its  response. 

2)  It  was  assumed  that  there  is  no  propagation  loss  in  the 
delay  line  in  the  absence  of  liquid  on  the  surface.  That 
is,  the  propagation  constant  m^  in  Eq.  (13)  is  purely 
imaginary.  However,  there  exists  a finite  propagation  loss 
in  any  material  due  to  the  internal  friction.  Therefore, 
the  propagation  vector  in  Eq. (13)  becomes  a complex 
number;  and  in  turn  both  the  phase  AB  and  the  attenuation 
AA  in  Eq. (28)  will  be  coupled  with  this  propagation  loss 
through  n^  of  Eq. (20). 

3)  Modes  of  propagation  other  than  SS  (O)  are  permissible  in- 
asmuch as  we  are  above  the  cut-off  frequency  of  2.2  MHz. 

The  magnitude  of  such  possible  modes  was  observed  to  be 
small.  The  echo  amplitude  as  seen  on  the  oscilloscope  only 
oscillates  by  a small  amount  as  we  tune  through  a continuous 
range  of  frequency  in  the  neighborhood  of  each  operating 
frequency. 

115 


r 


4)  The  smoothness  of  the  surface  can  be  of  critical  importance 


for  low  viscosity  liquids  such  as  ethanol  since  the  Shear 
wavelength  X - 2.5  X 10  3 cm  begins  to  approach  the 

O 

magnitude  of  the  surface  irregularities . 

In  practice,  since  the  experimental  calibration  data  (Table  3) 
does  show  fluctuations  from  the  theoretical  curve,  it  was  decided  to  use 
a separate  proportionality  constant  for  each  frequency.  These  constants 
cure  taken  as  the  average  of  the  measured  calibration  values.  Thus,  for 

R =•  K,  (f ) Afl,  (64! 

X = Ki  Cj ) A 8 <«> 

The  proportionality  constants  K^ff)  and  ^(f)  were  obtained  from  DOP  and 
ETOH  at  each  frequency  and  are  shown  in  Table  2. 


TABLE  2 

EXPERIMENTAL  CALIBRATION  COEFFICIENTS 


Table  3.  Experimental  (Raw)  Data 


8232.083 

- C.R. 

ETOH(25*C) 

(40  dB  preamp) 

b°As 

3.4  MHz 

263.728 

3.250 

263.627 

3.250 

263.526 

3.250 

6.0  MHz 

263.070 

5.142 

263.442 

5.156 

263.813 

5.170 

8.4  MHz 

321.181 

6.427 

321.195 

6.435 

321.210 

6.442 

(20  dB  preamp) 

11.2  MHz 

325.776 

0.9578 

325.736 

0.9611 

325.696 

0.9644 

14.5  MHz 

262.165 

1.605 

262.245 

1.605 

262.325 

1.605 

14.6  MHz 

248.012 

1.643 

248.237 

1.638 

248.463 

1.633 

19.5  MHz 

294.359 

2.071 

294.447 

2.071 

294.535 

2.070 

25.0  MHz 

276.835 

0.6500 

277.235 

0.6531 

277.634 

0.6561 

27.1  MHz 

5.835 

0.3792 

5.842 

0.3786 

5.849 

0.3780 

8228.112 
8228.138  - C.R. 


ETOH ( 20°  C ) 
(40  dB  preamp) 


3.4  MHz 

257.365 

3.267 

257.044 

3.272 

256.722 

3.276 

6.0  MHz 

251.096 

5.289 

251.060 

5.288 

251.023 

5.286 

8.4  MHz 

304.365 

6.546 

304.762 

6.536 

305.158 

6.527 

(20  dB  preamp) 

11.2  MHz 

302.003 

0.9873 

302.038 

0.9874 

302.072 

0.9874 

14.5  MHz 

230.574 

1.642 

230.930 

1.643 

231.285 

1.643 

14.6  MHz 

216.187 

1.692 

216.329 

1.685 

216.471 

1.677 

19.5  MHz 

252.530 

2.074 

252.531 

2.074 

252.531 

2.074 

25.0  MHz 

224.650 

0.6550 

224.650 

0.6558 

224.646 

0.6565 

27.1  MHz 

-51.689 

0.3788 

-51.740 

0.3802 

-51.764 

0.3816 

(308.26) 
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' - V 


'I 


ET0H(15*C) 
(40  dB  preamp) 

8223.475 

-C.R, 

b%n' 

av 

h;. 

3.4  MHz 

249.567 

3.279 

249.722 

3.290 

249.876 

3.300 

6.0  MHz 

236.890 

5.383 

236.816 

5.373 

236.742 

5.362 

8.4  MHz 

283.667 

284.566 

6.715 

6.680 

284.117 

6.698 

(20  d3  preamp) 

11.2  MHz 

274.634 

1.008 

274.652 

1.008 

274.669 

1.008 

14.5  KHz 

193.780 

1.689 

193.934 

1.686 

194.088 

1.682 

14.6  MHz 

179.523 

1.710 

179.580 

1.708 

179.635 

1.706 

19.5  MHz 

202.150 

2.081 

202.127 

2.082 

202.104 

2.082 

25.0  MHz 

161.295 

0.6577 

161.290 

0.6554 

161.284 

0.6530 

27.1  KHz 

-119.982 

0.3829 

-120.072 

0.3839 

-120.162 

0.3848 

(239.928) 

8219.156 

-C  R 

ETOH(10*C) 
(40  dB  prearar) 

242.156 

3.4  MHz 

3.301 

242.251 

3.304 

242.346 

3.306 

6,0  MHz 

222.972 

5.482 

223.040 

5.473 

223.108 

5.464 

8.4  MHz 

26/a.  368 

264.243 

6.756 

264.305 

6.785 

(20  dB  preamp) 

6.813 

11.2  MHz 

248.230 

1.023 

248.276 

1.022 

14.5  MHz 

248.321 

1.021 

159.302 

1.720 

159.175 

1.719 

14.6  MHz 

159.047 

1.717 

143.679 

1.740 

143.570 

1.735 

143.461 

1,730 

19.5  MHz 

155.032 

2.067 

155.016 

2.068 

155.000 

2.069 

25.0  MHz 

102.450 

0.6539 

102.429 

0.6536 

102.408 

0.6533 

27.1  MHz 

-183.301 

0.3875 

-183.352 

0.3875 

-183.402 

0.3875 

(176.648) 

118 


8232.019  - 

■C.  R. 

\r 

DOP(25*C) 

Vo'?* 

h/a.!,' 

'Ocp 

Af\ 

(40  dB  preamp) 

3.4  MHz 

259.791 

3.003 

259.753 

3.011 

259.715 

3.019 

(3.874) 

(0.926) 

6.0  MHz 

258.489 

4.731 

258.498 

4.741 

258.507 

4.750 

(4.944) 

(0.919) 

8.4  MHz 

314.752 

6.187 

314.828 

6.166 

314.903 

6.144 

(6.367) 

(0.958) 

(20  dB  preamp) 

11.2  MHz 

318.957 

0.8560 

319.015 

0.8558 

319.072 

0.8555 

(6.721) 

(0.890) 

14.5  MHz 

252.650 

1.380 

252.442 

1.383 

252.234 

1.386 

(9.803) 

(0.862) 

14.6  MHz 

237.996 

1.393 

238.188 

1.397 

238. 380 

1.401 

(10.049) 

(0.853) 

19.5  MHz 

283.312 

1.644 

283.238 

1.645 

283.164 

1.645 

(11.209) 

(0.7943) 

25.0  MHz 

266.524 

0.5049 

266.496 

0.5041 

266.467 

0.5032 

(10.743) 

(0.772) 

27.1  MHz 

355.322 

0.2859 

355.912 

0.2850 

356,502 

0.2840 

(9.93) 

(0.753) 

DOP( 20CC ) 

8227.830 

-C.  R. 

(40  dD  prearao) 

3.4  MHz 

252.270 

3.026 

252. 182 

3.026 

252.094 

3.026 

(4.862) 

(0.9248) 

6.0  MHz 

243.911 

4.781 

244.173 

4.791 

244.435 

4.800 

(6.88?) 

(0.9060) 

8.4  MHz 

296.109 

6.003 

296.352 

5.985 

296.595 

5.967 

(8.410) 

(0.9157) 

(20  dB  Dreamp) 

11.2  MHz 

293.040 

0.8686 

293.045 

0.8692 

293.049 

0.8698 

(8.993) 

(0.8803) 

14. 5 MHz 

219.343 

1.413 

219.215 

1.413 

219.087 

1.413 

(11.715) 

( 0.8600) 

14.6  MHz 

204.759 

1.422 

204.797 

1.423 

204.834 

1.424 

(11.532) 

(0.8445) 

19.5  MHz 

238.305 

1.612 

238.333 

1.606 

238.360 

1.600 

(14.198) 

(0.7743) 

25.0  MHz 

211.294 

0.4950 

211.221 

0.4961 

211.148 

0.4971 

(13.428) 

(0.7564) 

27.1  MHz 

-65.547 

0.2785 

-65.265 

0.2792 

-64.983 

0.2798 

(13.528) 

(0.7344) 
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8223.496 

-C.  R. 

6k  " kk 

DOP(15*C) 

(40  dB 

preamp) 

3.4 

MHz 

244.787 

2.989 

244.812 

244.836 

2.999 

(4.91) 

6.0 

MHz 

230.023 

4.803 

230.075 

230.127 

4.837 

(6.741) 

8.4 

MHz 

276.643 

5.992 

276.943 

277.243 

6.164 

(7.174) 

(20  dB 

preamp) 

266.472 

11.2 

MHz 

0.8642 

266.451 

266.429 

0.8603 

(8.201) 

14.5 

MHz 

183. 032 

1.409 

183.048 

183.064 

1.405 

(10.886) 

14.6 

MHz 

168.380 

1.417 

168.402 

168.424 

1.410 

(11.178) 

19.5 

MHz 

190.120 

1.568 

190.178 

190.236 

1.569 

(11.949) 

25.0 

MHz 

151.889 

0.4762 

152.130 

152.370 

0.4780 

(9.16) 

27.1 

MHz 

-130.563 

0.2724 

-130.342 

-130.120 

0.2721 

(10.27) 

DOPI 

IlO'C) 

8219.114 

-C.  B. 

(40  d3 

preamp) 

236.811 

3.4 

MHz 

2.993 

236.873 

236.934 

2.997 

(5.738) 

6.0 

MHz 

215.229 

4.828 

215.483 

2.5.737 

4.834 

(7.557) 

8.4 

KHz 

254.678 

6.056 

254.655 

254.632 

6.026 

(9.650) 

(20  dB 

MHz) 

11.2 

MHz 

238.802 

0.8576 

238.776 

238.750 

0.8582 

(9.500) 

14.5 

MHz 

147.014 

1.393 

147.162 

14.6 

147.310 

1.399 

(12.013) 

MHz 

132.566 

1.398 

132.558 

132.551 

1.400 

(11.012) 

19.5 

MHz 

141.995 

1.514 

142.019 

142.043 

1.506 

(12.997) 

25.0 

MHz 

90.726 

0.4524 

90.662 

90.598 

0.4542 

(11.767) 

27.1 

MHz 

-195.499 

0.2629 

-195.193 

-194.887 

0.2625 

(11.841) 
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l I 


Ra-t.'o*.  -a/\ 

0.2994 

(0.9100) 

4.820 

(0.8970) 

6.078 

(0.9074) 

0.8623 

(0.8554) 

1.407 

(0.8345) 

1.414 

(0.8279) 

1.569 

(0.7536) 

0.4771 

(0.7280) 

0.2723 

(0.7093) 


2.995 

(0.9065) 

4.831 

(0.8827) 

6.041 

(0.8903) 

0.8579 

(0.8394) 

1.396 

(0.8121) 

1.399 

(O.8063) 

1.510 

(0.7301) 

0.4533 

(0.6935) 

0.2627 

(0.6779) 


: 


1 


8232.041  -C. 

5%  PS/DOP(25#C)  , 1,/ 

(40  dB  preamp)  "/f 

(byin'jaKr 

3.4 

MHz 

259.151 

2.963 

259.116 

2.963 

259.080 

2.962 

(4.511) 

(.9117) 

6.0 

MHz 

257.141 

4.621 

257.148 

4.613 

257.154 

4.604 

(6.294) 

(.8947) 

8.4 

MHz 

313.615 

5.766 

314.182 

5.763 

(20  dB 

314.749 

5.760 

(7.013) 

(.8956) 

preamp) 

11.2 

MHz 

317.642 

0.8482 

317.746 

0.8454 

14.5 

317.850 

0.8426 

(7.99) 

(.8796) 

MHz 

250.229 

1.392 

250.485 

1.390 

250.740 

1.388 

(11.76) 

( . 8660) 

14.6 

MHz 

236.202 

1.405 

236.271 

1.407 

236.340 

1.410 

(11.966) 

(.8590) 

19.5 

MHz 

280.097 

1.522 

280.110 

1.521 

280.122 

1.520 

(14.337) 

(.7344) 

25.0 

MHz 

266.824 

O.4532 

266.802 

0.4529 

266.779 

0.4526 

(10.433) 

(.6935) 

27.1 

MHz 

356.560 

0.2543 

356.516 

0.2533 

356.472 

0.2522 

(9.326) 

(.6690) 

52  PS/D0P( 20* C ) 8227. 

895  -C.  R. 

(40  dB 

preamp) 

3.4 

KHz 

252.190 

2.939 

251.941 

2.934 

251.691 

2.929 

(5.103) 

( .8967) 

6.0 

MHz 

2^3.782 

4.682 

243.792 

4.681 

243.802 

4.680 

(7.268) 

(.8852) 

8.4 

MHz 

295.115 

5.941 

295.044 

5.955 

294.973 

5.968 

(9.718) 

(.9111) 

(20  cLj 

preamp) 

292.216 

11.2 

MHz 

0.8417 

292.305 

0.8416 

14.5 

292.393 

0.8414 

(9.733) 

(.8523) 

MHz 

216.303 

1.385 

216.495 

1.388 

14.6 

216.687 

1.390 

(14.435) 

(.8448) 

MHz 

201.533 

1.396 

201.872 

1. 3Q0 

202.211 

1.384 

(14.457) 

(.8249) 

19.5 

MHz 

234.584 

1.465 

234.535 

1.467 

234.485 

1.468 

(17.996) 

(.7073) 

25.0 

MHz 

210.145 

O.4343 

210.187 

0.4361 

210.229 

0.4378 

(14.463) 

(.6649) 

27.1 

KHz 

-64.954 

0.2467 

-65.335 

0.2445 

-65.715 

0.2422 

(13.595 

(.6430) 

10* 

PS/DOP( 25*C ) 

8232.044  -C. 

. / 

(40  dB 

preamp) 

4*-f* 

feJtv  - 

3. '4 

MHz 

258.523 

2.883 

258.149 

2.882 

257.774 

2.880 

(5.478) 

(.8868) 

6.0 

MHz 

255.795 

4.535 

255.726 

4.538 

255.656 

4.541 

(7.716) 

( .8801) 

8.4 

MHz 

311.984 

5.708 

312.294 

5.704 

312.603 

5.700 

(81901) 

(.8864) 

(20  dB 

preamp) 

11.2 

MHz 

315.837 

0.8119 

316.019 

0.8119 

316.201 

0.8118 

(9.717) 

(.8445) 

14.5 

MHz 

246.896 

1.317 

246.835 

1.319 

246.774 

1.322 

(15.410) 

( .8218) 

14.6 

MHz 

232.734 

1.328 

232.779 

1.329 

232.824 

1.329 

(15.458) 

(.8113) 

19.5 

MHz 

275.604 

1.376 

275.660 

1.374 

275.716 

1.372 

(18.787) 

(.6634) 

25.0 

MHz 

264.323 

0.3976 

264.384 

0.3958 

264.444 

0.3939 

(12.851) 

( .6060) 

27.1 

MHz 

355.595 

0.2181 

355.388 

0.2180 

355.180 

0.2178 

(10.454) 

(.5758) 

1055 

PS/DOP( 

20»c)  8227.836  -C.  R 

• 

(40  dB 

oreamp) 

3.4 

MHz 

250.952 

2.927 

251.230 

2.906 

251.507 

2.885 

(5.814) 

( .8881) 

6.0 

MHz 

243.077 

4.620 

242.974 

4.610 

242.871 

4.600 

(8.086) 

(.8717) 

8.4 

MHz 

295.461 

5.777 

295.908 

5.760 

296.355 

5.743 

(8.854) 

(.8812) 

(20  dB 

preamp) 

11.2 

MHz 

289.339 

0.8121 

289.436 

0.8141 

289.533 

0.8160 

(12.602) 

( .8245) 

14.5 

MHz 

212.131 

1.3H 

212.081 

1.315 

212.031 

1.318 

(18.849) 

( .8003) 

14.6 

MHz 

197.827 

1.325 

198.037 

1.329 

198.247 

1.333 

(18.292) 

(.7887) 

19.5 

MHz 

229.125 

1.332 

229.164 

1.331 

229.202 

1.329 

(23.367) 

(.6417) 

25.0 

MHz 

206.771 

0.3850 

206.554 

0.3845 

206.336 

0.3840 

(18.096) 

(.5863) 

27.1 

KHz 

-68.020 

0.2104 

-68.314 

0.2095 

-68.608 

0.2086 

(16.574) 

(.5510) 

10* 
(40  dB 

8223.462  -C.  R. 

PS /OOP  (15*0  , l / . 

preamp) 

(4 AT 

((>•/., r)  £*■ 

3.4 

MHz 

243.202 

2.835 

243.172 

2.837 

243.141 

2.839 

(6.550) 

(.8623) 

6.0 

MHz 

227.821 

4.564 

227.477 

4.563 

227.132 

4.562 

(9.339) 

(.8492) 

8.4 

MHz 

272.586 

5.544 

272.85 

5.573 

273.114 

5.601 

(11.267) 

(.8320) 

(20  dB 

preamp) 

262.035 

(12.617) 

11.2 

MHz 

261.952 

0.7985 

0.7989 

262.117 

0.7993 

(.7926) 

14.5 

MHz 

176.579 

1.291 

176.588 

1.291 

176.596 

1.290 

(17.346) 

(.7657) 

14.6 

MHz 

161.968 

1.296 

162.352 

1.296 

162.736 

1.295 

(17.228) 

(.7588) 

19.5 

MHz 

181.210 

1.273 

181.161 

1.273 

181.111 

1.272 

(20.966) 

(.6114) 

25.0 

MHz 

147.093 

0.3600 

147.105 

0.3608 

147.117 

0.3615 

(14.185) 

(.5505) 

2?.l 

MHz 

226.585 

0.2043 

226.635 

0.2029 

226,684 

0.2014 

(13.293) 

(.5285) 

10* 

P3/COP ( 

10*0  8219 

.i4o  -c.  n. 

(40  dB 

pref«mp) 

234.272 

3.4 

MHz 

234.339 

2.757 

2.755 

234.245 

2.753 

(7.959) 

(.8338) 

6.C 

MHz 

212.682 

4.401 

212.688 

4.401 

212.693 

4,401 

(10.352) 

(.8041) 

8.4 

MHz 

251.973 

5.607 

251.848 

5.581 

251.585 

5.582 

(12.457) 

(.8225) 

251.986 

5.553 

(20  dB 

preamp) 

0.7696 

11.2 

MHz 

235.027 

0.7700 

235.032 

235.037 

0.7693 

(13.244) 

(.7530) 

14.5 

MHz 

141.557 

1.238 

141.558 

1.238 

141.559 

I.238 

(17.617) 

(.7202) 

14.6 

MHz 

127.069 

1.248 

126.957 

1.246 

126.844 

1.243 

(16.613) 

(.7181) 

19.5 

MHz 

133.258 

1.194 

133.306 

1.193 

133.354 

1.192 

(21.71) 

(.5768) 

25.0 

MHz 

86.866 

0.3322 

86.709 

0.3319 

86.551 

0.3316 

(15.72) 

(.5078) 

27.1 

MHz 

-197.774 

0.1937 

197.976 

0.1942 

-198.177 

0.1946 

(14.624) 

(.5011) 
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8232.070  - C.R 


20%  PS/DOP (25*C)  . 

)av. 

( bo/ais), 

(40  dB  preamp) 

& A 

3 . 4 MHz 

255.338 

2.674 

255.633 

2.665 

255.927 

2.656 

(7.994) 

(.820) 

6.0  MHz 

252.842 

4.145 

252.978 

4.144 

253.113 

4.142 

(10.462) 

(.8037) 

8.4  MHz 

309.481 

5.173 

309.206 

5.134 

308.934 

5.095 

(11.989) 

(.7978) 

(20  dB  preamp) 

11.2  MHz 

312.380 

0.7168 

312.267 

0.7164 

312.153 

0.7159 

(13.469) 

(.7453) 

14.5  MHz 

242.907 

1.153 

242.734 

1.155 

242.561 

1.157 

(19.511) 

(.7196) 

14.6  MHz 

228.377 

1.158 

228.427 

1.158 

228.476 

1.157 

(19.810) 

(.7070) 

19.5  MHz 

269.758 

1.108 

269.665 

1.107 

269.572 

1.105 

(24.782) 

(.5345) 

25.0  MHz 

260.455 

0.307 

260.317 

0.3068 

260.179 

0.306 

(16.918) 

(.470) 

27.1  MHz 

351.724 

0.166 

351.830 

0.166 

351.935 

0.166 

(14.012) 

(.439) 

8227 

.867  - C.R. 

20%  PS/DOP (20* 

c) 

(40  dB  preamp) 

3.4  MHz 

247.842 

2.612 

247.622 

2.617 

247.441 

2.622 

(9.422) 

(.7708) 

6.0  MHz 

238.075 

4.065 

238.025 

4.076 

237.975 

4.086 

(13.035) 

(.7708) 

8.4  MHz 

287.809 

5.121 

288.142 

5.086 

288.475 

5.050 

(16.620) 

(.7021) 

(20  dB  preamp) 

11.2  MHz 

285.769 

0.6935 

285.718 

0.693 

285.667 

0.693 

(16.320) 

(.702) 

14.5  MHz 

208.294 

1.112 

208.283 

1.111 

208.272 

1.109 

(22.647) 

( .676) 

14.6  MHz 

193,663 

1.118 

193.502 

1.122 

193.340 

1.125 

(22.827) 

(.666) 

19.5  MHz 

223.383 

1.045 

223.483 

1.047 

223.583 

1.048 

(29.048) 

(.505) 

25.0  MHz 

203.875 

0.284 

203.909 

0.283 

203.942 

0.283 

(20.741) 

(.4321) 

27.1  MHz 

288.299 

0.157 

288.865 

0.155 

289.153 

0.153 

(19.460) 

(0.407) 
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20* 

PS/DOP( 15*C ) . 8223»  525 

-C.  R. 

a . 

(40  dB 

preamp) 

ft -^-4 

^#2..V0 

0 

3.4 

MHZ 

238.560 

2.530 

238.790 

2.528 

239.020 

2.525 

(12.906 

(.7684) 

6.0 

MHz 

223.028 

3.934 

223.106 

3.935 

223.183 

3.935 

(13.710) 

(.7324) 

8.4 

MHz 

266.773 

4.856 

266.737 

4.854 

266.701 

4.852 

(17.380) 

(.724?) 

(20  dB 

preamp) 

11.2 

MHz 

257.258 

0.6557 

257.620 

0.6595 

257.982 

0.6633 

(17.032) 

(.6543) 

14.5 

MHz 

172. 560 

1.043 

172.38? 

1.041 

172.213 

1.039 

(21.54?) 

(.6174) 

14.6 

MHz 

157.992 

1.049 

158.016 

1.049 

158.040 

1.049 

(21.564) 

(.6141) 

19.5 

MHz 

175.617 

0.9533 

175.677 

0.9539 

175.736 

0.9544 

(26.45) 

(.4566) 

25.0 

MHz 

145.024 

0.2598 

145.166 

0.2602 

145.307 

0.2606 

(16.124) 

(.3970) 

27.1 

MHz 

224.689 

0.1410 

224.671 

0.1419 

224.653 

0.1428 

(15.072) 

(.3696) 

20* 

PS/DQP( 10* 

c)  8219’ 

150  -C.  R. 

(40  dB 

preamp) 

3-4 

MHz 

229.769 

2.264 

229.941 

2.268 

229.912 

2.271 

(12.310) 

(.6864) 

6.0 

MHz 

208.619 

3.505 

208.666 

3.508 

208.713 

3.5H 

(14.374) 

(.6410) 

8.4 

KHz 

247.315 

4.344 

247.219 

4.379 

247.122 

4.414 

(17.086) 

(.6454) 

(20  dB 

preamp) 

11.2 

MHz 

230.958 

O.585O 

230.950 

O.5816 

230.941 

0.5783 

(17.326) 

(0.5691) 

14.5 

MHz 

138.167 

0.9080 

137.934 

0.9073 

137.700 

0.9065 

(21.241) 

(.5278) 

14.6 

MHz 

123.288 

0.9099 

123.275 

0.9114 

123.261 

0,9129 

(20.295) 

(.5253) 

19.5 

MHz 

129.345 

0.8291 

129.357 

0.8285 

129.368 

0.8278 

(25.659) 

( .4006) 

25.0 

KHz 

85.422 

0.2159 

85.925 

0.2170 

86.459 

0.2178 

16.504 

(.3320) 

85.895 

0.2174 

27.1 

MHz 

-194.935 

0.1234 

-195.079 

.1231 

-194.504 

0.1246 

(11.727) 

(.3176) 

-195.797 

0.1212 
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40* 
(40  dB 

, / # % 8232.105  -c.  R. 

PS/DOP(25*C)  / / 

preamp) 

(i-l  4- 

RtCtjo  e. 

3.4 

MHz 

246.189 

1.871 

246.385 

1.878 

246.580 

1.884 

(17.242) 

(.5??K) 

6.0 

MHz 

243.488 

2.708 

243.192 

2.712 

242.895 

2.715 

(20.250) 

(.5260) 

8.4 

MHz 

297.284 

3.076 

296.844 

3.066 

296.403 

3.056 

(24.35D 

(.4765) 

(20  dB 

preamp) 

304.216 

11.2 

MHz 

0.4125 

304.273 

0. 412*5 

304.330 

0.4124 

(21.463) 

( .4292) 

14.5 

MHz 

237.339 

0.6173 

237.054 

0.6175 

236.769 

0.6177 

(25.191) 

(.3847) 

14.6 

MHz 

223.373 

0.6138 

223.607 

0.6175 

223.840 

0.6l60 

(24.630) 

(.3754) 

19.5 

MHz 

268.482 

0.5369 

268.504 

0.5374 

268.525 

0.5278 

(25.943) 

(.2595) 

25.0 

MHz 

266.099 

0.1327 

266.385 

0.1371 

266.673 

0.1415 

(10.850) 

(0.2099) 

27.1 

MHz 

O.863 

0,0701 

361.630 

0.0685 

2.430 

0.0670 

( 4.212) 

(.1809) 

1.596 

0.0683 

40* 

PS/DOP( 

20®C ) 8227. 

877  -C.R. 

(40  dB 

prear.it>) 

236.617 

3.4 

MHz 

1.746 

236.497 

1.744 

236.376 

1.741 

(20.547) 

(.5330) 

6.0 

MHz 

227.231 

2.464 

227.303 

2.468 

227.375 

2.451 

(23.757) 

(.4667) 

8.4 

MHz 

277.569 

2.753 

277.738 

2.842 

277.906 

2.905 

(27.024) 

(.4348) 

(20  dB 

preamp) 

278.615 

- 

11.2 

MHz 

O.3615 

278.834 

0.3630 

279.053 

0.3644 

(23.204) 

(.3676) 

14.5 

MHz 

204.239 

O.5349 

204.465 

0.5341 

204.689 

0.5332 

(26.465) 

(.325D 

14.6 

MHz 

191.209 

0.5342 

191.149 

0.5367 

191.089 

0.5392 

(25.180) 

(.3185) 

19.5 

MHz 

225.697 

0.4568 

225.689 

0.4562 

225.681 

O.4556 

(26.842) 

(.2200) 

25.0 

MHz 

215.101 

0,1143 

215.375 

0.1146 

215.649 

0.1148 

(9.275) 

(.174?) 

27.1 

MHz 

305.359 

0.0596 

304.365 

O.O585 

303.075 

0.0591 

(3.895) 

(.1539) 

304.661 

0.0567 
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PART  IV 


EXPERIMENTAL  RESULTS 

The  series  of  experimental  results  reported  herein  form  the  basis 
of  a study  of  viscoelastic  properties  of  low  molecular  weight  polymer  solu- 
tions based  on  polystyrene  in  di-actyl-phthalate  plasticizer.  The  study 
encompassed  a wide  concentration  range  from  5%  to  40%  of  2200  molecular 
weight  polystyrene  over  a temperature  range  from  10° C to  25° C. 

A.  Dynamic  Measurement  of  Solvent  and  Solution 

In  this  kind  of  experiment  it  is  important  to  take  into  account 
the  history  of  sample  handling  in  order  to  minimize  the  effect  of  long 
term  drifts.  It  is  therefore  useful  to  specify  the  complete  chronology. 

The  measurements  were  done  in  the  following  temperature  sequence  .- 
25°C,  20°C,  15° C and  10° C. 

The  cell  was  cleaned  with  DOP,  followed  by  a rinse  with  ETOH. 

After  each  sample  change  or  new  temperature  setting,  it  is 
necessary  to  wait  at  least  two  hours  before  taking  data,  in  order  for  the 
liquid  to  reach  thermal  equilibrium.  Emptying  and  filling  were  always 
carried  out  at  the  temperature  of  25° C. 

The  amplifier  at  the  output  of  the  delay  line  was  set  to  40  db 
gain  for  frequencies  of  3.4  MHz,  6.0  MHz,  8.4  MHz  and  to  20  db  gain  for 
11.2  MHz,  14.5  J«z,  14.6  MHz,  19.5  MHz,  25.0  MHz  and  27.1  MHz.  These  set- 
tings were  chosen  so  as  not  to  overload  the  sample  and  hold  circuit. 

The  chronology  is  summarized  in  Table  4. 
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TABLE  4 


SAMPLE  CHRONOLOGY 


Date 

Liquid 

Temperature 

Counter  Reading 

5/14/76 

ETOH 

25°  C 

8232.083 

5/14/76 

ETOH 

20°C 

8228.120 

5/14/76 

ETOH 

15°  C 

8223.470 

5/14/76 

ETOH 

10°  C 

8219.160 

5/15/76 

DOP 

25°  C 

8232.020 

5/15/76 

DOP 

20°  C 

8227.830 

5/14/76 

DOP 

15°  C 

8223.488 

5/16/76 

DOP 

10°  C 

8219.120 

5/17/76 

5%  PS /DOP 

25°C 

8232.036 

5/17/76 

5%  PS/DOP 

20°  C 

8227.890 

5/17/76 

5%  PS/DOP 

15°  C 

8223.494 

5/17/76 

5%  PS/DOP 

10°  C 

8219.175 

5/18/76 

10%  PS/DOP 

25°  C 

8232.045 

5/18/76 

10%  PS/DOP 

20°  C 

8227.836 

5/19/76 

10%  PS/DOP 

15°  C 

8223.462 

5/19/76 

10%  PS/DOP 

10°  C 

8219.140 

5/19/76 

20%  PS/DOP 

25°  C 

8232.075 

5/20/76 

20%  PS/DOP 

20°  C 

8227.867 

5/20/76 

20%  PS/DOP 

15°  C 

8223.525 

5/20/76 

20%  PS/DOP 

10°  C 

8219.150 

5/21/76 

40%  PS/DOP 

25°  C 

8232.100 

5/21/76 

40%  PS/DOP 

20°  C 

8227.875 

5/21/76 

40%  PS/DOP 

15°  C 

8223.504 

5/21/76 

40%  PS/DOP 

______  __ 

10°  C 

8219.120 

Table  5 


Ethanol 

Temp. 

C c ) 


Summary  of  Data 


(degree) 


Amplitude 

Ratio 

hQ^lS 


Calc. 

Newt. 

XRef 


f = 3.4  MHz 


25 

263.4 

3.26 

304 

20 

256.6 

3.27 

317 

15 

249.5 

3.29 

337 

10 

242.2 

3.30 

352 

f =6.0 

MHz 

25 

263.1 

5.17 

404 

20 

250.0 

5.28 

421 

15 

236.6 

5.37 

448 

10 

222.5 

5.48 

468 

f a 11. 

2 MHz 

25 

325.3 

9.62 

552 

20 

300.1 

9.90 

576 

15 

274.3 

10.07 

612 

10 

247.4 

10.24 

639 

f = 14. 

5 MHz 

25 

261.7 

16.07 

628 

20 

228.3 

16.46 

655 

15 

193.9 

16.85 

697 

10 

158.3 

17.21 

727 
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Temp. 

1"°) 

(degree) 

Amplitude 

Ratio 

V«15 

Calc. 

Newt. 

Ref 

£s  14.6 

MHz 

25 

247.5 

16.38 

630 

20 

213.6 

16.89 

657 

15 

178.8 

16.89 

699 

10 

142.4 

17.37 

730 

f * 19.5  MHz 


25 

293.5 

20.71 

729 

20 

248.3 

20.77 

760 

.15 

202.0 

20.81 

808 

10 

153.4 

20.69 

843 

fa  25.0 

MHz 

25 

277.0 

6.55 

825 

20 

219.8 

6.55 

860 

15 

161.3 

6.55 

915 

10 

100.2 

6.55 

955 

f = 27.1 

MHz 

25 

364.8 

3.77 

859 

20 

303.2 

3.80 

895 

15 

240.0 

3.84 

953 

10 

174.0 

3.88 

994 
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Temp.  ACf  Attenuation  Calculated  Measured  Measured  RX/X 

(°C)  (degree)  (nepers)  Newtonian  R X 


065  44262  25283  10362 


cn  <n  r*  vo  oo  't  id  cn  id 

O'  (0Q5  rOro^M1  <H  Tj<  (II  tt  <( 

CN  9) 


ID  j 

CN 

CN 

M 

105 


TJ  C 

0)  (0 
•U  M 

<0  c 

ro 

*<* 

vO 

00 

rH 

CM 

00 

r* 

O' 

H 0 

ro 

<T\ 

in 

in 

ro 

H 

rH 

Stl^ 

04 

in 

o 

VD 

o 

r* 

• r* 

CM 

o> 

ro 

ro 

CM 

Pm  ro 

rH  (I) 


10.6  -.214  5689  5460  5229 


Attenuation  Calculated  Measured  Measured 


250  6495  6107  5904 


L A ft 

^<P  Attenuation  Calculated  Measured  Measured  RX/^ 

(degree)  (nepers)  Newtonian  R X 


530  69045  33637  7189 


17.7  -.421  7506  7560  7215 


Attenuation  Calculated  Measured  Measured 


895  79794  31081  5667 


At  14.5  MHz  X=  358 

H=  22,370 


Temp. 

(°C) 

2^4 

(degree) 

Attenuation 

(nepers) 

Calculated 

Newtonian 

XL 

DOP-w.r , 

.t.  ETOH 

25 

9.2 

-.136 

3678 

20 

10.1 

-.153 

4089 

15 

10.7 

-.181 

4603 

10 

11.7 

-.210 

5254 

NBS  OIL 

I 

25 

6.1 

-.108 

2334 

NBS  OIL 

L 

25 

15.5 

-.264 

5407 

5%  Polystyrene  {2200  M.W.)  in 

DOP  w.r.t.  ETOH 

25 

11.5 

-.145 

4253 

20 

12.6 

-.171 

4797 

15 

13.6 

-.200 

5597 

10 

15.0 

-.249 

6473 

10%  Polystyrene 

25 

15.5 

-.205 

5278 

20 

16.3 

-.236 

5902 

15 

17.1 

-.266 

6816 

10 

17.9 

-.333 

7920 

Measured 

R 


3379 

4077 

4745 

5424 

2415 

5905 

REF. 

3871 

4480 

5171 

6297 

5213 

5934 

6647 

8176 
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At  14.5  MHz  X=358 

11=22,370 


Temp. 

(•c) 

. aA 

Attenuation 
(degree)  (nepers) 

Calculated 

Newtonian 

*L 

Measured 

R 

20%  Polystyrene (2200  M.W.)  in  DOP  w.r.t.  ETOH 

REF. 

25 

19.8 

-.331 

8263 

8032 

20 

20.6 

-.396 

9748 

9513 

15 

21.5 

-.484 

11537 

11524 

10 

21.5 

-.646 

13860 

15178 

40%  Polystyrene 

25 

25.7 

-.969 

30002 

22304 

20 

24.5 

-1.126 

39068 

25843 

15 

22.5 

-1.307 

51063 

29934 

10 

20.5 

-1.521 

68808 

34751 
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At  25.0  MHz 


• AA 

Temp.  Attenuation  ' Calculated 

(°C)  (degree)  (nepers)  Newtonian 


DOP-v.r.t.  ETOH 


25 

9.4 

-.254 

4006 

20 

9.5 

-.286 

4511 

15 

9.9 

-.326 

5131 

10 

10.3 

-.367 

5947 

NBS  OIL 

I 

25 

6.1 

-.175 

3065 

NBS  OIL 

L 

25 

11.0 

-.443 

7100 

5%  Polystyrene (2200 

M.W.)  in  DOP 

W.r.t.  ETOH 

25 

10.0 

-.369 

4760 

20 

11.0 

-.402 

5440 

15 

12.0 

-.440 

6435 

10 

13.3 

-.503 

7545 

10%  Polystyrene 

25 

13.2 

-.506 

5975 

20 

13.9 

-.537 

6890 

15 

14.8 

-.601 

8035 

10 

15.9 

-.686 

9445 
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At  25.0  MKz 


* 


Temp. 

CC) 

(degree) 

Attenuation 

(nepers) 

Calculated 

Newtonian 

. 

20%  Polystyrene  (2200  M.W.)  in  DOT 

w.r.t.  ETOH 

25 

17.7 

-.764 

10025 

20 

16.9 

-.829 

11940 

15 

16.5 

-.916 

14235 

10 

16.4 

-1.114 

17243 

40%  Polystyrene 

25 

12.0 

-1.579 

38675 

20 

in 

. 

-1.664 

50440 

15 

-1.7 

-1.963 

66135 

10 

-8.6 

-2.208 

89395 
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At  27.1  MHz 


Temp.  Attenuation 

(°C)  (degree)  (nepers) 


DOP  - w.r. 

t.  ETOH 

25 

8.9 

-.280 

20 

9.3 

-.313 

15 

10.0 

-.349 

10 

11.3 

-.393 

NBS  OIL  I 

- w.r.t  AIR 

25 

5.6 

-.178 

NBS  OIL  L 

- w.r.t.  AIR 

5%  Polystyrene (2200  M.W.)  in 

25 

8.9 

-.395 

20 

10.2 

-.423 

15 

11.0 

-.470 

10 

12.5 

-.549 

10%  Polystyrene 

25 

9.4 

-.516 

20 

10.9 

-.584 

15 

12.9 

-.637 

10 

15.0 

-.704 

calculated 

Newtonian 


4171 

4695 

5342 

6191 

3191 

.t.  ETOH  REF. 
4956 
5664 
6699 
7856 

6221 

7174 

8365 

9834 
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At  27.1  MHz 


Temp. 

(°C) 

a4 

(degree) 

A* 

Attenuation 

(nepers) 

Calculated 

Newtonian 

xL 

20%  Polystyrene (2200  M.W.)  in  DOP  w.r.t.  ETOH  REF 

25 

14.1 

-.826 

10437 

20 

14.9 

-.916 

12432 

15 

16.0 

-1.016 

14820 

10 

17.5 

-1.149 

17955 

40%  Polystyrene 

25 

5.3 

-1.743 

40267 

20 

-0.6 

-1.880 

52516 

15 

-8.5 

-2.019 

68856 

10 

-23.5 

-2.200 

93074 
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The  raw  experimental  measurements  were  found  in  Table  3 and 
Figs. 12,  13,  14  and  15.  Hie  summary  of  the  results  are  shown  in  Table  5. 

B.  Steady  State  Measurements  of  Solvent  and  Solutions 

1.  POP  Plasticizer 

The  steady  state  viscosity  of  the  DOP  (C,H.  - 1. 2 (COO  (CH.)_CH,)^) 

0 4 7 2 7 3 2 

plasticizer  used  in  these  experiments  has  been  measured  over  a range  of 
temperatures  (25° C,  20° C,  15° C and  10° c)  with  a capillary  viscometer;  and 
the  density  was  also  measured  over  the  same  temperature  range  with  volu- 
metric flask  and  (Mettler  H-20)  analytic  balance.  The  results  are  pre- 
sented in  Table  6. 

2.  Polystyrene/DOP  Solutions  (5%, 10%, 20%  and  40%) 

The  viscosity  and  the  density  of  the  solutions  were  measured  in 
the  same  manner  as  above  and  the  results  are  shown  in  Table  6 . 

C.  Sample  Description  and  Preparation 

The  anionically  polymerized  polystyrene  sample  was  obtained 
from  Pressure  Chemical  Company.  The  weight-average  molecular  weight  of 
this  polymer  is  2200.  It  is  in  powder  form. 

Hie  DOP  was  purchased  in  liquid  form  from  Eastman  Chemical 
Company,  and  bears  the  chemical  formula  C^H^-1,  2 [COO  (CH^ ) ^CH^  1 2 and  1:116 
molecular  weight  is  390.56.  The  steady  state  viscosities  at  various  tem- 
peratures are  shown  in  Table  6. 

All  solutions  were  prepared  by  weighing,  and  the  weight  fraction 
of  polystyrene  was  converted  to  g/cc  by  multiplying  by  the  solution  density. 

Solutions  using  polystyrene  and  DOP  were  dissolved  quite  readily 
with  gentle  continuous  swirling  over  an  average  period  of  two  days.  Hie 


■ — ] 

swirling  was  done  at  50° C with  the  aid  of  a small  motor  which  was  oper- 
ated at  a moderate  speed  to  avoid  any  possibility  of  mechanical  degradation. 

The  heating  was  necessary  since  polystyrene  in  powder  form  would 
not  dissolve  at  room  temperature  in  a 48  hour  period. 


Table  6.  STEADY  STATE  DATA 


TEMP. 

C 

DENSITY  *■ 

g/cc 

VISCOSITY 

(poises) 

in' 

{-On 

Ethanol (ETOH) 

25 

0.785 

0.0100 

0.088 

0.165 

20 

0.789 

0.0120 

0.097 

0.172 

15 

0.793 

0.0133 

C.103 

0.182 

10 

0.798 

0.0147 

0.108 

0.191 

DOP 

25 

0.971 

0.306 

0.545 

0.966 

20 

0.976 

0.376 

0.606 

1.074 

15 

0.977 

0.475 

0.682 

1.209 

10 

0.980 

0.619 

0.779 

1.380 

NBS  Oil  ] 

[ 

25 

0.828 

0.100 

0.346 

0.613 

NBS  Oil  L 

25 

0.868 

0.740 

0.801 

1.420 

5%  Polystyrene (2200  M.W.) 

in  DOP 

25 

0.975 

0.407 

0.630 

1.117 

20 

0.978 

0.518 

0.712 

1.260 

15 

0.982 

0.701 

0.830 

1.470 

10 

0.986 

0.931 

0.958 

1.700 

10%  Polystyrene (2200  M.W. 

) in  DOP 

25 

0.976 

0.602 

0.766 

1.360 

20 

0.979 

0.783 

0.875 

1.550 

15 

0.982 

1.042 

1.011 

1.790 

10 

0.986 

1.403 

1.176 

2.080 
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TEMP. 

C 

DENSITY 

g/cc 

'’’'VISCOSITY 

(poises) 

in 

20%  Polystyrene (2200  M.W. 

) in  DOP 

25 

0.985 

1.528 

1.227 

2. no 

20 

0.989 

2.114 

1.446 

2.560 

15 

0.992 

2.951 

1.711 

3.030 

10 

0.996 

4.236 

2.054 

3 ..640 

40%  Polystyrene (2200  M.W. 

) in  DOP 

25 

1.000 

19.880 

4.459 

7.900 

20 

1.010 

33.200 

5.791 

10.260 

15 

1.010 

56.700 

7.567 

13.410 

10 

1.010 

102.96 

10.197 

18.070 

NBS  Oil  D 

25 

0.760 

NBS  Oil  N 

25 

0.883 

9.162 

2.844 

Aroclor  1242 

25 

1.381 

0.6044 

0.913 

*■  CANNON  VISOVETE*  #3ooD/SJ 
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PART  V 


DISCUSSION  OF  RESULTS 

A.  Analysis  of  the  Data 

Experimental  measurements  are  summarized  in  Tables  7 and  8.  In 
Table  7 the  ratio  T]  /l]  (0)  of  the  solvent  viscosity  7)  in  the  absence  of 

S S 7 

the  polymer  to  the  solution  viscosity  in  the  low  frequency  limit  11(0), 

is  compared  to  the  ratio  T)/ (10)/7]  (0),  the  solution  viscosity  at  10  f«z 

to  the  solution  viscosity  in  the  low  frequency  limit.  This  comparison 

is  made  at  different  concentrations  and  temperatures.  In  this  form,  the 

value  of  7] ' (10)/T|  (0)  = 1 in  the  limit  of  no  relaxation  at  10  mz,  and  can 

approach  Tl^/Tl  (0 ) in  the  high  frequency  limit.  Similarly,  in  Table  8,  the 

same  results  are  presented  normalized  by  71  , that  is,  7]  ( 0 ) /T)  and  7] 7 (10 ) /T1  . 

s s s 

In  this  form,  the  value  of  7T| ' (10 ) /T]  =*  — in  the  limit  of  no  relaxation, 

s s 

and  can  approach  1 in  the  high  frequency  limit.  In  fact,  if  we  call  the 
value  of  the  high  frequency  limiting  solution  viscosity  T^,  then  the 
limiting  values  become  T)  /T)  (0)  for  Table  7 and  — for  Table  8. 

The  most  significant  conclusion  of  this  present  data  is  that 
T)7  (10)/7)  (0)  shows  little  evidence  of  relaxation  of  viscosity  even  at  con- 
centrations of  20%  and  temperature  of  10°C.  This  leaves  us  with  the  pri- 
mary conclusion  that  any  relaxation  mechanism  responsible  for  the  solu- 
tion viscosity  must  occur  at  a frequency  greater  than  10  MHz,  or  equiva- 

_0  *7 

lently,  the  relaxation  time  t must  be  shorter  than  1.6  X 10  sec  (1/2tt10  Hz). 
Table  7 indicates  that  only  when  concentration  reaches  40%  or  beyond  does 
relaxation  become  apparent,  that  is,  the  relaxation  time  of  the  solution 
approaches  the  time  scale  of  our  measurement. 
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NORMALIZED  VISCOSITIES  WITH  RESPECT  TO  SOLVENT  VISCOSITY 


In  order  to  account  for  the  enhanced  viscosity  of  the  solution 


TJ (0)  over  that  of  the  solvent  T|g  at  low  frequencies,  we  examine  the  predic- 
tion of  several  known  mechanisms  of  viscosity  enhancement.  We  shall 
briefly  discuss  each  of  these  possible  mechanisms,  Then,  in  turn  calculate 

the  contribution  T|(0)  ” T)  to  the  viscosity  enhancement.  We  summarize  the 

s 

results  in  Table  9.  Each  calculation  is  based  on  the  experimental  observa- 
tion (Table  7 or  8)  that  the  longest  relaxation  time  predicted  by  each 

_8 

model  must  be  shorter  than  1.6  X 10  sec  as  long  as  the  concentration 
does  not  exceed  20%. 

Briefly,  predictions  of  the  theories  for  known  mechanisms  are 


as  follows: 


1.  The  rigid  rod  model  of  Kirkwood  and  Auer 


13 


Y\-Y]  - ± C RT  rr  • r - /L 


_ ® 7s  L3 


*T 

where  c is  concentration  in  gm/cc  of  the  solution,  M is  the 

molecular  weight,  T1  is  the  solvent  viscosity,  L is  the  length 

of  the  molecule  and  b is  the  distance  between  the  submolecules. 

14  15 

2.  The  prolate  ellipsoid  model  of  Cerf  and  Scheraga 


Xl — Vi  - £1  c RT.  • '7L 

I *5  m T > L TX77 


1 


where  p is  the  axial  ratio. 

17 

3.  The  rigid  dumbbell  model 
y»_u  C RT  _ 

*1-1*=  -pr*  ; 


where  § is  a friction  coefficient  for  a rigid  dumbbell  of 
length  L. 
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4.  Hie  head-spring  model  with  no  hydrodynamic  interaction  among 


heads  (Rouse  model ) : 

l-is=  i- <4s 


Cn'f'fiJ 


5. 


where  p =*  1,  2,  3,  . . . . , and  where  N is  the  number  of  sub- 
molecules, f is  the  monomeric  friction  coefficient  and  c 
' o 

is  the  root-mean- square  end-to-end  length  of  the  submolecule. 
The  bead-spring  model  with  dominant  hydrodynamic  interaction 
(Zimm  model) : 


% 


o.t 06 

X*.  -fe.  l 


where  are  numerical  coefficients  (eigenvalues),  whose 
first  few  values  are  4.04,  12.79,  24.2  and  37.9,  T)  is  the 
solvent  viscosity,  p is  the  integer  1, 2, 3,  ...  . At  finite 
concentrations  the  value  of  t in  the  various  cases  becomes : 


(1>  hi  (VlOM  /c  RT 

where  m = 5/4  for  rod  (Kirkwood  and  Auer),  m = 25/21  for 
very  long  ellipsoids  (Scheraga)  and  m = 1 for  rigid  dumb- 
bells (Ullman) 


Hi)  /e  RT 


where  s = 1.645  for  Rouse's  model  and  s * 2.369  for  Zimm' s 
model. 
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From  the  calculated  values  of  T)  (0 ) — T)  in  Table  9,  it  is  clear 

s 

that  none  of  these  models  are  able  to  account  for  the  experimentally 
measured  viscosity  enhancement.  In  the  absence  of  any  other  mechanism  to 
account  for  polymer  contribution  to  the  total  viscosity,  we  conclude  that 
the  primary  effect  of  adding  the  polymer  to  the  solvent  is  to  affect  the 
relaxation  time  of  the  solvent.  Namely,  the  relaxation  time  of  the  sol- 
vent is  moved  to  longer  times,  and  hence  at  low  frequencies  the  solvent 
viscosity  itself  is  increased. 

We  have  no  specific  model  to  account  for  this  increase.  One 
possible  mechanism  could  be  that  the  role  of  the  low  molecular  weight 
polymer  is  to  reduce  the  available  free  volume  and  thereby  slow  down  the 
solvent  re-orientation  process. 

In  this  discussion,  viscosity  is  proportional  to  the  relaxation 
time.  The  relaxation  time  is  related  to  the  free  volume  available  for 
rearrangement  of  solvent  molecules,  we  postulate  here  that  the  primary 
effect  of  polymer  molecule  is  to  decrease  the  free  volume  available  for 
solvent  molecule  rearrangement.  Thereby  relaxation  time  is  increased, 
so  as  the  solution  viscosity  and  so  as  the  solvent  viscosity. 

Another  possibility  is  that  the  polymer  molecule  plays  some 


role  in  correlating  motions  of  solvent  molecules  and  thereby  extending 
the  relaxation  time  to  longer  times  and  possibly  broadening  the  solvent 
relaxation  spectrum. 


Table  9 (cont'd 


B.  Assumptions  and  Accurac 


Hie  development  of  the  equations  is  based  on  the  following 
assumptions. 

1.  The  excitation  of  transducer  is  a steady-state  sinusoidal 

motion,  whereas  a pulse-echo  technique  is  used  in  this  experiment.  The 

total  response  approaches  the  steady-state  sinusoidal  vibration  when  the 

22  11 

pulse  length  is  comparable  with  10  cycles  and  according  to  Hunston's 
results,  a pulse  as  short  as  14  cycles  at  2 MHz  has  been  used  success- 
fully. In  the  experiment  present  here,  the  20  p,sec  pulse  duration  used 
at  a frequency  of  10  MHz  corresponds  to  200  cycles. 

2.  A fused  quartz  delay  material  is  assumed  to  be  homogeneous 
and  isotropic  to  an  acoustical  wave  in  the  frequency  range  of  our  interest. 
When  the  wavelength  (.03  cm)  of  the  propagating  acoustic  wave  is  comparable 
with  the  grain  size  of  fused  quartz,  the  above  assumptions  of  the  homo- 
geneity and  isotropy  become  invalid.  However,  the  wave  propagating  in  the 
fused  quartz  at  10  MHz  is  at  least  two  orders  larger  than  the  grain  size. 

3.  The  mode  for  the  delay  line  is  infinite  in  width  and  length. 
This  assumption  is  based  on  the  fact  that  the  length  (7.62  cm)  and  the 
width  (5.08  cm)  are  large  with  respect  to  wavelength  (.03  cm)  in  the  delay 
line.  Hie  transducer  is  also  designed  to  excite  only  the  principal  shear 
mode  SS (0) . However,  with  a thickness  of  the  delay  medium  = 0.152  cm, 

other  plate  modes  could  propagate,  as  seen  in  Meeker's  analysis: 

U>  b 

s 3 .2 


and  therefore,  the  cut-off  frequency  for  SS(1)  becomes 


f,  = 22  * io‘  Hr 
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where  2b  and  Vg  are  the  thickness  of  the  delay  medium  and  the  velocity, 
respectively.  Therefore,  it  is  possible  to  have  SS(1)  mode  propagation 
at  any  frequency  higher  tnan  2.2  MHz.  However,  frequencies  at  which 
measurements  are  taken  have  been  selected  so  that  higher  order  mode  contri- 
butions are  not  apparent. 

4.  It  is  also  necessary  to  consider  reflection  from  the  end  of 
the  liquid  if  the  thickness  of  the  liquid  (distance  from  fused  quartz  to 

-5 

plastic  case)  is  not  large  compared  to  the  damping  length  (1  X 10  cm). 

In  our  case  the  thickness  is  0.0375  cm,  which  is  approximated  as  infinite 
with  respect  to  damping  length. 

5.  Smooth  surface:  The  fused  quartz  surface  is  much  easier  to 
polish  that  others,  i.e.,  steel,  aluminum.  The  polishing  is  of  critical 
importance  since  for  the  viscoelastic  liquids  the  wavelength  or  damping 
length  begins  to  approach  the  magnitude  of  the  surface  irregularities,  and 
for  low  viscosity  liquids  such  as  ethanol  (X  = 2.5  x 10  ^ cm).  Surface 
smoothness  of  the  delay  medium  was  not  well  esablished,  and  the  small 
streaks  of  random  lines  are  observed  under  a 20  x magnification  microscope. 

6.  Possible  residual  coating  of  the  surface  by  the  polymer  solu- 
tion can  occur  during  the  drying  of  the  cell.  In  order  to  minimize  this 
problem,  the  repeated  cleaning  of  the  cell  with  low  alcohols  were  carried 
out. 

7.  The  temperature  of  the  delay  line  must  be  accurately  controlled 
and  measured,  even  though  the  temperature  coefficient  of  acoustic  velocity 
is  smaller  in  fused  quartz  than  in  many  other  materials.  We  can  see  from 
the  example  of  Fig. 8,  that  at  19.5  Wiz,  the  phase  shifts  by  12°  per  °C 
temperature  change.  Therefore,  it  is  necessary  to  control  temperature  to 

: 
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within  . 01° C in  order  to  measure  the  phase  shift  to  within  a few  tenths  of 
a degree.  It  is  also  important  that  equilibrium  is  reached  before  any 
measurements  are  made. 

8.  Hie  phase  and  amplitude  response  of  p. Z.T.  transducer  will 
be  influenced  by  the  total  mechanical  impedance  of  the  delay  line;  this 
impedance  will  change  slightly  by  the  presence  of  liquid  on  the  delay 
medium,  called  a "loading  effect".  For  the  experiment  presented  here  the 
loading  effect  has  not  been  included  in  the  analysis. 


PART  VI 


CONCLUSIONS  AND  SIGNIFICANCE 

1.  An  ultrasonic  delay  line  device  capable  of  measuring  visco- 
elastic properties  of  liquids  has  been  constructed  and  tested.  The  fre- 
quency range  is  from  3 to  30  Miz.  The  viscosity  range  is  down  to  0.5  centi- 
poise,  and  has  been  used  to  beyond  100  poise. 

2.  The  method  of  "walking  frequency"  has  been  developed  to 
transform  a high  frequency  signal  into  1. 25  Hz  in  order  to  overcome  the 
limitation  of  digitizing  a high  frequency  periodic  signal  in  ADC,  and 
read  into  the  computer  at  a rate  of  one  point  every  40  y,sec. 

3.  A series  of  measurements  were  performed  on  low  molecular 
weight  polystyrene  solutions  over  the  concentration  range  5%  to  40%. 

4.  A basic  problem  in  the  molecular  theory  of  viscoelastic 
fluids  was  examined  experimentally  with  this  apparatus.  The  problem, 
called  the  "71^  problem",  asks:  "Why  does  the  viscosity  contribution  of 
the  polymer  71(0)  - T|s  not  drop  to  zero  in  the  high  frequency  limit  as  seems 
indicated  by  macromolecular  theories?" 

While  some  investigators  have  attempted  to  include  corrections 
to  the  macromolecular  theories,  we  conclude  from  these  experimental  results 
that  the  apparent  discrepancy  is  due  to  the  neglect  of  a polymer  solvent 
interaction  mechanism,  rather  than  due  to  any  basic  defect  in  the  theories 
for  macromolecular  backbone  relaxation. 
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PIG. 16.2 PHASE  vs  TEMPERATURE 
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FIC.ITi-  RELATIVE  AMPLITUDE  vs  TEMPERATURE 
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FIG. |7. 5 RELATIV3  AMPLITU" 
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FIG.  | g I ATTENUATION  VS  SOUARE  ROOT 

OF  FREQUENCY 
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FIG.  18.2  ATTENUATION  VS  SQUARE  ROOT 
OF  FREQUENCY 
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FIG.  13.2  PHASE  DIFFERENCE  VS  SQUARE 
ROOT  OF  FREQUENCY 
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FIG.  10.3  PHASE  DIFFERENCE  VS  SQUARE 
ROOT  OF  FREQUENCY 
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FIG. 19.4  PHASE  DIFFERENCE  VS  SQUARE 
ROOT  OF  FREQUENCY 
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FI6 . 20.1  RELATIVE  ATTENUATION  VS  SQUARE 
ROOT  OF  FREQUENCY 


APPENDIX  A 

DATA  ACQUISITION..  AND  CROSS-CORREIATION  PROGRAM 

\ ‘ • * \ 2A  \ " 

READ  DATA 


♦INITIALIZE  REGISTER 


ENTER 

IAI1 

INHIBIT  IU 

13000 

2503 

0000 

STX 

PLACE, 1 

01 

0614 

0223 

LXK 

0,7 

02 

0770 

0000 

LXK 

0,6 

03 

0760 

0000 

LXK 

0.5 

04 

0750 

0000 

LXK 

0.4 

05 

0740 

0000 

♦INITIALIZE  HARDWARE  CLOCK 

RESET 

LDK 

3000  RESET  CON. 

06 

4000 

3000 

OUT 

1024 

07 

2504 

1024 

LDA 

CONTROL  SCOFE 

10 

0004 

0215 

OUT 

1001 

11 

2504 

1001 

♦ZERO  DATA  REGION  21,000  to  37,777 

LDZ 

12 

0500 

0000 

STA 

21000,4 

13 

3242 

1000 

INX 

1.4 

14 

2640 

0001 

TXH 

17000,4 

15 

244? 

1000 

BTR 

*~3 

16 

3007 

7775 

♦INITIALIZE  IU 

LDA 

MASK 

17 

0004 

0210 

LNR 

20 

2500 

0302 

LDA 

IU  CON. (BRU  PATCH  3) 

21 

0004 

0211 

STA 

265 

22 

3200 

0265 

START 

CLOCK 

LDA 

CONTROL  CLOCK 

23 

0004 

0203 

OUT 

1024  A/D  ch.  b 

24 

2504 

10?4 

BRU 

PATCH  1 

25 

1404 

0212 

READ 

DATA 

LDA 

20,000 

26 

0002 

0000 

SRA 

10i  q 

27 

0501 

4052 

INX 

1,6 

30 

2660 

0001 

BRU 

POINTER, 6 

31 

1464 

0000 

•IUjBRU 

to  READ 

DATA 

POINTER 

BRU 

PATCH  2 

32 

1404 

0213 

M 

READ  bO 

33 

1404 

0042 

W 

- bl 

34 

1404 

0045 

N 

" b2 

35 

.1404 

0050 

W 

” b? 

36 

1404 

0053 

•f 

- b4 

37 

1404 

0056 

M 

" b5 

40 

1404 

0061 

H 

" b6 

41 

1404 

0064 

M 

" b7 

42 

1404 

0067 

N 

- b8 

43 

1404 

0072 

M 

" b9 

44 

1404 

0075 

M 

SET  A/D  ch.  a 

45 

1404 

0100 

M 

COMPLETE  PULSE  TRAIN 

46 

1404 

0130 

m 

M M H M 

47 

1404 

0127 

M 

READ  al4  ' 

• ‘ 50 

I4a4 

0101 

n 

RESET  A/D  ch.  b 

51 

1404 

0104 

230 


i 


BRU 

COMPLETE  P'JLSE 

TRAIN 

13052 

14C4 

0124 

99 

INX  "7" 

53 

1404 

0106 

H 

COMPLETE  PULSE 

TRAIN 

54 

1404 

0122 

M 

TEST  AVERAGE  COMPLETE 

55 

1404 

0111 

H 

COMPLETE  PULSE 

TRAIN 

56 

1404 

0120 

W 

H 99 

H 

57 

1404 

Oil? 

M 

« « 

•« 

60 

1404 

0116 

M 

H M 

M 

61 

1404 

0115 

M 

W H 

N 

62 

1404 

0114 

H 

M N 

99 

63 

1404 

0113 

M 

M 99 

M 

64 

1404 

0112 

•1 

99  99 

99 

65 

1404 

0111 

M 

99  99 

N 

66 

1404 

0110 

99 

99  99 

M 

67 

1404 

0107 

H 

M •• 

99 

70 

1404 

0106 

M 

40 th  POINTER  ENTER 

71 

1404 

0104 

READ 

NOTHING 

NOP 

72 

2620 

0000 

NOP 

73 

2620 

0000 

NOP 

74 

2620 

0000 

READ 

DATA 

bo 

ADD 

2?200,7 

75 

1172 

2200 

STA 

22200,7 

76 

3272 

2200 

NOP 

77 

2620 

0000 

BRU 

* -1 

100 

140? 

7777 

READ 

DATA 

bl 

ADD 

23400,7 

101 

1172 

3400 

STA 

23400,7 

102 

3272 

3400 

NOP 

103 

2620 

0000 

BRU 

* -1 

104 

1407 

7777 

READ 

DATA 

b2 

ADD 

24600,7 

105 

1172 

4600 

STA 

24600.7 

106 

3272 

4600 

NOP 

107 

2620 

0000 

BRU 

* -1 

110 

1407 

7777 

READ 

DATA 

b3 

ADD 

26000,7 

111 

1172 

6000 

STA 

26000,7 

112 

3272 

6000 

NO? 

113 

2620 

0000 

BRU 

* -1 

114 

140? 

7777 

READ 

DATA 

b 4 

ADD 

27200,7 

115 

1172 

7200 

STA 

27200,7 

116 

3272 

7200 

NOP 

117 

2620 

0000 

BRU 

* -1 

120 

1407 

7777 

READ 

DATA 

b5 

ADD 

30400,7 

121 

1173 

04C0 

STA 

30400,7 

122 

3273 

0400 

NOP 

123 

2620 

0000 

BRU 

* -1 

124 

1407 

7777 

231 


READ 


DATA  b6 

ADD 

31600,7 

125 

1173 

1600 

STA 

31600,7 

* 126 

3273 

1600 

NOP 

127 

2620 

0000 

BRU 

* -1 

. 130 

140? 

7777 

READ 

DATA  b? 

ADD 

33000,7 

131 

1173 

3000 

STA 

33000,7 

132 

3273 

3000 

NOP 

133 

2620 

0000 

BRU 

* -1 

134 

140? 

7777 

READ 

DATA  b8 

ADD 

34200,7 

135 

1173 

4200 

STA 

34200,7 

136 

3273 

4200 

NOP 

137 

2620 

0000 

BRU 

* -1 

140 

140? 

7777 

READ 

DATA  b9 

ADD 

35400,7 

141 

1173 

5400 

STA 

35400,7 

142 

3273 

5400 

NOP 

143 

2620 

0000 

BRU 

* -1 

144 

140? 

7777 

SET 

A/D  ch.  a 

LDA 

CONTROL  CLOCK 

145 

0004 

0061 

OUT 

A/D  ch.  a 

146 

2504 

1004 

NOP 

147 

2620 

0000 

BRU 

* -1 

15° 

1407 

7777 

READ 

DATA  a 14 

ADD 

36600,7 

151 

1173 

6600 

STA 

36600,7 

162 

3273 

6600 

NOP 

153 

2620 

0000 

BRU 

* -1 

154 

140? 

7777 

RESET 

A/D  ch,  b 

LDA 

CONTROL  CLOCK 

155 

0004 

0051 

OUT 

A/D  ch.  b 

156 

2504 

1024 

NOP 

157 

2620 

0000 

BRU 

* -1 

160 

140? 

7777 

INX  "7" 

INX 

1.7 

161 

2670 

0001 

TXH 

1161,7 

162 

2473 

661? 

BTR 

* +5 

163 

3004 

0005 

LXK 

0.7 

164 

0770 

0000 

INX 

1.5 

165 

2650 

0001 

TEST 

AVERAGE 

COMFLETE 

TXH 

456.5 

166 

2453 

7322 

BTS 

FINI 

167 

3404 

0025 

NOP 

170 

2620 

0000 

BRU 

* -1 

171 

140? 

7777 

NOP 

172 

2620 

0000 

TXH 

40,6 

173 

2463 

7740 

BTR 

* +2  WAIT  LOOP 

174 

3004 

0002 

4oth 

POINTER 

ENTER 

LXK 

0,6 

175 

0760 

0000 

232 


COMPLETE 

PULSE 


TRAIN  NOP 

13176 

2620 

0000 

BRU 

* -1 

177 

1407 

7777 

PINI  I All 

INHIBIT  IU 

214 

2503 

0000 

LDA 

NOP 

215 

0004 

0014 

STA 

265 

216 

3200 

0265 

LMO 

1*8 

217 

0506 

0000 

LMR 

220 

2500 

0302 

LDK 

3000 

221 

4000 

3000 

OUT 

222 

2504 

1024 

LDP 

223 

1504 

0001 

PATCH  1 LXK 

0.3 

237 

0730 

0000 

LDZ 

INITIALIZE  PULSE 

240 

0500 

0000 

OUT 

1001 

241 

2504 

1001 

PAI 

242 

2502 

0000 

NOP 

243 

2620 

cooo 

BRU 

* -1 

244 

1407 

7777 

PATCH  2 ADD 

21000,7 

245 

1172 

1000 

STA 

21000,7 

246 

3272 

1000 

NOP 

247 

2620 

0000 

BRU 

* -1 

250 

1407 

7777 

PATCH  3 INX 

1.3 

251 

2630 

0001 

TXH 

1.3 

252 

2433 

7767 

BTR 

* +3 

253 

3004 

0003 

LDA 

"BRU  to  READ  DATA" 

254 

0001 

3230 

STA 

265 

255 

3200 

0265 

NOP 

256 

2620 

0000 

BRU 

* -1 

257 

1407 

7777 

STABLE 

PLACE  BSS 

224 

7777 

7777 

CONTROL  SCOPE 

DATA 

225 

4000 

0000 

CONTROL  CLOCK 

DATA 

226 

0014 

2400 

MASK 

DATA 

227 

7775 

7777 

IU  CONSTANT 

DATA  BRU  to  READ  DATA 

230 

1401 

3026 

CON  "NOP" 

231 

2620 

0000 

DATA  "BRU  to  PATCH  3" 

232 

1401 

3251 

233 


D/A  DISPLAY 


*D/A  DISPLAY  MAIN 


STX1 

PLACE 

17200 

0614 

0023 

LXK 

0.7 

01 

0?70 

0000 

SAVE  X7 

STX 

TABLE  COUNTER(TEMPl) , 

7 02 

0674 

0020 

LDA 

DISPLAY  SEQ, 7 

03 

0071 

7224 

STA 

EXEC.  #1 

04 

3204 

0014 

LDA 

DISPLAY  SEQ  +1,7 

05 

0071 

7225 

STA 

EXEC.  j}2 

06 

3204 

0013 

SPB 

REPACK  TABLE 

07 

3304 

0033 

NOP 

10 

2620 

0000 

NOP 

11 

2620 

0000 

SPB 

D/A  DISPLAY 

12 

3304 

0054 

NEXT 

DISPLAY 

LDX 

TEMP  1,7 

13 

1674 

0007 

I NX 

2.7 

14 

2670 

0002 

TXH 

X7  * 15 

15 

2473 

7763 

BTR 

* -14(1) 

16 

3007 

7764 

LPR 

17 

3504 

0004 

BSS 

EXEC.  #1 

17220 

7777 

7777 

BSS 

EXEC.  it?. 

221 

7777 

7777 

BSS 

TEMP  1,  TABLE  COUNTER 

22 

7777 

7777 

BSS 

PLACE 

23 

7777 

7777 

"DISPLAY  SEQUENCE 

DISPLAY 

SSQ. 

TABLE 

1 

DATA  al4 

24 

0073 

6600 

25 

2473 

6617 

TABLE 

2 

DATA  bO 

26 

0072 

2200 

27 

2473 

6617 

TABLE 

3 

DATA  bl 

30 

0072 

4600 

31 

2473 

6617 

TABLE 

4 

DATA  b2 

32 

0072 

4600 

33 

2473 

6617 

TABLE 

5 

DATA  b3 

34 

0072 

6000 

35 

2473 

6617 

TABLE 

6 

DATA  b4 

36 

0072 

7200 

37 

2473 

6617 

TABLE 

7 

DATA  Sin.  F 

40 

0073 

5400 

41 

2473 

6617 

EXEC. 

1 

LDA  TAELE  7 

EXEC. 

2 

TXH  N,7 

234 


•REPACK  TABLE  AND  APD  AT 


STX1  PLACE 

17242 

0614 

0023 

LXK  0,7  Y COUNTER 

243 

0770 

0000 

LXK  0,7  AT  COUNTER 

244 

0760 

0000 

(1) 

XEC  17220  to  (LDA  TABLE, 

7)  245 

0407 

7753 

SLA  4 BITS 

246 

4500 

2044 

ANA  MASK 

247 

2004 

0030 

BRU  PATCH (17300) 

250 

1401 

7300 

STA  10000,7  DISPLAY  TABLE, 7 251 

3271 

0000 

INX  AT, 6 

252 

2660 

0002 

INX  1,7 

253 

2670 

0001 

EXC  17221  to  (TXH  N,7) 

254 

0407 

7745 

•STORE 

BTR  REPEAT  (1) 

Y=D  IN  THE  REST  OF  TABLE 

255 

3007 

7770 

(2) 

LDA  X6 

256 

0000 

0006 

STA  DISPLAY  TABLE, 7 

257 

3271 

0000 

INX  AT, 6 

260 

2660 

00  03 

INX  1,7 

261 

2670 

0001 

TXH  1161,7 

262 

2473 

6617 

BTR  * (2) 

263 

300? 

7773 

LPR 

264 

3504 

0001 

BSS  PLACE 

•D/A  DISPLAY  PROTION 

265 

7777 

7777 

LXK  0,7 

266 

0770 

0000 

LDA  DISPLAY, X7 

267 

0070 

0000 

OUT  D/A 

270 

2504 

1003 

INX  1,7 

271 

2670 

0001 

TXH  ll6l,7(#  of  display 

pts ) 272 

2473 

6617 

BTR  * -4 

273 

3007 

7774 

•CHECK 
ON  TO 

DEMAND  IF  NOT  SET  THEN  REPEAT 
NEXT  FCN 

DISPLAY, 

IF  NOT  GO 

JND 

2?4 

2504 

0000 

LPR 

275 

3500 

0001 

BRU  * -10 

276 

1407 

7770 

CONI  MASK 

277 

7777 

0000 

PATCH 

JND 

300 

2506 

0000 

LDA 

301 

0004 

0003 

ADD 

302 

1100 

0006 

BRU 

303 

1401 

7251 

CON=3777  0000 

304 

7777 

7777 

235 

L 


♦INTO  LOCATIONS  35400  to  36,561 


STX  PLACE,  1 

13350 

06l4  0025 

FMS 

51 

7401  4000 

LXK  0,7 

52 

0770  0000 

LDZ 

53 

0500  0000 

ST  A,  Q 

54 

3200  0010 

DST  TEMP 

55 

6304  0016 

SPB  to  SIN  F 

56 

3300  2340 

FIX 

57 

7400  0010 

NOP 

60 

2620  0000 

STA  35400,7 

61 

3273  5400 

DLD  CON  1 (floating  pt . ) 

62 

4104  0007 

FAD  TEMP 

63 

7004  0010 

DST  TEMP 

64 

6304  000? 

I NX  1,7 

65 

2670  0001 

TXH  1161,7 

66 

2473  6617 

BTR 

67 

3 00?  7767 

LPH 

70 

3504  0005 

DATA  CON  1 

71 

1752  4467 

DATA  CON  1 

72 

3537  1300 

TEMP 

73 

7777  7777 

TEMP 

74 

7777  7777 

BSS  PLACE 

75 

7777  7777 

MINI  EXECUTIVE  PROG. (SELECT  DATA) 


STX1 

PLACE 

13300 

06l4 

0012 

LXK 

0,4 

01 

0740 

0000 

STX 

TEMP, 4 

02 

0644 

0011 

SPB 

CROSS -CORR  AND 

PRINT 

03 

3301 

3500 

LDX 

TEMP.4 

04 

1644 

0007 

INX 

4.4 

05 

2640 

0004 

TXH 

24,4(also  indicates 

the  # of  lines 

to  be 

printed) 

06 

2443 

7754 

BTR 

* -5 

07 

3007 

7773 

LPR 

10 

3504 

0002 

NOP 

11 

2620 

0000 

BSS 

PLACE 

12 

7777 

7777 

BSS 

TEMP 

13 

7777 

777  7 

■CORRELATION  AND  PRINT 

STX1 

PLACE 

13500 

06 14 

0073 

LXK 

0,7 

01 

0770 

0000 

LXK 

0,6 

02 

0760 

0000 

LXK 

234,5 

03 

0750 

0234 

TA3LE 

LDZ 

04 

0500 

0000 

ST  A 

13563.6 

05 

3264 

0056 

INX 

1.6 

06 

2660 

0001 

TXH 

10,6 

07 

2463 

7770 

BTR 

* -3 

10 

3007 

7775 

;UM3 

LDZ 

11 

0500 

0000 

X2C , 4 

(LDO  A. 7) 

12 

0441 

3400 

XEC , 4 

(MPY  C,7) 

13 

0441 

3402 

DAD  ,4 

(TEMP  1) 

14 

5104 

0047 

DST,  4 

(TEMP  1) 

15 

6304 

0046 

LDZ 

16 

0500 

0000 

XEC,  4 

(LDQ  B,7) 

17 

0441 

3401 

XEC,  4 

(MPY  C ,7) 

20 

0441 

3402 

DAD,  4 

(TEMP  2) 

21 

5104 

0044 

DST,  4 

(TEMP  2) 

22 

6304 

0043 

LDZ 

23 

0500 

0000 

XEC,  4 

(LDQ  A. 7) 

24 

0441 

3400 

XEC,  4 

(MPY  C^/2,5) 

25 

0441 

3403 

DAD,  4 

(TEMP  3) 

26 

5104 

0041 

DST,  4 

(TEMP  3) 

27 

6304 

0040 

LDZ 

30 

0500 

0000 

XEC,  4 

(LDQ  B,7) 

31 

0441 

3401 

XEC, 4 

(MPY  C*/2,5) 

32 

0441 

3403 

DAD,  4 

(TEMP  4) 

33 

5104 

0036 

DST,  4 

(TEMP  4) 

34 

6304 

0035 

237 


INX 

1.7 

13535 

2670 

0001 

INX 

1.5 

36 

2650 

0001 

TXH 

1161.5 

37 

2453 

6617 

BTS 

to 

RESET  5 

40 

3404 

0021 

TXH 

1161.7 

41 

2473 

6617 

BTR 

to 

REPEAT (LOC.  13511 

CALC.  SUM) 

42 

r- 

0 

0 

r\ 

7747 

•DONE, 

SUMS  AND 

PRINT 

BRU 

to 

RECYCLE 

43 

1401 

3320 

BRU 

* 

44 

1404 

0000 

LXK 

0.3 

45 

0730 

0000 

DLD 

TEMP  1,3 

46 

4134 

0015 

DRA 

47 

4500 

4404 

SPB 

to 

(FIX  PT.  BINARY  to 

FIX  PT.  DECIMAL) 

50 

3300 

1537 

LDK 

K=l6l0 

51 

4000 

1610 

SPB 

to 

PRINT 

52 

3300 

1640 

INX 

2.3 

53 

2630 

0002 

TXH 

10 , 3( Indicates  how  many 

# of  groups  In  a row) 

54 

2433 

7770 

BTR 

* -7 (LOC.  13546) 

55 

3007 

7771 

LDA 

CON 

56 

0004 

0016 

SPB 

to 

PRINT 

57 

3300 

1640 

LPR 

60 

3504 

0013 

•RESET 

•i  ^ it 

LXK 

0.5 

61 

0750 

0000 

BRU 

* -21 (LOC . 13541) 

62 

1407 

7757 

•DATA 

BSS 

TEMP  1 AlCi 

63 

7777 

7777 

64 

7777 

7777 

3SS 

TEMP  2 B1C1 

65 

7777 

7777 

TENT  3 AlCl^/2 

66 

7777 

7777 

BSS 

6 7 

7777 

7777 

TEMP  4 B 1C  1^2 

70 

7777 

7777 

BSS 

71 

7777 

7777 

72 

7777 

7777 

BSS 

PLACE 

73 

7777 

7777 

DATA 

CON=400  13576 

74 

4001 

3575 

DATA 

LFLFCR=0240  5015 

75 

024o 

6415 

DATA 

CRLF  0^0320  5000 

76 

0320 

5000 

238 


OCD>  otu>  ow>  oco>  ow> 


DATA  SELECTION  TABLE 


A LDQ 

b0,7 

13400 

4272 

2200 

B LDQ 

bl,7 

01 

4272 

3400 

C MPY 

b9.7 

02 

5573 

5400 

MPY 

b9.5 

03 

5553 

5400 

A LDQ 

b2,7 

04 

4272 

4600 

B LDQ 

b3,7  b9=C 

05 

4272 

6000 

C MPY 

b9.7 

06 

5573 

5400 

MPY 

b9.5 

07 

5553 

5400 

A LDQ 

b4, 7 

10 

4272 

7200 

B LDQ 

b5.7 

11 

4273 

0400 

C MPY 

b9.7 

12 

5573 

5400 

MPY 

b9,5 

13 

5553 

5400 

A LDQ 

b6,7 

14 

4273 

1600 

B LDQ 

b7 ,7 

15 

42?3 

3000 

C MPY 

b9,7 

16 

5573 

5400 

MPY 

b9,5 

17 

5553 

5400 

A LDQ 

b8,7 

20 

4273 

4200 

B LDQ 

al5.7 

21 

4273 

6600 

C MPY 

b9.7 

22 

5573 

5400 

MPY 

b9.5 
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APPENDIX  B 


Sample  Printout  Format  of  Cross-Correlated  Values 

After  the  cross-correlation  has  been  carried  out  in  CPU, 

the  teletype  prints  out  the  terms  in  the  following  format: 


012964 

012964 

067403 

067402 

012964 

012965 

067402 

067403 

012965 

012964 

067402 

067402 

012965 

012964 

067402 

067401 

012965 

-019439 

067402 

-004932 

and  these  terms  are  identified  respectively  as  follows: 
b#cces<^  b«ccos^>o  b*csin<^  b0csin<j> 

II  II  II  II 

II  II  II  II 

II  II  II  II 

" a|#ccos^>4  " a^csint^ 

From  the  above  terms  we  obtain^- ^^and  bg/a^s  which  are  the 
phase  difference  between  the  zeroth  echo  and  the  input  signal, 
and  the  amplitude  ratio  of  the  zeroth  echo  to  input  signal, 
respectively. 


APPENDIX  C 


THEORY  OF  THE  FABRY- PEROT  INTERFEROMETER 
INTRODUCTION 

The  Fabry-Perot  interferometer  is  a multiple  beam  interferometer  char- 
acterized by  interference  between  many  rays.  The  basic  advantage  over  two 
beam  interferometers  such  as  the  Michelson  type  is  greatly  increased  resolu- 
tion, roughly  dependent  on  the  number  of  interfering  beams. 

The  construction  of  a Fabry-Perot  consists  of  two  plane  partially  re- 
flecting surfaces  oriented  parallel  to  each  other.  These  surfaces  may  be 
the  opposite  faces  of  a quartz  plate  or  other  similar  low  thermal  expansion 
material.  Such  a fixed  gap  Fabry-Perot  is  known  as  an  e talon.  The  other 
major  design  incorporates  two  opposed  partially  reflecting  mirrors  a distance  d 
apart,  where  d may  be  varied.  In  both  designs  the  gap  between  the  surfaces  is 
constant  over  the  whole  area  to  1/50  to  1/200  of  the  wavelength  of  green 
light,  depending  on  the  specifications  of  the  particular  instrument.  Our 
instrument,  a Burleigh  Instrument  RC  10  is  of  the  opposed  mirrow  type  with 
the  gap  distance  accurate  to  1/200  wavelength. 


THEORY 


A treatment  of  the  basic  Fabry-Perot  can  be  treated  by  considering  the 
infinite  number  of  partial  reflections  in  the  cavity.  Referring  to  Fig. 1, 
the  path  length  difference  between  two  successive  reflections  E^  and  E2  is 
given  by  A ■ 2d  cos  6 where  d is  the  distance  between  the  plates  and  6 is  the 
angle  of  incidence  of  the  plane  wave.  This  path  difference  introduces  a 
correponding  phase  difference 


6 - kA 


4rrd  cos  9 
1 


(1) 


241 


where  X is  the  wavelength  of  light  between  the  plates. 

If  the  complex  amplitude  of  the  incident  wave  is  taken  as  Eq,  then  the 
amplitudes  of  the  transmitted  waves  are 


E = E t 
1 o 


„ „ J*  2 i6 

E = E t r e 
z o 


„ 24  2i6 

E = E t r e , etc. 
3 o * 


where  r and  t are  the  amplitude  reflection  and  transmission  coefficients.  A 
constant  phase  factor  corresponding  to  a single  pass  through  the  plates  and 
common  to  all  terms  has  been  left  out.  If  the  transmitted  waves  are  brought 
together,  as  by  a focussing  lens  their  amplitudes  will  add,  yielding 

2 2 i6  4 2i6 

Et  = EQt  (1  + r e + r e + ...  ) (2) 

The  term  within  the  parentheses  is  a geometric  series  and  can  be  evaluated, 
thus  we  get 


. 2 i6  , _ i6 

1 - r e 1 - Re 


2 2 

where  R=r  and  T=t  are  the  fraction  of  intensity  reflected  and  transmitted 
at  each  interface.  The  quantity  of  interest  is  the  fraction  of  intensity 
transmitted  through  the  Fabry-Perot  and  is  given  by 

- £ (4) 

o — E E*  1 + R-2R  cos  6 
2 o o 

2 

With  the  substitutions  T = 1 - R (no  mirror  losses)  and  cos  6 = 1-2  sin  6/2 
we  finally  arrive  at  the  Fabry  Perot  transmission  equation 


4R  sin  (6/2) 
(1  - R)  2 


This  is  the  well-known  Airy  function.  Its  plot  as  a function  of  6 for 
various  values  of  R is  shown  in  Fig. 2.  Equation  (5)  is  periodic  in  its 
argument  6 and  becomes  a maximum  for  6 = 2nn,  m = 0, 1, 2,  ...  . Equating  this 
expression  with  the  definition  of  6 from  Eq.  (1),  we  get  the  maxima  condition 
for  Fabry- Perot  transmission 

2d  cos  0 = m\  m=0, 1,  2,  ....  (6) 

There  are  three  quantities  of  interest  referring  to  the  Airy  function. 

One  quantity,  the  contract  C,  is  the  ratio  of  transmission  maxima  to  trans- 

2 

mission  minima.  Transmission  maxima  occurs  for  sin  6/2  = 0,  and  transmission 
2 

minima  for  sin  6/2  = 1.  Substituting  into  Eq.  (5)  we  have 

£)  -i.  £).  - rrV- 

1 - R 


therefore 


C = 


VV,.,„ 

Wmin 


1 + 


4R 

1 - R 


(7) 


A plot  of  contrast  C versus  reflectivity  R is  shown  in  Fig. 3.  Notice  how 
rapidly  contrast  rises  with  R.  In  a real  instrument,  however,  there  are 
other  factors  wore  important  than  mirror  reflectivity  which  limit  the 
contrast,  such  as  deviations  in  true  surface  flatness  and  variations  in  gap 
width  over  the  aperture.  The  theoretical  contrast  of  our  instrument  with 
R = 0.985  is  17,512.  The  actual  contrast  available  is,  however,  closer 
to  1,000. 

The  second  quantity  of  interest  is  the  Free  Spectral  Range  (AXFSR) . 

Like  the  grating  spectrometer  the  Fabry  Perot  suffers  from  overlapping  of 
orders.  The  range  of  wavelengths  which  can  be  displayed  in  the  same  spectral 
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order  without  overlapping  is  called  the  Free  Spectral  Range.  To  derive 
this  assume  any  X and  X’  are  entering  the  Fabry-Perot  such  that  the  (m+l)*^ 
order  of  X'  exactly  overlaps  the  nth  order  of  X\  Clearly  then,  X'-X  pro- 
vides an  upper  bound  on  the  wavelength  shift  we  can  have  without  overlapping 
orders.  Since  the  maximum  condition  (cos  0 - 1)  is  met  for  both  X and  X'  at 
d (see  Fig.  4)  we  have 

2d  = (m+1 ) X ' (8a) 

2d  = mX  (8b) 

Solving  for  m in  Eq.  (8b)  and  substituting  in  (8a)  we  get 

AXFSR  = X - X'  = (9) 

The  third  quantity  of  interest  is  known  as  the  Finesse.  It  is  defined  as 
the  ratio  of  the  Free  Spectral  Range  (FSR)  to  the  full  width  at  half  height 
(FWHH)  of  the  transmission  function  of  the  Fabry  Perot.  In  the  ideal  case 
the  transmission  function  is  the  Airy  function  and  the  Finesse  is  a function 
of  the  mirror  reflectivity  alone.  This  is  known  as  the  Reflectiviry  finesse. 

The  FWHH  may  be  computed  by  setting  It/io  = 1/2  in  the  Airy  formula 
[Eq.(5)].  We  find  that 

sl„2  5/2  . 

or 

sin  6/2  -i  ^ 

There  are  two  values  of  6,  <$i,  and  62  which  satisfy  this  condition,  see  Fig.  5. 
Since  and  62  are  nearly  equal  to  2imr  we  can  use  the  small  angle  approxima- 
tion, yielding 

5i  - + 2niw’  62  = + 2m71‘ 
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Their  difference  6^  - 62  is  given  by 

. . 2(1  - R) 

«1-  «2  'yir 

From  the  definition  of  6,  Eq.(l)  and  referring  to  Fig.  5,  we  have 


(10) 


4nd  cos  9 
*1  = (A  + AA/2 


and  2 


4Trd  cos  8 
(A  - AA/2) 


thus  we  have 


6f  - 62  = |-  (A^d  cos  0)  ^1  - j ^ (4fid  cos  6)  + -|^ 


Since  ~ « A 


AA 


6l  - 62  = — r2  (4^d  cos  0) 


(11) 


Equating  this  expression  with  Eq . (10)  and  solving  for  AA  we  get  as  the  FWHH 

1 - R A2 


AA 


VTT  2nd  cos  0 


(12) 


The  Finesse  is  simply  the  Free  Spectral  Range  divided  by  this  quantity 

FSR  A7  VRZrrd  cos  0 


since  cos  0 = 1 


R -AA  2d  (1-  *)P~ 

F = *_yr 

R 1 - R 


(13) 


The  actual  instrument  finesse  is  more  than  a function  of  mirror  reflect- 
ivity. The  other  major  factors  limiting  the  net  finesse  are:  (1)  lack  of 
parallelism  and/or  planeness  of  mirrors,  (2)  the  size  of  the  pinhole  following 
the  interferometer  (determines  the  angular  spread  of  the  light  entering  the 
instrument)  and  (3)  diffraction  losses  arising  from  the  finite  aperture  of  the 
interferometer.  Each  of  these  factors  will  broaden  the  instrument  functions 
independently  and  the  net  instrument  function  will  be  a convolution  of  these 
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factors  with  the  Airy  function.  If  we  approximate  each  of  these  profiles 
near  a maximum  as  a Gaussian  curve,  we  may  obtain  a simple  expression  for 
the  net  finesse.  The  convolution  of  two  Gaussians  is  another  Gaussian  with 
a half  intensity  width  which  is  the  Pythagorian  sum  of  the  component  half- 
widths. Since  each  Finesse  is  inversely  Proportional  to  the  halfwidth,  we 
have 


N 


i 


where  Fp  is  the  net  instrument  finesse  and  the  F^  are  the  component  finesses. 

A brief  description  of  the  major  component  finesses  follows. 

Flatness  finesse:  F„  = M/2  for  X/M  plate 

r 

where  M is  the  fractional  wave  length  deviation  from  planeness  over  the  mirror 

aperture.  Normally  specified  for  X = 5461  A. 

Pinhole  finesse:  defined  from  d ^ | X/2  (15) 

cos  6 y Fp  ) 


therefore  Fp  = ^7^  , see  Fig. 
the  next  section. 

Diffractive  finesse: 


6.  This  will  be  treated  in  more  detail  in 


(on  axis) 


(16) 


where  D is  the  diameter  of  the  limiting  aperture.  This  factor  j s usually 
negligible  except  for  very  small  apertures  (less  than  1 mm) . 
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THE  FABRY- PEROT  AS  A SPECTORMETER 


There  are  two  ways  a Fabry- Perot  may  be  used  as  a spectrometer.  The 

first  way  which  may  be  called  the  static  or  photographic  method  fixes  the 

distance  d and  measures  X as  a function  of  e.  The  second  method  may  be 

termed  the  photoelectric  or  scanning  method.  In  this  method  e is  fixed 

usually  to  0°  and  1 is  measured  as  a function  of  d. 

The  photographic  method  uses  the  setup  shown  in  Fig. 7.  The  pattern  in 

plane  1 produced  by  an  extended  monochromatic  source  is  a series  of 

concentric  rings.  The  Airy  function,  Eq.  (5),  is  a maximum  for  2d  cos  ® = ml 

Now  is  axially  symmetric  and  the  lens  L2  converts  9 into  r by  r = f20. 

The  condition  for  a maximum,  cos  9 then  becomes  cos  -r  = • In  the 

Za  1 2 Z.O. 

small  angle  approximation  this  becomes 


r _ mA 
2f 2~  2d 


(17) 


This  is  the  basic  equation  for  the  photographic  method.  If  we  wish  to  plot 

the  Airy  function  as  a function  of  r (i.e.,  the  location  of  the  rings)  we 

must  have  a starting  point.  Let's  say  there  is  a central  maximum  for  the 

incident  wavelenght  -y  of  the  order  m so  ?d  = mQl^.  WH1  subscript  the  pth 

A ° 

ring  from  the  center  as  r„.  and  let  m = m 4-  p.  In  this  way  we  can  solve  for  the 

PA  o 


radius  of  the  pth  ring. 


2 - ,.2„  (*°  ~ PHA  , 

PA  2f2[1  2d  1 


since  m X.  * 2d,  we  get 
o A 


f2 


PA 


(18) 


This  plot  is  shown  in  Fig. 8. 
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Unfortunately  a central  maximum  in  a photograph  of  the  spectra  is  not  always 
the  case,  more  often  the  case  is  that  the  first  maximum  is  a ring  of  a 
certain  radius.  This  more  general  case  is  treated  as  follows.  The  Fabry- 
Perot  maxima  condition,  Eq.  (6),  can  be  written 


cos  I 


m 

2d 


(19) 


For  m and  d held  constant,  the  quantity  cos  0/X  is  a constant.  Thus  if  two 

wavelenghts  X and  X are  entering  the  Fabry-Perot  they  will  have  maxima 
A D 

of  the  same  order  m = m at  the  angles  0 and  0 where 

o AU  BO 


cos  0 


AO 


COS  0 


BO 


a c (2Q) 

Here,  the  first  subscript  refers  to  the  wavelength  and  the  second  refers  to 

the  relative  order  (0  = m , 1 = m +1,  etc.).  Since  we  are  dealing  with 

o o 

0 **' 

small  angles  we  may  use  the  approximations  cos  0=1-—  and  0 = r/f. 
Condition  (20)  then  becomes 


1 - 


AO 

2f2 


1 - 


2f 


BO 

2 


or 


since  — « 1 
2f 


1 E0 

\(^4> 


i - 


AO 
2f  2 

2 
AC 

2f' 


multiplying  and  ignoring  higher  order  terms 

2 2 


v(, , !»| . Iso) . x 

A'  2£  2£2/ 
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Finally 


P*  (rAO  " 4>  “ XB  " XA 


«iis  equation  allows  us  to  calculate  X if  the  other  parameters  are  known. 

The  effective  focal  length  of  the  system  f is  difficult  to  measure,  however, 
because  the  imaging  system  usually  contains  magnifying  optics  and  a camera 
in  addition  to  the  primary  lens  (see  Fig. 7). 


We  wish  to  eliminate  f from  Eq.  (21).  The  maxima  condition  (19)  may  be 

. . . cos  6 X cos  0 . 

rewritten  as  *=  — . For  fixed  X = X.  and  fixed  d,  is  a constant 

m la  A m 


for  all  6 and  integers  m.  Thus 


cos  eA1  cos  eA0 
m + 1 

o o 


Using  the  small  angle  approximation  as  before 


tn  +1 
o 


These  are  two  equations  in  two  unknowns  m and  f.  Using  the  left  side  of  (22) 

o 

we  have  „ r2  » / 


multiplication  yeilds 


"o(l  - t5)  " - ^i) 


VaI  VaO  rAO  1 

° 2f2  ° 2f  2f2 


Simplifying,  we  now  have  an  expression  for  mfi 


2 2 
rAO  ” rAl 


Substituting  this  in  the  right  side  of  Eq.  (22)  we  have 


/ 2 2 » A 

*rA0  “ rAL^  ” 2d 
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Ignoring 


because  it  is  smaller  than  1,  we  may  solve  for  f 


AO 


2f 


r2  d , 2 2 . 


AO  A1 


This  is  the  quantity  which  we  were  looking  for. 

Substituting  this  result  into  Eq.  (21)  we  get 


XX 

B A 

2d 


2 2 

AO  ~ rB0 
2 2 

rA0  rAl 


XB  ' XA 


which,  since  X X, 
B A 


2 2 

rA0  ~ rB0 
2 2 

rA0  ~ rAl , 


= X„  - XA 
B A 


but  _A  = AXFSR  so  finally 
2d 


XB  ' \ = (AAFSR) 


( 2 2 

' r - r 
AO  BO 


vrAO  ” rAl 


(24) 


(25) 


A graphical  interpretation  is  given  in  Fig.  9. 

Data  that  we  have  taken  from  Brillouin  scattering  measurements  on 
several  different  liquids  using  the  photographic  method  is  located  in  the 
data  section  of  this  report. 

The  photoelectric  or  scanning  method  of  Fabry-Perot  spectroscopy  obtains 
a spectra  by  varying  the  optical  path  length  through  the  plates  with  0 held 
constant.  The  condition  for  maximum  transmittance  through  the  Fabry-Perot, 
Eq.  (6),  can  be  written 

2nd  cos  0 = mX. 

The  index  of  refraction  n has  not  been  included  in  previous  discussions  be- 
cause it  was  considered  constant  and  nearly  equal  to  one.  ThaB  we  see 
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scanning  may  be  done  by  either  changing  n or  directly  changing  the  distance  d 
between  the  plates.  The  index  of  refraction  may  be  changed  by  enclosing  the 
Fabry-Perot  in  a pressure  tight  container  and  changing  the  pressure  of  the 
enclosed  gas.  Pressure  scanning,  as  this  method  is  called,  relies  on  the  fact 
that  over  a considerable  range  of  pressures,  the  refractive  index  of  a gas  is 
directly  proportional  to  the  pressure. 

The  distance  d can  be  directly  changed  by  employing  piezoelectric  trans- 
ducers. These  transducers  are  placed  symmetrically  behind  one  mirror  holder. 
Applying  voltage  to  these  transducers  in  parallel  serves  to  translate  that 
mirror,  thus  changing  d.  This  method  has  numerous  advantages  over  pressure 
scanning.  The  scanning  rate  can  be  much  faster  and  the  possibility  exists  of 
applying  feedback  voltages  to  the  transducers  to  correct  for  drift  thus  maxi- 
mizing the  instrument  finesse.  Since  our  instrument  employs  piezoelectric 
scanning  we  will  restrict  our  discussion  to  this  method. 

Unlike  the  photographic  method,  the  scanning  method  fixes  0 and  measures  X 
as  a function  of  d.  An  advantage  of  this  method  is  that  the  spectra  may  be 
displayed  on  a scope  or  chart  recorder,  thus  showing  the  line  shape  as  well  as 
the  line  position.  This  method  has  a additional  advantage  in  that  the  results 
of  each  scan  may  be  recorded  in  the  memory  of  a computer.  If  the  results  of 
each  scan  are  coherently  added  to  the  previous  and  accumulated  over  time  the 
signal-to-noise  ratio  of  the  output  spectra  will  improve  by  \/n  where  N is  the 
number  of  scans.  This  procedure  is  known  as  signal  averaging. 

Assuming  we  are  operating  within  the  FSR  in  the  output  plane  P of 
Fig.  10  different  wavelengths  will  produce  rings  of  different  radii  so  by 
cutting  an  annular  hole  in  a screen  before  the  detector  only  a certain  band 
of  wavelengths  will  be  detected  such  that  the  maximum  condition 

2d  cos  0 = mX 
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1 


is  met,  where  0 is  the  angle  of  the  annulus.  A band  of  wave  lengths  will  be 

transmitted  because  the  annulus  has  a finite  width  which  is  equivalent  to  a 

spread  in  0.  From  Eq.  (6),  then  there  will  be  a corresponding  spread  in  X. 

If  a single  wavelength  is  the  imput  this  will  manifest  itself  as  a line 

2 

broadening.  According  to  Eq.  (25)  r is  directly  proportional  to  X. 

Figure  11  shows  the  intensity  profile  of  the  aperture  and  the  Fabry-Perot, 

2 

both  plotted  against  r . The  output  from  the  Fabry-Perot  has  an  underlying 
periodic  structure  dependent  on  X.  Thus  we  can  define  an  annular  screen 
finesse  in  a similar  way  to  the  reflectivity  finesse — the  ratio  of  the  dis- 
tance between  peaks  to  the  width  of  the  aperture  function. 

fa  - § ■ 4^4 

a - a 
o 1 


but  from  Eq.  (24) 


2 2 . *\ 

rA0  rAl  d 


d (a"  - af) 
o 1 

In  the  limiting  case  when  a^  = 0 the  annulus  becomes  a pinhole,  thus  we 
have  a pinhole  finesse 


P da2  dD2 

o 

Where  0 is  the  pinhole  diameter.  Since  the  net  finesse  is  given  by 


i_  + -i-  + _L 

F2  F2  F2 

R F P 
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it  is  apparent  that  we  want  the  pinhole  finesse  to  be  as  large  as  practicable, 
so  that  it  is  not  a major  source  of  line  broadening.  Normally  the  pinhole 
finesse  should  be  two  to  three  times  the  desired  operating  finesse. 

There  is  a tradeoff  however  between  the  size  of  the  pinhole  and  the 
transmitted  intensity.  Referring  to  Fig.  12  we  see  that  if  the  angular 
dispersion  of  the  light  is  greater  than  the  annular  extent  of  the  pinhole 
then  some  light  is  lost.  But  if  the  situation  is  like  that  shown  in  Fig.  12b 
the  light  is  not  lost  and  the  resultant  transmitted  intensity  is  much  greater. 
The  conclusion  here  is  that  the  angular  spread  of  the  input  beam  should  be 
matched  to  the  angular  extent  of  the  pinhole  for  the  most  efficient  detec- 
tion. This  cannot  always  be  accomplished  due  to  finite  size  of  object. 

In  our  particulat  case  we  will  have  an  operating  finesse  of  about  70 
if  possible,  which  dictates  a pinhole  finesse  of  140  to  210.  Referring  to 
Eq.  (28)  we  see  that  a pinhole  finesse  of  200  required  a 250  p pinhole.  This 
is  an  angular  (full  diam)  deviation  of  9.8  X 10  ^ r for  a 25.4  cm  fl  lens. 

The  pinholes  we  have  and  their  associated  pinhole  finesses  are  as  follows 
where  f = 25.4  cm: 


(cm) 

D (p) 

FP 

.05 

500 

52 

.02 

200 

323 

.01 

100 

1292 

.005 

50 

5168 

Therefore  the  200  P pinhole  is  the  most  appropriate  size.  Since  the  size 
of  our  image  is  about  400  u we  can  meet  the  angular  matching  condition  be 
using  a single  lens  of  focal  length  f'  = (2)  (25.4)  = 50.8  cm  as  shown  in 
Fig.  13. 
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Fig.l.  Fabry-Perot  Geometry  for  Derivation  of  Transmission  Equation 


Fig. 6.  Geometry  for  Pi 


Fig. 7.  Typical  Experimental  Setup  for  Photographic  Method 
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12.  Pinhole  Aperture  and  Light  Intensity  into  Pinhole  as 
Function  of  8 (8  = 
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